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HAUPTAUFSÄTZE 


Zur Theorie der kompressiblen Potentialströmungen I 


Neue Linearisierung der Grundgleichung 
der ebenen adiabatisch kompressiblen Potentialströmung. 


Von H. Behrbohm und M. Pinl (Messerschmitt-A.G.) in Augsburg. 


Das Variationsproblem der ebenen adiabatisch kompressiblen Potential- 
strömung ist nicht selbstadjungiert. Seine Behandlung sowie die des zu- 
gehörigen adjungierten Variationsproblems führt auf die beiden Euler- 
schen Differentialgleichungen E (p,q.r,s,UD=0 und E(y,9,r.5,D=0 für 
die Extremalflächen z und z, wobei z=z (x,y, als Geschwindigkeitspotential 
der Strömung zu deuten ist. Ist » (a, ß,y) eine geeignet verallgemeinerte 
Minkowskische Stützfunktion für z und ebenso © (a,ß,y) eine solche 
für 2, so ergeben sich für E=0 und E=0 die Linearisierungen 
1 er e4 En 
Daat w33 + TV ayy—=0 bzw. wa t+W33+ Ger = 0 


Die Linearisierung 
Verlauf der 


mit bekannten Koeffizienten C (a, p,y) und C (a, ß, y). 


von E=0 gelingt, obwohl die Gestalt dieser Gleichung im 
Untersuchung gar nicht ermittelt wird. 


$ 1. Einleitung. 
Die Integration der Differentialgleichung für das Geschwindigkeitspotential einer ebenen 
adiabatisch kompressiblen Strömung ist nicht nur für die weitere Entwicklung der Strömungs- 
theorie kompressibler Medien von entscheidender Bedeutung, sondern auch wegen ihrer Zu- 
sammenhänge mit der Theorie der hypergeometrischen Reihen einerseits und der der Minimal- 
flächen andererseits mathematisch von besonderem Interesse. Beide Zusammenhänge waren 
bereits A.S. Tschaplygin bekannt [1]') und wurden seither in der Literatur mehrfach weiter 
verfolgt [2], [3], [4]. [5], [6], 17]. [9], [10], [11], [12], [13]. Methodisch war man dabei stets 
bestrebt, die zunächst quasilineare Grundgleichung zweiter Ordnung für das Geschwindigkeits- 
potential durch Legendres Transformation in eine lineare Gleichung zweiter Ordnung über- 
zuführen, genau wie man diese Behandlung aus der Theorie der Minimalflächen her seit 
langem kennt. 
Bekanntlich existiert nun noch eine zweite Methode, die Differentialgleichung 
iI+P)r—2pgqs+i+pf)t=0 


1) Die Zahlen in eckigen Klammern beziehen sich auf das Schriftenverzeichnis am Schluß der Arbeit. 


(1,1) 
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der Minimalflächen zu linearisieren, die von der Legendreschen Transformation wesentlich 
verschieden ist und eine Bedingung für die Minkowskische Stützfunktion der Minimal- 
flächen darstellt. Dieses Linearisierungsverfahren geht von der Differentialgleichung der 
Minimalflächen in der Gestalt 


d p ® 4 = . 
dzyiFpFtE'iyvVirprg en 

aus und gewinnt die gesuchte lineare Gleichung in der Form: 
dw dw 0° Be LE NER Fe 


ET ET 


Entsprechend gehen wir in der vorliegenden Untersuchung von der Grundgleichung des Ge- 
schwindigkeitspotentials in der Form 

N) Ö G 

2 SR,.- 2 - —_— 2 Sn : z —, nn ERERÄTREEATIE (1,4) 
YAtulp+q) Y yAtulp®+qQ) 


aus und gewinnen die gesuchte lineare Gleichung in der Gestalt: 


%”o,%o 1 N) 


3 + PR Era Son ME BE a Er (15). 
ET + (u — 1)“ a de 


Offenbar geht (1,4) für A=u=1 in (1,2) und ebenso (1,5) in (1,3) über. Dieser Über- 
gang von Formeln unserer Untersuchung in solche, die man aus der Theorie der Minimal- 
flächen kennt, ist für die ganze Studie charakteristisch, wo immer die Koeffizienten 4, u auf- 
treten. Auf die geometrische Begründung dieser Erscheinung gehen wir in $ 8 genauer ein. 
Das Variationsproblem der Minimalflächen ist bekanntlich selbstadjungiert. Demgegenüber 
ist das zur ebenen adiabatisch kompressiblen Potentialströmung gehörige Variationsproblem 
nicht selbstadjungiert. Es erscheint bemerkenswert, daß es ohne Kenntnis der quasilinearen 
adjungierten Differentialgleichung gelingt, die zu (1,5) adjungierte durch Verwendung einer 
(verallgemeinerten) Minkowskischen Stützfunktion linearisierte Differentialgleichung auf- 
zustellen. 

Da die Lösungen des einen Variationsproblems durch Quadraturen zu ermitteln sind, 
wenn man «die des adjungierten kennt, und da die beiden Differentialgleichungen nicht gleich- 
zeitig ihren Typus wechseln, besteht die Möglichkeit, die eine zu verwenden, wenn die andere 
versagt und umgekehrt. Ferner arbeiten wir dauernd im Raum der kartesischen Koordi- 
naten x, y, z, dessen x, y-Ebene mit der Strömungsebene zusammenfällt. Damit wird es im 
allgemeinen leichter sein, den Rückweg in die Strömungsebene zu finden, als dies bei anderen 
Methoden bisweilen zutrifft. 


$ 2. Die quasilineare Grundgleichung. 


Es sei iw der Geschwindigkeitsvektor der ebenen adiabatischen Potentialströmung vom 
Betrag w. Für die Dichte o des strömenden Mediums besteht bei adiabatischer Zustands- 
änderung die Gleichung 

1 
x» — lw’\.-ı 
e=em)=o,l! 2 = ’ 
worin wir mit x den adiabatischen Exponenten?) und mit o, und ä; die dem Anfangs- 
zustand mw = w, entsprechenden Werte der Dichte o und der Schallgeschwindigkeit @ bezeichnen. 
Nach Voraussetzung ist der Geschwindigkeitsvektor w wirbelfrei und der Stromdichtevektor 


ke} 


> quellenfrei, die Strömung überdies reibungslos und stationär: 
0 


mn \n 


rotw—0, divs=Idivimtwgradt=0. ...... al). 
0% 0% 

Wir führen kartesische Raumkoordinaten «x, y, z ein, identifizieren die Stroemebene mit 
der x, y-Ebene und deuten das Geschwindigkeitspotential 2=2(x, y) als Fläche im Raum der 
Koordinaten x, y, 2; seine Existenz folgt aus (2,1). Zur Elimination der Dichte o aus div 8 
verwenden wir die aus der Bernoullischen Druckgleichung hervorgehende Relation [6], [8] 


. .. * - —% f7 
*2) Für strömende Luft gilt «=1,405, d.h. u = - -— = — 0,2025. 





YııM 


ch. 


941 
ch 


er 





YVılaa 


Z. augew. Math. Mech. 





Bd. 21 Nr.4 Aug. 1941 Behrbohm u. Pinl, Neue Linearisierung der Grundgleichung. 195 
1 a? 
5 grad (m?) = — — grad o 


und erhalten für das Geschwindigkeitspotential 2 mit Rücksicht auf (2,1) 


d2 dz 1 i 
w=grad 2 = 152. Fr —=/{p,gy, divm 5! grad(w)=0 . . . . 22). 
Zur Elimination der Schallgeschwindigkeit aus (2,2) verwenden wir die Definitionsgleichung 
der Schallgeschwindigkeit zusammen mit der Bernoullischen Druckgleichung und der adia- 
batischen Zustandsgleichung [6], [8] und erhalten: 
„—1 1—x 


A=d@=i+ulpf+g), A=-si tm, u=— 


Setzen wir diese Größen in die Gl. (2,2) ein, so entsteht die quasilineare Grundgleichung: 





0’ 2 0’ 2 0’ 2 


+ up + M- pr —2pgs ++ up +) -Elt=0, r=),;, s= 23). 








‚i=1; 
0x0 y dy? 





Für die Diskriminante der Differentialgleichung (2,3) gilt 
(P+P?—-@)=0 für (PP) Zä 
Entsprechend ist die Grundgleichung im Unterschallgebiet von elliptischem, im Überschall- 
gebiet von hyperbolischem und im Gesamtgebiet von gemischtem Typus. 
$ 3. Die homogene Funktion » (a, ß, y) 
Wir betrachten die Vektoren 
e=f{a,b,c}, a= {a, ß,y}, sau -D.... .. DB 


mit den Komponenten 


1 1 1 
a=—(2u—1)pA °*, b=—(2u—l1)gA °", -Ba®#r... rn A 
=/+1—-W(pP+Q), 
En a R_ 
a=—(2u—1)pA °*, B=-Pu-l)gA °, yalku-NA °’" ....683). 


Da a=a, b=ß, c=*=y, sind die Vektoren e und a linear unabhängig. Dagegen sind a und n 
proportional 
1 


a=—2u—1)A ®*n. 
Ferner gilt 


. 21 3 «£ 
=—-py, B=-ay, +), i+ u 5 =A=|(5 i) 


und wir erhalten also die Beziehung 


yatı-ı? 


ea) + Hr 1-0. 2.2... 8) 


Die geometrische Bedeutung dieser Relation sowie auch die der Vektoren e und a behandeln 
wir in $ 8. Denken wir uns y aus (3,4) als Funktion von a und f berechnet und als solche 
in (3,4) wieder eingetragen, so entsteht die Identität: 


„2 un 


& (a, 8,7 (a) = 57, pale (a +) + Ar] -1=0 daß} 9) 


Leiten wir nach a und ß ab (und bezeichnen Ableitungen durch angehängte Indizes), so entstelıt 


dy yru-n 2(u—1) zu 19y 
9. +9, 5 —=2 a  lames 2) +Ay2]+2 wast=0, 
ne Hau — ja t (2 u -jp«? u (+ PP) +Ay’]+2 ig 1%“ da 
„2u-n 2 (u — 1) yau-ı \ dy 
ie ==» u m y + En au ; % 0 . 
33 +9, ° = 2uPß 27 u 1% He 12 ) [u a +Pp)+iy? ?]+2 /. Zu 12°| op 
3) Das Zeichen = bedeutet hier und an parallelen Stellen des Textes „identisch gleich‘; die Klammern {...} 


geben die Argumente an, in welchen die Identität gilt. 
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j u : ET e er. 07 .. une N 
Aus diesen Gleichungen gewinnen wir für die Ableitungen da’ 38 die wichtigen Beziehungen 
< ‘ [ 
oy a 1 oy ) 1 was a? —+ P? Le 
= ——— — -=— pP — 2 C=4+(u I (3.5). 
da y C(a,ß,y) op y C(a,ß.y) y® 
Jetzt betrachten wir das Skalarprodukt 
ag=zas+ßy+tye=w(aP,Y)) - - -» -» » 2.0.0.3) 


Die so definierte Funktion & (a, ß,y) ist in a, ö,y homogen von erster Ordnung. Daher, sowie 
mit Rücksicht auf (3,4), besteht das lineare homogene System 


ale —- a )+Py-we)+tyr @-w)=0, a, P, y 


‘ 


nn ia dy_ „9% 
x [22] u u z V)=V, N a "2 ERIIER, } / am 
a da da =) “a By, ng (3,7). 
0y 


7 
y—@3 +37“ ee | 1, dB 


y für a=ß=0 beschränkt und verschwindet nicht. Daher gilt 
0, 0, y10,0) 
00) +0. 00.0, (0.0 +0 07 97 a 0.0) --0 
y +0, C(00,7(0,0)= Ä=+0, —- =;7 ‚1, 0, =y(0,0)==0. 
y (0,0) 4 (0,0, 7 (0,0) = 4 Ti r tler FE r( 


Demnach besitzt das System (3,7) in einer gewissen Umgebung der Stelle a=ß=0 die 
einzige Lösung: 


Nach (3,4) bleibt 


C=Oa: MO, 2=W,, d. h. x= grad ET nk + 


$ 4. Linearisierung der Grundgleichung mit Hilfe von w» (a, ß,y). 
Wie wir bereits in der Einleitung hervorgehoben hatten und wie man leicht nach- 
rechnet, kann die Grundgleichung (2,3) in der Gestalt 
h p 0 1 


FIT "T- Si, 
Vitule+e) 9 Vita +) 
geschrieben werden. Verwenden wir die Ausdrücke (3,2) und (3,3) für die Komponenten der 
Vektoren (3,1), so bekommen wir einfacher: 


Ei. 
0x 


EEE Re er ir 


Somit ist bdz=— ady ein vollständiges Differential. Zufolge 
db yda=d(bz—ay)—(bdx—ady) 
ist nun auch zdb - yda ein vollständiges Differential und es gilt äquivalent mit (4,1) 
WERE EIER u eu. Ma 

Hierbei ist darauf zu achten, daß die angeschriebenen partiellen Ableitungen nach a und b 
bzw. nach a und 5 nicht nur explizit, sondern durch ce bzw. y hindurch implizit auftreten. 

Nach (4,2) genügt also w(a, ß,y) wegen (3,8) der linearen Differentialgleichung zweiter 
Ordnung: 


PSPRGEN.  AEUBEN 0 4,3 
[77] [477 Day 3, = . . . . . . . . . . . ‚>. 
aa 33 da “, opß 3 | ) 
Da &,; homogen von nullter Ordnung ist, gilt 
i a 
away + Pay tyay=0, By = — — Way — — Wir 
67 " y nv 


/ / 
und wir erhalten mit Rücksicht auf (3,5) für ®& aus (4,3) die lineare Gleichung 





1 j 2 2 
Va 71 + 33 + C Vyy = 0; C =A — (u 1) . 2 . . . . . . (4,4). 


4 











Ihre geometrische Bedeutung behandeln wir in $ 8. Gemäß ihrer Herleitung ist diese 
Gleichung an das Bestehen der Relation (3,4) gebunden. Wir können uns aber leicht von 


dieser Bindung befreien. Sei nämlich «a, £, y ein Wertetripel, für das (3,4) nicht besteht, 


und sei etwa 


20 21 
2 “u ‚ 2 ı 


la ++ rg, 


En Bu-ım[tn =o=0.... (4). 





a? +) 





YUM 


b 


B 
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Dann genügt das Tripel 


were. wen: Ve SR 2 5 2 FE 7 
o } o } o 
der Relation 
ed ar? + p*® ’ 
a re et ee ii 


also der Bedingung (3,4). 
Ist nun ® (a*, 5*, y*) eine (von erster Ordnung homogene) Lösung von 


1 
Dart + @3r 3 + ar? + p*: BR er 
i+l -YD—a 


so ist, weil die linke Seite von (4,8) in «a*, 5*, y* von minus erster Ordnung homogen ist, 
wegen o==0 auch m (a, $,y) Lösung von 
1 
Daatws3 + 2 3m =0 . . .. er 
Artlm 1° = 
Damit ist die Gültigkeit von (4,4) von der Bindung (3,4) der Variablen «a, ß, y be- 
freit; (4,4) gilt im ganzen a, ,y-Raum mit etwaiger Ausnahme des Nullpunktes. Auf diese 
Weise ist das Hauptziel der Arbeit erreicht: durch Hinzunahme einer weiteren un- 
abhängigen Veränderlichen kann die in der Gestalt (23) zunächst quasi- 
lineare Grundgleichung in der Gestalt (44) linearisiert werden. Jede 
von erster Ordnung homogene Lösung »=o(a,ß,y) dieser Gl. (4,4) führt dureh 
Gradientenbildung 


z=@,(n,ß,y), y9=ws(a,ß,y),, 2=wr(a,ß,y) 


und anschließende Bindung von a, ß, y durch die Relation (34) zu einem 
zweiparametrigen Wertetripel«, „, z, aus dem nach Elimination der beiden 
Parameter eine Lösung der Grundgleichung (23) in der Gestalt z2=z(x,y) ge 
funden werden kann. 

Es ist nun möglich, nach Einschaltung von Quadraturen die Grundgleichung (2,3) noch 
auf eine andere Weise zu linearisieren. Das soll im folgenden in den Paragraphen 5, 6 und 7 
geschehen. 


$ 5. Die Grundgleichung als Integrabilitätsbedingung. 


Wir betrachten das partielle Differentialsystem erster Ordnung: 


ee Ken eb „. ..,. 2.5 MM 
zu =c bq, Yyza 4, Zy = —a en et re 

für die unbekannten Funktionen 
=), yeah ses) - - -» : .-.- GM 


Nach (3,2) und (3,3) gilt 





de __dc_ 1 o(By) 1 Ö |; _(l +P) 
da da 2u-l da 2u—l Az a y | 
1 | a?’ +? 0y a) 
—— ) ’ n 2 -9- 
g- ir 1 ee 2 (u ) „fh 
sodaß sich nach (3,5) 
dc a 
—— — con 
da y l 
ergibt. Analog zeigt man das Bestehen von 
dc FRE: P Er 
ob Y Kr 


Also ist 
dc dc dc dc 
Ze PL + 3b b,=pa, + gb, (y= da Wu + 3b b,=pay+ gb,. 


Hiernach und zufolge (4,1) sind die Integrabilitätsbedingungen 
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xy Lyx =—b, P—cx+ bx gq=a,p+ b, 4—-%x=0, 


Yu Yu AP Cy—Adg=ayp+byg—cy=0, 
2xy— 2yx >= b, +a, =0 


erfüllt und also mit dem Bestehen der Grundgleichung (2,3) äquivalent. 
Die Funktionen (5,3) erscheinen so durch Quadraturen bis auf additive Konstante be- 
stimmt. Ihre geometrische Bedeutung behandeln wir in $ 8. Wir berechnen die Funktional- 


IN) _ - . — — 5 5 
determinante TE „= #x9y—%yYx und erhalten aus (5,1) und (5,2) 
2 5 — bp, ap-—c 1-1 
(9) _ =c(e-ap-b99=BA "...... #5. 


(N) ec bq, ag 
Aus (5,3) eliminieren wir x, y und bekommen 
a N u 1 
2=2(%,4), pP=r=32 1-5-7J5- 
Zur Berechnung von p und g erhalten wir mit Rücksicht auf (5,1) und (5,2) das lineare System 


- >, + y =b=—bpp+tl(ap — og, 


zZ 
1 


x 


wi 
I 


27%, +7 y=-a=(—bgy)p+tagg, 


1 
und daraus (BA *"=0) 





- —ale—ap—bg) a - b(ie—ap—bg) b L- 

ad re © er rer (6,5). 
8$ 6. Die homogene Funktion 5 (a, B, 7). 

Wir betrachten die Vektoren 
=659), i=ladr, Kelia-D..:...:.:.06 
mit den Komponenten 

u (6,2) 
a=—pc, B=—ge, =, f (6,3). 


Da a=a,b=P,c=+y, sind die Vektoren € und ä linear unabhängig. Dagegen sind @ und ı 
proportional 


ee m 
= ring n. 
Ferner gilt nach (3,4) und (5,5) 
a 2u—l e b 2u—l 2 ee DB _ £ 
7 -7— a bzw.p=5 IP, 1=3- 778 - (6,4). 
In Analogie zu (3,6) bilden wir jetzt das Skalarprodukt ; 
ar=sar+ßy+Yz=nlaBY); I=SEHNN.: .:. 0.0.0. 165) 
Aus (6,2) und (6,3) erhalten wir für die Parameter a, 5, 7 die Darstellung 
1 1 > BR 


2u ” 2u 


.. F 4 > ze —] ai de 4 ıp 
a=a(p=—-pyA °, Bed, d=—gA ‚ 7=Y9)= Fr (6,6). 


Insbesondere gilt 
1 


2u 


2 (0,0) —3(0,0)=0, 7(,0)=._*- N 7.727 
2u—l 


Die nach (6,5) definierte Funktion ® (a, ß,y) ist in a, ß,y homogen von erster Ordnung. 
Berücksichtigen wir diese Eigenschaft und betrachten (6,5) zufolge (6,6) als Identität in p 
und 4, so ergibt sich in Analogie zu (3,7) das System 





XUM 





wVYııaa 
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da _ B I 
Op (x a) + AN p Bike. ‚@ o)—=0, | 
da JR REN m (6,8). 
dg 9)+- 9% zu, 3 2 —- )=V, | 
a@-)+ PY-)+ 7E-m)=0 
Aus (6,6) folgt: 
ee 1 1 
da 24 Ju < z Ry/] Er Ju \ 
OP. A [A — (u —1)p?] Fan — I)pgy AA (69). 
= 1 2. 1 
da 31 oß + ;, 
er MH — )pq A . da _— A “ [A (2u 1)q?] (6.10). 
Insbesondere gilt 
= 1 > ” u 
da 14 dB va Ay; 142. 
ze. ee‘ =0, —v, el 2 
dp oo pP ve 4 oo Ö 00 


Somit erhält die Determinante des homogenen Systems (6,8) an der Stelle p=q=0 den Wert: 


3a aB 3% 
3 


op’ op’ dp 45 
> _ .. & E- 
da op 0Y =-3 9. 
: : Zu—| 
öq’ dqg’ dq ’ 
a P) ß . y p=4= u 


Daher besitzt das System (6,8) in einer gewissen Umgebung der Stelle y=q=% die einzige 
Lösung 


z=o, y=@w;, 2=@ bzw. T=grad . (6,11). 
8 7. Das Analogon zu (4,4). 
Nach (5,4) besteht X, yy — £yYx=F0 und daher die Möglichkeit der Auflösung von (5,3) 
nach x, y in der Gestalt 
zZ=2(tYy), yzyay, z=er(ele,Yy),Yy(la,Yy)), vanylelke, m). ya, m) 
Aus diesen Identitäten folgt durch partielle Ableitung nach x bzw. y 
dad dr0y g_Ir9% dxd y 
 dxdxr dydx’ dzdy dydy' 
a aydz , oydy dydr , dydy 
ddr dydx’ dr0y 9ydy' 
oder mit Rücksicht auf (5,1) und (5,2 
öx dr 
1 = -bp,z t(ap uFr% 0 = .. ser | 
(7,1) 
0x 0x oy | 
)—(e — — — , =(—! _ 
0=(ce—bq) N =+ ag, z 1 =(e u) +a 15 - 
mit der Systemdeterminante 
1 
wr“ c 
(ap—ce)(bq—e)—abpgq=BA '=- er 
me 0x dx dy dy ni RO iR n 
Ö y —, — 4 rw B & : 1 (6, = 
Lösen wir (7,1) nach 32’ au’ I2' 3y auf und verwenden (5,5), (6,2), (6,3) und (6,4), so 
bekommen wir 
dx PT 5 02 0zdx, dzdy — ” 
 ___ — = — = u t . . . . . . 2 . 
0x bp, oz "PT 9% dxdz dydx x 1,5) 
BE: 0: Bin De Dee Bi - — 
—_—_c— = E— ‚B). 
y e—bg, 3409 %y = 3297 Iydy a (7,3) 
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Die Analogie dieser Relationen zu den = ® 1) und (5,2) ist evident. Sie gehen aus- 
einander hervor, wenn die Größen «, Y, 2, P, ‚b, ce bzw. &, 9, 2, p, 9, a, b, € mit Quer- 
strichen versehen werden und dabei die Eikecen Regel #—x, y=y, 3=3,P9=Pp, ... usw. 
beachtet wird. Den inneren Grund dieses Sachverhaltes behandeln wir in $ 8. 

Jetzt betrachten wir die Integrabilitätsbedingung 2, —=2,;. Sie bedeutet nach (7,2 
und (7,3) das Bestehen der Relation 

7 +b=0, 
woraus wir analog wie in $ 4 auf das Bestehen der Gleichungen 


rm u ie... ee 


schließen. Wiederum hat man bei den partiellen Ableitungen nach @ und 5b bzw. nach a 
und ? zu berücksichtigen, daß diese Argumente nicht nur explizit, sondern durch € bzw. 7 
hindurch auch implizit auftreten. Nach (7,4) und (6,11) genügt also die Funktion © (a, ß, y) 


in einer gewissen Umgebung der Stelle p=q=0 bzw. «= ß=0 der linearen Differential- 
gleichung zweiter Ordnung 


Tr RR 
at w33r+ ni, sr  ° 
14 
da "art 33 


. ei 0y Od, m ; "a 
Zur Berechnung der Koeffizienten = und ,„ differenzieren wir y nach » und 4 und erhalten 
on ( 


für diese Ableitungen aus (6,6), (6,9) und (6,10) 
37 


‘ 


o0yY0da 0y02 
Zu 


Dy Vy 2 u 
d» .n dBöp' d.h. [A- Zu —-1)p?] ze 2u—1l)py 33 =» B u 22 3 
und 
07 __Ody da, 0y092 
dq dadg Opdg’ 


Ar] 2(1 
d.h. —(2u 1) va ura euDalya=alß Au Dal. 


Daraus erhalten wir mit Hilfe von (6,2) und (6,3) 
rn 212 u 
O7 (=, 

da 


2;/> 
. 0 DE I SARU 
le: 2u—l, y 


‘ 


Da ©; homogen von nullter Ordnung ist, gilt 
a0; +B03,+YyWw757—0, 
Ferner genügt c der quadratischen Gleichung 
10 — (2u—V)ye+(- u) (@+ P°)= 
Daher verwandelt sich (7,5) in Analogie zu (4,3) in 


2/(2u 


ER 
Dat 933 — ra ri Day 


oder mit Rücksicht auf (6,2) und (6,3) i 





2 — Fri 09] De u ou = 15 
Vaat+ 033 Go — a zu )y+ylzu TI, 44(1— u) (a? + P°) (76). 











Damit ist das Analogon zu (4,4) gefunden. Auch hier erscheint die Gültigkeit dieser Gleichung 
zunächst durch die zwischen a, 3; y herrschende Verknüpfung eingeschränkt. Durch Elimi- 
nation von p und g aus (6,2) und (6,3) erhält man zunächst eine Beziehung zwischen a, 3, © 
Ersetzt man dann noch © durch die n (7,6) angegebene Funktion von a, ß, 7, so ergibt sich 
die erwähnte Beziehung zwischen a, 3, y (in selır komplizierter Form). Mit ihrer Hilfe kann 
man die in $ 4 in den Gl. (4,5) bis (4,9) verwendete Schlußweise benutzen, um auch die 
Gl. (7,6) von der speziellen Bindung der Variablen a, 5, y zu befreien. 

Mit (7,6) ist eine w eitere Linearisierung des Integrationsproblems der 
Grundgleichung (23,3) gewonnen, insofern als die Berechnung des Vektors x 
nurmehr Quadraturen erfordert, wenn (aus (7,6)) der Vektor £ bekannt ist. 
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88. Beziehungen zur Theorie adjungierter Variationsprobleme. 


Mit der Auffassung der Grundgleichung (2,3) als Eulersche Differentialgleichung des 
Variationsproblems 


2u—1 


ou ddrdy=0; deßtule +): .....2.. 68, 


die auf H. Bateman und G. Braun zurückgeht‘), war nicht nur die Möglichkeit der An- 
wendung invarianter Methoden auf unseren Problemkreis geboten, sondern auch der Versuch 
einer geometrischen Deutung des Integrationsproblems nahegelegt. Tatsächlich hat sich er- 
geben): die Integralflächen der Grundgleichung (2,3) sind als Extremalflächen des Variations- 
2u—1 
problems (8,1) die Relativminimaltlächen mit der Relativoberfläche: \\A+u(p +] °" dazdy 
und der Eichfläche e. Der Vektor n ist offensichtlich Normalvektor auf den Extremalflächen x. 
Da a zu n proportional ist, kann die homogene Funktion o=ar (aber nicht er!) als ver- 
allgemeinerte Minkowskische Stützfunktion gedeutet werden. Die Funktionen (5,3) sind 
die Komponenten des Ortsvektors { der Extremalflächen des zu (8,1) adjungierten Variations- 
problems. Der zu £ gehörige Normalvektor ist n, die zugehörige Eichfläche & Wiederum ist a 
proportional zu n und die homogene Funktion » (a,3,y) spielt die Rolle einer verallgemei- 
nerten Minkowskischen Stützfunktion der Flächen £. Deuten wir e und & als Relativ- 
normalen der Flächen r und x, so gilt: die Relativnormale e von x ist Normale zu r, die 
Relativnormale e von £ umgekehrt Normale zu x. Alle diese Beziehungen zwischen den beiden 
Variationsproblemen sind involutorisch. Für A= u«=1 wird das Variationsproblem (8,1) selbst- 
adjungiert und mit dem der gewöhnlichen Minimalflächen im Raum der Koordinaten x, y, 2 
identisch. Dann gilt 


A=-i+p+Q, B=1, C=1, C=1. 


Sämtliche Vektoren a, a, e, &, n, n werden proportional zueinander und die beiden lineari- 
sierten Grundgleichungen (4,4) und (7,6) werden miteinander und mit der Potentialgleichung 
in drei unabhängigen Veränderlichen identisch. Die Relativnormalen gehen in die gewöhn- 
lichen Normalen über, die Zuordnung paralleler Tangentialebenen von Eichfläche und Ex- 
tremalenfläche geht in die von Tangentialebenen der Extremalen und Tangentialebenen ihres 
sphärischen Bildes über. Der Variation der Relativoberflächen längs der Relativnormalen 
entspricht die Variation der gewöhnlichen Oberfläche längs der gewöhnlichen Normalen usw. 
Ferner gilt jetzt an Stelle von (3,4) 


a? +? E- y° = 1 
und ebenso an Stelle der im Text des $ 7 erwähnten komplizierteren Beziehung zwischen a, 3, y 
®+%+ ’=1. 


Die Skalarprodukte (3,6) und (6,5) bestimmen dann die Hesseschen Normalformen der 
Gleichungen der Tangentialebenen, und die homogenen Funktionen ® und ® reduzieren sich 
auf Minkowskische Stützfunktionen im gewöhnlichen Sinne, welche den Abstand der 
rn . k N a ’ dc dc 
Tangentialebenen vom Koordinatenursprung darstellen. Die Relationen a IE und 
dc 


die entsprechenden (im Text nicht angeführten) für Js sind ein Ausdruck der parallelen 


dc 
a’ob 
Zuordnung der entsprechenden Tangentialebenen von x und e (bzw. von I und ®). Die auf- 
einander abwickelbaren Extremalflächen x und r fallen im allgemeinen (auch im selbst- 
adjungierten Falle) nicht zusammen. Bezüglich ihrer weiteren Eigenschaften verweisen wir 
auf das Schriftenverzeichnis [14], [15], [16]. 


$ 9. Parabolische Ausartungen, singuläre Lösungen, Charakteristiken. 


Wir betrachten die zu den linearisierten Grundgleichungen (4,4) und (7,6) gehörigen 
„charakteristischen quadratischen Formen“: 
3 
9= N aryıy, 9=-L Ur Yıy 
i,k=1 i,k=1 
mit den Determinanten 


4) Vgl. [4], [5], [6], f17]. 
8) Vgl. [13], [14], [15], [16]. 





i y R er F . a 2 th. Mech. 
Behrbohm u. Pinl, Neue Linearisierung der Grundgleichung. u 1941 





vu 8 
0, 2 ; 
la;r|= Km ee  ° 


0, 


Im reellen Gebiet gilt: 


ei 1 
C =o0 u ri <oO. 
Das bedeutet: die Differentialgleichung (4,4) ist für reelle Parameterwerte a, f, y von ge- 
mischtem, die Differentialgleichung (7,6) für reelle Parameterwerte a, 3, y von hyperbolischem 
Typus. Insbesondere gilt 

a =0 für ” =0 bzw. —— ER =o0. 

o. Ay? +(u 1) (a? +P?) = i+ (u Y(p+g)— 
Der Koeffizient C verschwindet nicht, solange 


2, —1 
P+T+i, >0 bw.p+tg+ 2a es 


Der Koeffizient 5 verschwindet nicht für beschränkte p und q mit 


y#0 bzw A=+0 bzw. a=#+0. 

Nach Definition des Geschwindigkeitspotentials z ist p®+ q? mit dem Quadrat des Betrags w 

des Geschwindigkeitsvektors w identisch. Somit gilt: die Differentialgleichung (4,4) 

2% +(—-D)w} 
“+1 

ö s j . ; 2a —1 w; 2 1 —1)w? 

Geschwindigkeitsbereich ren ) wa <m<" re Z mo 


| 
; 2% In; \ 
schem Typus. Der parabolische Grenzfall w’= a te Im entspricht dem 





ist im Geschwindigkeitsbereich 0<w’< von elliptischem, im 


von hyperboli- 


physikalischen Grenzfall der Maximalgeschwindigkeit 


Vito * (a*—=kritische Schallgeschwindigkeit), 


Ü ma ax 


wo Druck, Dichte und örtliche Schallgeschwindigkeit d zum Wert Null ab- 
gesunken sind. 

Der Grenzfall identisch verschwindender Schallgeschwindigkeit ist von besonderem 
mathematischen Interesse. Aus A=0 folgt durch Differentiation nach x und y: 


Zupr+2ugs=0, Zups+2ugt=0, d.h. rt—-s#®’=0. 
Infolgedessen reduziert sich die Grundgleichung (2,3) auf n 
-pr— 2pqVYrt—-gQt=0 bzwp’r=gqt, 


d. h. auf eine Identität, da p r=qg’t=—pgqs. Die Torsen A=0 erscheinen mithin 
als singuläre Integrale der Grundgleichung (23). Mit Hilfe eines vollständigen 
Integrals der Differentialgleichung A=0 (z. B. einer geeigneten zweiparametrigen Ebenen- 
schar) ist ihre Darstellung integrallos leicht zu erhalten und hängt noch von einer willkür- 
lichen Funktion einer Veränderlichen ab. 

Der verhältnismäßig einfachen Form der Koeffizienten der Gl. (4,4) und (7,6) steht die 
im Vergleich zu (2,3) erhöhte Anzalıl der unabhängigen Veränderlichen gegenüber. Doch 
bieten sich damit auch wieder mehrere Wege, die Lösung des Problems zweiter Ordnung 
durch Methoden zu fördern. die nurmehr die Lösung von Problemen erster Ordnung voraus- 
setzen. Zunächst treten die physikalisch wichtigen Unstetigkeiten von Lösungen nur längs 
sogenannter charakteristischer Mannigfaltigkeiten auf, deren Gleichung erster Ordnung sich 
nach (9,1) in unserem Falle auf 

3 


ER] en Ser 


1 ge 


reduziert. Als geometrisches Hilfsmittel bietet sich jetzt die zugeordnete ternäre Riemann- 
sche Metrik 
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3 
LTE TIERRT N +dR®+C(aß,y)dy, Nardsdn=d®+d®+laAnd? 


,k=1 
dar, worin \alk\ die zu \a;, reziproke Matrix darstellt. Ihre Nullinien 


d®=da’+dß#+Cdy=0 bzw d#®—=da’+dß?’+ ldy 


sind die charakteristischen Strahlen („Bicharakteristiken“) der Gl. (4,4) bzw. (7,6). Anderer- 
seits vermittelt aber die gleiche Riemannsche Metrik auch den Zugang zur Konstruktion 


der Hadamardschen Grundlösungen 


u (a, Bud)... T._ u (a, B, 1; @°, f ß; y’ 
vr i Vf 


V= DER 


der Gl. (4,4) bzw. (7,6), deren wesentlicher Bestandteil, der geodätische Abstand 


r 


=yYT bzw. =VfF, 


der Differentialgleichung erster Ordnung 


ra t+3 +Ca,ß,y)ü =1, bzw. +3 +C@aBr)G =1 


genügt. Eine eingehendere Untersuchung der hier angeführten Probleme erster Ordnung soll 
einer späteren Arbeit vorbehalten bleiben ®). 
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Über die Stabilität einer Laminarströmung in einem geraden 
Rohr mit kreisförmigem Querschnitt. 


Von J. Pretsch in Göttingen. 


Die Methode der kleinen Schwingungen, die von W. Tollmien und 
H. Schlichting für die Untersuchung der Stabilität ebener laminarer 
Strömungen entwickelt worden ist, wird auf die rotationssymmetrische 
Laminarströmung durch ein Rohr übertragen. Mit ihrer Hilfe wird die 
Stabilität der Hagen - Poiseuille-Strömung gegenüber zweidimensionalen 
Störungen bewiesen. 


I. Aufgabenstellung. 


Um das Problem des Umschlags einer Laminarströmung in die turbulente Strömungs- 
form von der theoretischen Seite her zu klären, haben W. Tollmien, H. Schlichting und 
H. Görtler in einer Reihe erfolgreicher Arbeiten die Stabilität einer ebenen laminaren 
Grundströmung gegenüber kleinen zweidimensionalen Störungen untersucht. Im einzelnen 
wurden behandelt die stationäre Reibungsschichtströmung an der ebenen Platte [1,3,5]') und 
die stationäre Reibungsschichtströmung bei selır großen Reynoldsschen Zahlen längs 
ebener [6] und gekrümmter [7] Wände, ferner die instationäre Reibungsschichtströmung, die 
sich an einer Platte [2] bzw. an der Innenwand eines Hohlzylinders [4] ausbildet, wenn in 
ruhender Flüssigkeit die Platte plötzlich aus der Ruhe heraus mit gleichförmiger Geschwin- 
digkeit geschleppt wird bzw. der Zylinder gedreht wird. In einer kürzlich erschienenen Arbeit 
hat H. Görtler [8] diese Stabilitätsuntersuchungen auf dreidimensionale Störungen an kon- 
kaven Wänden ausgedehnt. 

In der vorliegenden Mitteilung soll die Methode der kleinen Schwingungen auf die 
Laminarströmung in einem geraden Rohr mit kreisförmigem Querschnitt übertragen werden, 
bei der die Hauptströmungsrichtung parallel zur Rohrachse verläuft. 

Als Beispiel einer solehen Rohrströmung werden wir in dieser Mitteilung nur die be- 
kannte Hagen-Poiseuille-Strömung behandeln und ihre Stabilität gegenüber kleinen zwei- 
dimensionalen Störungen sehr einfach beweisen. Diese spezielle Strömungsform für sich allein 
ist schon früher von Th. Sex1 [9] untersucht worden, jedoch ist der Beweis ihrer Stabilität 
dort nicht zu Ende geführt worden. Für das ebene Problem der Strömung zwischen zwei 
ebenen Wänden wurde die Frage der Stabilität der parabelförmigen Geschwindigkeitsverteilung 
einerseits von H. Solberg [10] in einer leider noch nicht abgeschlossen vorliegenden Arbeit, 
andererseits von W.Me. F.Orr[11] auf der Grundlage der Energiedissipationsgleichung und 
von J. L. Synge [12] unter Zugrundelegung der Wirbeldissipationsgleichung behandelt, ohne 
daß sie jedoch bislang als entschieden angesehen werden könnte. 

In einer weiteren Mitteilung werden wir die von P. Szymanski [13] berechnete in- 
stationäre Rohranlaufströmung untersuchen, die sich ausbildet, wenn die Flüssigkeit im Rohr 
aus der Ruhe heraus unter dem Einfluß eines konstanten Druckgefälles beschleunigt wird, 
bis sie die parabelförmige Geschwindigkeitsverteilung angenommen hat. Diese Strömung er- 
hält man nur in Querschnitten, die genügend weit stromab vom Eintrittsquerschnitt liegen, 
so daß die Parabel wirklich den stationären Endzustand der Geschwindigkeitsverteilung 
darstellt. 

Des weiteren soll die stationäre Rohreinlaufströmung untersucht werden, die sich in 
der Einlaufstrecke vom Eintrittsquerschnitt bis zu genügend großen Entfernungen stromab 
entwickelt, in denen die Parabelverteilung erreicht wird. Diese stationäre Rohreinlaufströmung 
läßt sich in Analogie zu der von H. Schlichting [14] berechneten laminaren Kanaleinlauf- 
strömung zwischen zwei parallelen Wänden aus der Prandtlschen Grenzschichtgleichung 
ermitteln. 

Mit der Untersuchung der Stabilität der aufgeführten Rohrströmungen hoffen wir, einen 
theoretischen Beitrag zur Frage der Turbulenzentstehung im Rohr zu liefern. 


II. Die allgemeine Störungsdifferentialgleichung. 

Wir wollen die Differentialgleichung des rotationssymmetrischen Störungsproblems 
zunächst in voller Allgemeinheit aufstellen. 

Es bezeichne r den Abstand von der Rohrachse und & die Entfernung vom Eintritts- 
querschnitt, « die Geschwindigkeitskomponente in Richtung der Rohrachse, v die Geschwindig- 
keitskomponente normal zur Rohrachse und p den Druck. Dann lauten die hydrodynamischen 
Bewegungsgleichungen für eine rotationssymmetrische Strömung: 


1) Die Zahlen in eckigen Klammern verweisen auf das Schrifttumverzeichnis auf S. 211. 
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du du, du 10p, ee: Iou es | 


es TV ern TI 
of dx or 00x or” rOör 92 


Ov a. dv 10p RT 5. v  ÖO’v\ | 
— tu v = —— Fr + — — — — 
ot "dx or o or "DRr'rdr PR dee) 


und die Kontinuitätsgleichung hat die Form: 


o(ru) Or v)_y 
0x or 


Eliminieren wir in (1) den Druck, so erhalten wir die Bewegungsgleichung: 


du ev ( du 0° si v(* u "®v .. (? u 0 “ | 


otdr Atdr' \oxdr 9x or dx0r) r\ör dx 
O’ 10”u 1du O’ u 2.0) 
= »l[‘ 16-8 2 =: +2 er = 3 
or roör’ rodr oxdr dx 
Es sei nun U und V die Axial- und Radialkomponente der rotationssymmetrischen Grund- 
strömung, deren Stabilität gegenüber kleinen Schwingungen geprüft werden soll. Die über- 
lagerte zweidimensionale Störungsbewegung habe die Stromfunktion: 


ylar,t)=— plz, r) eis(«-et) (4), 


wo a die räumliche Kreisfrequenz und der Realteil von e=c,+ ie; ihre Phasengeschwindig- 
keit bedeutet. Die Störungsgeschwindigkeiten sind dann gegeben durch: 


10y_ gg 


Bi ie(x—ct) 
u = = e 
ror r 
1 N) y N 1 R Ö q 


— e'‘ (x — ct +iagp 
rox r oa / 


+ — — 
wenn ’ Differentiation nach r bedeutet. 
Zu der grundlegenden Bewegungsgleichung des Störungsproblems gelangen wir nun, 
. . Tr Y n. Di 
wenn wir die durch Überlagerung der Grundströmung und der Störungsbewegung hervor- 
gehende Strömung mit den Geschwindigkeitskomponenten: 


nr RE A EE 


in die Bewegungsgleichung (3) einführen. 

Unter Beschränkung auf die in g linearen Glieder und unter Weglassen der Terme 
von unterlegener Größenordnung [vgl. Anhang I] erhalten wir die allgemeine Störungsdifferential- 
gleichung zunächst in der Form: 


er oe” ji TEN ) er „)" 


vr 


39" 39 
F 4 . z 20" +20® e 


’ 
-a?g’+ 


r Ox 


2a’o\ 9U , 
)+ q 


, f ; U’\0y Udg’ og” 
A 2 2 j 
(2. ce—3a?l U + iz, rY sr L ox 





Für r>x gelıt diese Gleichung, wie zu erwarten ist, in die allgemeine Störungsdifferential- 
gleichung des ebenen Stabilitätsproblems über. 

Führen wir dimensionslose Koordinaten ein, indem wir die Geschwindigkeiten auf die 
Geschwindigkeit U, beziehen, die etwa in Rohrmitte bei Erreichen der parabelförmigen Ge- 
schwindigkeitsverteilung angenommen wird, die Entfernung x vom Eintrittsquerschnitt auf 
die Länge ! des durchströmten Rohres und den Abstand r von der Rohrachse auf eine Be- 
zugslänge der Reibungsschicht, etwa die Verdrängungsdicke in einem festen Querschnitt x, 
oder zu einer gegebenen Zeit #,: 

R 
. Ü (&o, fo, ?) 
5 |j1- tun), 
. m (2, , t,) 
0 
wo R den Rohrradius und U, die Geschwindigkeit in Rohrmitte bedeutet, so läßt sich (7) mit 


Be ya Bi a ee 
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Rey — led (9) 
auch in der Form schreiben: 
w-alr-£-#0)-{- 2) 
Ö). 
[ver er ee 2 
2 : 2 (20 -3a®U u +7)5? Zr + vr 








Diese lineare partielle Differentialgleichung 4. Ordnung stellt die allgemeine Bewegungs- 
gleichung der gestörten rotationssymmetrischen Strömung dar. 
Die drei Summanden auf der rechten Seite der Gl. (10) sind von derselben Größen- 


r 


ö, h B 7 
ordnung, da V, %: und z von der Größenordnung Re* sind. 


Setzt man das Geschwindigkeitsprofil der Rohrströmung U(r) bzw. die Amplitude der 
Störungsbewegung p(r) von x unabhängig voraus, so erhält man eine einfachere gewöhn- 
liche Differentialgleichung, in der der zweite bzw. dritte Summand auf der rechten Seite 
von (10) fehlen. 

Es ist nun sehr bemerkenswert, daß für die Stabilitätsuntersuchung genau wie beim 
ebenen Problem die beiden durch die x-Abhängigkeit der Grundströmung U und der Störungs- 
amplitude 9 bedingten Glieder in (10) bedeutungslos sind, wie wir im einzelnen in Anhang II 
und III auseinandersetzen werden. Für die Stabilitätsuntersuchung maßgebend ist also bei 
der rotationssymmetrischen Rohrströmung wie bei ebener Strömung nur die Form der ört- 
lichen Geschwindigkeitsverteilung, die ihrerseits zwar durch den Verlauf des Druckgradienten 
mittelbar bestimmt wird, jedoch ist hier wie dort sowohl die unmittelbare Wirkung des 
Druckgradienten wie der Einfluß der «-Abhängigkeit der Störungsamplitude vernachlässigbar 
gering. 

Da sich die partielle Differentialgleichung (10) auf eine gewöhnliche Differentialgleichung 
zurückführen läßt, empfiehlt es sich, die Geschwindigkeiten auf die örtliche Geschwindig- 
keit U„ in Rohrmitte und den Abstand von der Rohrachse auf die örtliche Verdrängungsdicke 





e U 
-\l1-7,)ar 
zu beziehen, so daß man mit 
Un 6" 
Re* — = (11) 
diese gewöhnliche Differentialgleichung schließlich in der Form erhält: 
T [7 ” P) m U’ 
(Ü - 9 + - a’) (1 -)9 
i 5) po” 3 op" 39 ‚ 2 a? r (12) 
er =) piN a r r- nt ich Ale —-+a'p 





mit den Randbedingungen, daß = und 2 an der Rohrwand (r= KR) verschwinden und in der 
Rohrachse (r =0) regulär bleiben. 
Die allgemeine Lösung der Differentialgleichung (12) setzen wir in der Form an 
„nn, rar Hi: ir use er 
und wenden uns nun der Berechnung der Partikulärlösungen zu, die derjenigen beim ebenen 
Problem sehr ähnlich verlaufen wird. 
II. Die reibungslosen Lösungen 1,,%:. 


Setzen wir im folgenden durchweg voraus, daß die Reibung der Flüssigkeit sehr klein, 
also «aRe* sehr groß ist, so vereinfacht sich die Störungsdifferentialgleichung (10) oder (12) 
weiter zu der reibungslosen Störungsgleichung: 
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„’ I 
( U e) (o” Ex - . a? r) “ (v” —) = 0. P n . . . . (14), 


deren Integrale @,, 9, wir zunächst angeben wollen. Auf ihren Zusammenhang mit den 
Lösungen der allgemeinen Störungsdiffentialgleichung (12) werden wir in Kap. IV eingehen. 
Die Differentialgleichung (14) hat zwei singuläre Stellen, nämlich einen Pol 1. Ordnung 
an der Stelle U=c, der wir den „kritischen“ Achsenabstand r=r, zuordnen, und einen 
Pol 1. Ordnung in der Rohrachse r=0. Wir werden aber sehen, daß nur die Stelle r=r, 
auch eine singuläre Stelle der Lösungen dieser Differentialgleichung ist, wogegen an der 
Stelle r—=0 die Lösungen y, und 9, regulär sind; die Stelle r=0 ist also nur eine schein- 
bar singuläre Stelle (apparent singularity) oder ein Nebenpunkt der Differentialgleichung (14). 
Wir führen nun als neue Veränderliche 


2 


Pen, = (15) 

ein, so daß U— c=0 für o=0 und U—c>0 für o >V wird. 

Die reibungslose Störungsgleichung geht dann über in die Form: 
3 d? p e d? U { 
(U-.e) (1 - 075 *% op 1077 Y=V EI EEE A 
wo 
r 

a, — — € (17) 


gesetzt ist. Die singuläre Stelle r=r. wird transformiert in die Singularität o—=0 und der 
Nebenpunkt r=0 in den Nebenpunkt o=1. 
Nach allgemeinen Sätzen über lineare Differentialgleichungen können wir ein Fundamental- 
system von (16) in der Umgebung der Singularität o=0, r=r, angeben: 
9,=0 PB, (0) \ 


= B,()+ Ayılno| 
Bevor wir auf den analytischen Charakter der Lösungen @,, 9, näher eingehen, wollen wir 
ihre Gestalt für die spezielle Klasse von Geschwindigkeitsverteilungen: 


(18). 


U r en 
EN RE, Me ei [2 
Di) Bel... 
berechnen, durch die wir die Geschwindigkeitsprofile jeder der in Kap. I beschriebenen 
Rohrströmungen annähern werden. 

Die reibungslose Störungsgleichung nimmt dann die Form an: 








d’p 
1 — (1 — 0o)"] d-)j;—ta ol+tnr(a —-Di—-o'o=0. . . . . (MD). 
Schreiben wir sie in der Gestalt: 
e d? ; ! 
"tt > En el a nn Katie Are 
i=1 
so ist die erste Lösung 9, gegeben durch: 
Te 7 a a 
u Fed 


‚r—( 
wobei die e, aus den £, durch das Gleichungssystem: 
v—] 
.6+Vr + Run =b; sel... ..:.2.22:20.:. 


u 
gewonnen werden. 
Die Lösung Y, erhält man zu: 


” 1 ei z 
9,=—-9,\>3d(r)=1+ y b,0’—2e, E Be. ; . 2:5, 5 
1 z 
m du ” } 
u) e } 4 R R > 9 
p; —1+ y b,0"—2e, , (Inlo]—ia) füro<V. . . (8). 
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Die Koeffizienten der Reihenentwicklungen %,, 9, Sind für n=2,3,4,5 in Zahlentafel 1 
bis 4 zusammengestellt. Für n=1 erhält man die Lösungen einfacher unmittelbar aus der 
Störungsgleichung (14), die sich auf die Form reduziert: 


(we)(g” — — a p)=0 De a Air, 2 
Ihre Lösungen sind neben der singulären Lösung U=e: 
2ir,. 
BETT J, (ia r) 3 R A r N P ’ : . r . . (27), 
0, = EEE N: A N RE, |. © 


wo J, die Besselsche und N, die Neumannsche Funktion vom Index 1 bezeichnet. Die 


ö j 2 
Lösung 9, fällt dabei im allgemeinen Integral fort, weil in Rohrmitte @, wie 7 gegen un- 


endlich geht. 

Um über die Konvergenz der Reihenentwicklungen ,, 9, eine Aussage machen zu 
können, erinnern wir uns zunächst, daß neben o=0 noch o=1 ein singulärer Punkt der 
Differentialgleichung (20) ist. In o=1 sind nun aber die Lösungen (22), (24) für @,, 9; 
regulär, so daß der Punkt o=1 kein singulärer Punkt der Lösungen ist, wie bereits erwähnt 
wurde. Infolgedessen konvergieren die Reihenentwicklungen 9,, 9, über o=1 hinaus und 


? , ; ehe u. 
zwar bein=1 bis 0>», bei l<n<s5 bis o„=2sin =. 


Da in Rohrmitte o=1 ist, konvergieren die Lösungen ,, 9, dort also für n<5 stets, 
R: 


2 
x 
es 


wenn auch ihre Konvergenz mit wachsendem » schlechter wird. An der Rohrwand o=1 


re 
R vi+o, 
ausreichen wird, da die kritische Schicht immer nahe der Rohrwand liegt. 

In dem singulären Punkt o=0 ist p, regulär, wogegen 9, dort eine logarithmische 
Singularität besitzt. Welcher Zweig des Logarithmus beim Übergang von o>0 zu o<0 für 
die Darstellung eines bestimmten Integrals genommen werden muß, läßt sich nur durch Be- 
rücksichtigung der Zähigkeitswirkung in der kritischen Schicht entscheiden und liefert die 
sogenannte Übergangssubstitution, die durch das imaginäre Glied in (25) ausgedrückt wird. Die 
Differentialgleichung für diese reibungsbedingte Übergangssubstitution leiten wir in Anhang II 
ab; sie ist genau die gleiche, die W. Tollmien beim ebenen Störungsproblem erhielt, indem 
er U und x von x unabhängig annahm. Wie in Anhang II gezeigt wird, ändert sich diese 
Gleichung aber nicht, wenn man U und von der Bogenlänge abhängig voraussetzt, so daß 
wir die von W. Tollmien angestellten Überlegungen nicht von neuem durchzuführen brauchen, 
sondern seine Ergebnisse ohne weiteres für das rotationssymmetrische Problem übernehmen 
können. 

Da das imaginäre Glied in (25) bei der Stabilitätsuntersuchung eine ausschlaggebende 
Rolle spielt, empfiehlt es sich, bei seiner Berechnung statt des Wertes 2e,, der aus der 
Näherungsverteilung folgt, durch die das Geschwindigkeitsprofil ersetzt wird, einen genaueren 
Wert zu nehmen, den man durch unmittelbare Differentiation der vorgegebenen Geschwindig- 
keitsverteilung erhält. Entwickelt man U(o) in der Nähe des kritischen Punktes o==0 in 
eine Taylorreihe, so erhält man hierfür die Beziehung: 


dagegen konvergieren die Lösungen nur für > ‚ was indessen für unsere Zwecke 











[d? U 
do? r, |[U” 1 
Gi S BR ha! An IC 
& e, = d U = > U’ ur or D . . . . . . . . . . (29). 
2 d 0 k=V0 





IV. Die Reibungslösung g,. 


Um die allgemeine Lösung (13) der allgemeinen Störungsdifferentialgleichung (10) dar- 
zustellen, brauchen wir noch zwei weitere Partikularlösungen %, und ,. Wir gewinnen sie 
aus einer Diskussion der allgemeinen Störungsgleichung, indem wir die reibungsbedingten 
Glieder jetzt mitberücksichtigen. 

Führen wir die neue Veränderliche ein: 


TE en ee 
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mit 


a Re (3) Kai +0 3 
= ar.Re(,,), ; ee 


S 





so erhalten wir, wie in Anhang II und III ausführlich abgeleitet wird, aus (10) unter Ver- 
nachlässigung der in e linearen Glieder die Differentialgleichung: 


r d! nn) d? N 

real Br DE Re ner nn ta Er 
Das ist genau dieselbe Differentialgleichung, die W. Tollmien beim ebenen Störungsproblem 
erhielt. Wie die Differentialgleichung für die Übergangssubstitution ändert auch diese Differential- 
gleichung sich nicht, wenn man voraussetzt, daß U und @ nicht nur vom Wandabstand, 
sondern auch von der Bogenlänge «= abhängen. 

Sorit ist auch beim rotationssymmetrischen Störungsproblem in allen Teilen nach- 
gewiesen, daß die x-Abhängigkeit der laminaren Grundströmung wie der Störungsamplitude 
bei der Stabilitätsuntersuchung nach der Methode der kleinen Schwingungen keine Rolle spielt. 

Die Lösung Y, von (32), die von der speziellen Form der zu untersuchenden Geschwin- 
digkeitsverteilung U (r) ganz unabhängig ist, können wir ohne weiteres vom ebenen Problem 
übernehmen, wobei wir auf die eingehende Darstellung von H. Schlichting [5] zurück- 
greifen wollen. 

Die Lösung 9, dagegen kommt wie beim ebenen Problem deswegen nicht in Betracht, 
weil sie mit wachsendem n, d. h. für kleines e, über alle Grenzen wächst. Wir haben also 
in dem Ansatz (13) für das allgemeine Integral C,—0 zu setzen. 

Über den unmittelbaren Zusammenhang und den physikalischen Charakter der beiden 
Lösungspaare 9,, 9, und ,, 9, gelten sinngemäß die Ergebnisse, die W. Tollmien [1] 
durch eine ausführliche Diskussion der allgemeinen Störungsgleichung des ebenen Problems 
in Wandnähe erhielt. Wir wollen deshalb auf eine Wiederholung dieser Erörterung ver- 
zichten. Die Lösung , braucht nur in unmittelbarer Wandnähe berücksichtigt zu werden, 
da sie nach außen rasch abklingt; sie stellt eine Art „äußerer“ Reibungsschichtlösung dar. 
Im ganzen übrigen Bereich der Strömung kann man das allgemeine Integral 9 aus den reibungs- 
losen Lösungen ,, 9, allein aufbauen, wobei in der kritischen Schicht eine „innere“ Reibungs- 
schichtlösung %,, hinzuzunehmen ist. 


V. Die Randbedingungen. 

Nachdem wir die Partikularlösungen der allgemeinen Störungsdifferentialgleichung ge- 
funden haben, stellen wir nunmehr die Randbedingungen auf. 

Da die Flüssigkeit an der Wand (Index w) haftet, verschwinden dort die Störungs- 
geschwindigkeiten u*, v*. Damit lauten die ersten beiden Randbedingungen, da die Abhängig- 
keit der Störungsamplitude @ von der Bogenlänge x sich als vernachlässigbar erwiesen hat, 
nach (5): 

rt inet -. ir. eerte h 


DEE. 0 in ae 


In Rohrmitte (Index m) muß rn endlich bleiben; da @, dort bereits abgeklungen ist, lautet 
die dritte Randbedingung: 

DEE: ea par 
Aus dem Gleichungssystem (33), (34), (35) erhalten wir die Gleichung unseres Eigenwert- 
problems: 

Yw Prw Pıw 

a a en era a ea 

Yım p m 0 


Diese komplexe Gleichung ist mit zwei reellen Gleichungen äquivalent, in denen die Ver- 
änderlichen a, Re*, c, und ce; enthalten sind. Beschränken wir uns auf den Fall indifferenter 
Störungen (c;—0), so erhält man aus ihnen nach Elimination von c, eine Beziehung zwischen der 


2; NEE 
Reynoldsschen Zahl Re* und der Wellenlänge A= r der Störung. Das ist die Gleichung 


der Indifferenzkurve, durch die die stabilen und instabilen Störungszustände getrennt werden, 
und deren niedrigste Reynoldssche Zahl die sogenannte Stabilitätsgrenze angibt, bei deren 


Überschreiten der Umschlag von der laminaren in die turbulente Strömungsform beginnt. 
14 
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Die Auflösung der Determinantengleichung (36) liefert die Beziehung: 


n N BR ’ 
Ysw _PımYPawT PamYıw 


= ’ Br, (37). 
Pw PımYfıw—PemYıw 
Ersetzen wir die Ableitungen nach r durch solche nach o, so wird aus (37): 
( "| (@ 2) 
( . . 
Y%,w __ dem" do 5 dla ’ 
De u wer "Ta en won . (38). 
des) d e.\ Pe ”.) 4 P5 ee) 
do vw \do/m\do/w \do Ar w 
Den Ausdruck 
l sw Yw a 
<; gr ”) — — dr) a A 5 
d 0 /u \d 7 /w 


hat H. Schlichting [5] in einer bequem auswertbaren Form angegeben und tabuliert. 
Führen wir seine Bezeichnungen ein: 


; D (N) 


u a 
Nw 
und 
( q ) (“ er) 
1 do u do „Fım } 
ee (dm) (des) (de) (dr) isst re 
‚do /m do ), = m do w 
wo 
ER RUE. 42 
a. vi 


ist, so erhalten wir die Gleichung der Indifferenzkurve schließlich in der Form: 
<= en ein tie (43). 

Zur Ermittlung dieser Indifferenzkurve trägt man den Imaginärteil von F(n.) über dem 
Realteil von F (n..) auf und zeichnet in dieses Polardiagramm zu vorgegebenem c, die Funktion E (a) 
ein. Die Schnittpunkte (43) geben zunächst das zu einem bestimmten « gehörende 7. und 
dann mit (30), (31) die gesuchte Kurve «a (Re*). 

Die Berechnung solcher Indifferenzkurven für die Geschwindigkeitsverteilungen der 
stationären Rohreinlaufströmung und der instationären Rohranlaufströmung bleibt einer 
weiteren Mitteilung vorbehalten. 


VI. Die Stabilität der Hagen-Poiseuille-Strömung. 


Um eine einfache Anwendung der in den Kap. II bis V entwickelten Stabilitätstheorie 
der Rohrströmung zu bringen, wollen wir hier auf sehr einfache Weise die Stabilität der 
parabelförmigen Geschwindigkeitsverteilung beweisen. 

Beim rotationssymmetrischen Stabilitätsproblem spielt die Veränderliche o [Gl. (15)] die- 
jenige Rolle, die der Wandabstand (y— y.) beim ebenen Problem innehat. Man vergleiche 
dazu die reibungslosen Störungsgleichungen beider Probleme und deren Lösungen und außer- 
dem die Definition von n. 

Beim ebenen Problem ist jede lineare Geschwindigkeitsverteilung U=cy stabil, weil 
U 
dy’ 
kurve umschlossene Bereich der a, Re*-Ebene, der die instabilen Störungszustände umfaßt, 
zu Null zusammenschrumpft. Deshalb ist die ausgebildete Couetteströmung stets stabil. 

Beim rotationssymmetrischen Störungsproblem ist analog jede Geschwindigkeitsverteilung 


im ganzen Strömungsbereich verschwindet. Das bewirkt, daß der von der Indifferenz- 


wu, a ; : ; 
stabil, für die dog m ganzen Strömungsbereich verschwindet, und das ist gerade bei der 


parabelförmigen Verteilung der Fall. Denn für diese ist: 

ME. BR, 

Dee ee 
VII. Zusammenfassung. 

Die mathematische Methode zur Untersuchung der Stabilität einer Laminarströmung 
durch ein kreiszylindrisches Rohr erweist sich derjenigen für eine ebene Laminarströmung 
sehr ähnlich. Auch hier ist die Abhängigkeit der Grundströmung und der Störungsamplitude 
von der Bogenlänge (Entfernung vom Eintrittsquerschnitt) von vernachlässigbarem Einfluß 
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auf die Bewegungsgleichung. Die reibungsbedingten Lösungen der Störungsdifferentialgleichung 
sind mit denen des ebenen Problems bei passender Transformation der Veränderlichen iden- 
tisch, die reibungslosen Lösungen werden für symmetrische Geschwindigkeitsverteilungen der 
Form U=1-—-r"(n=1,2,3, 4,5) berechnet, durch die die Geschwindigkeitsprofile der in- 
stationären Rohranlaufströmung und der stationären Rohreinlaufströmung für deren Stabilitäts- 
untersuchung in einer weiteren Mitteilung angenähert werden sollen. 

Als einfaches Ergebnis der entwickelten Methode erhält man, daß eine Laminarströmung 
durch ein kreiszylindrisches Rohr stets dann stabil ist, wenn die zweite Ableitung der Ge- 
schwindigkeit nach dem Quadrat des Achsabstandes r? im ganzen Strömungsbereich ver- 
schwindet, was bei der Hagen-Poiseuille-Strömung U=1— r? zutrifft. 
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Zahlentafel 1 bis 4. Reihenkoeffizienten der Lösungen 9,, 9, für die Grundströmung 
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Zahlentafel?2. n=3. 
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b,=-6+20}; = z&; berg ygütrze; 
9.38 56 9 3 ,, 12 =» 
eautrzutrze ET 67 Ka ELLE 
19 _ 106027 ,_, 9601 154 , 
97 40 7 58500 
39 399141 , 2391 43889 
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Zahlentafel4 n=5. 


























a, =%; B,= Ita; B»=-4—-1ua}; 
16 16 " 
B=—4tal; B=y ta; BR=5 ta; 
’ 4 8 56 
R=g ta; =-5 tu; Bergrin; 
2 
o,=1; e=—?2 ,=2+y4; 
5 1 4 2 = -; 3 
= 1-74; Sztrgatrgn Bgm 
4 934 “ — Er; 8159 4 278 
Sp htrgmatz er 
37 12031 398 
got zo tr a top a; 
_ 181 ,, 23219 ., 851, 408 
ot Tr tr a 
14 31 112 > 
,=—-10+2ua; ,=14+5a; ,=—-7— 7 atr3W; 
46 , 451 2 2 1068 , Zu. 
75 135 at 595 a; =57 675 . Aa 70875 ° 754 
4 9266 8179 4916 
y, == —- sum 75 a sonsE; 
75 34075 165372 33075 
„3681 , daeserz | WO , 2 
s=75 926100 * 13891500 7 10418635 315 ° 
„38147 196566042 ., 6872333 „, 12061 
= 750° 3778300 250047000 7 56260575 © 7 893035 % 
Anhang |. 
Ableitung der allgemeinen Störungsdifferentialgleichung (10). 
Setzen wir abkürzend: 
RER EBEN u ae re a A 
so folgt aus den Gl. (5) und (6) für v und v zunächst: 
du -7 ‚(e"r — p') 
u.” z riece -“ 
Mr _ 99V iace 2 : ‚a 
dtdxr Dtdr r + ia‘ da 
“u _aU e®/, dgp\ ef, „_, 09” 
a Re r (iav | =) 
”r_aV ’p 99,00) 
dr? da? r + sie, x? u". . 
® e? 
—=u"-° er — arg" +2rgN) en N‘ 
0 v Bi e? 99’ 9; dp _ 2 ')- a -8i dp 2 
drdz ag (tie "ei won een e) 
0° u m e® IV r7 6 _Mm 3 s. " Bu. 
Jr =U a —3rg2"+6r°E"— 6rt pt) 
SCHE LU SFR 20) (side LE" gr 
ed dm ;(2 - + . r a YET, ”. ) 
dv DV , e:jd* op ‚309 
dei da 5: tiere tatg) 
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Setzt man (I,2) in Gl. (3) von Kap. Il ein, so erhält man unter Beschränkung auf die in g linearen 


Glieder zunächst: 


DU iace ir i V iace? (? PL 9a dg ie 
” > q — )— ä — a°g 
- 7 a Bun FIT dtde r \dr dx / 
u, U dg"\ U e* dp” OU’g'e? 
RI PRDR 14 jap de”) _ 
F dr r? 9 Tyr r e dr a2 r 
U Er 2 +3ia p_ 9.a2 04 iade + u Ko 
dr? r \d r? dr? or ‚ dr r 
Ve: U” e? (dg dV’ 
ru". Pz „Pr — 272 "| 9r 4 (8? N ia )—VT 
fr es Y a: 7 Yuan 52 
Ve? /d:g’ 0g V e* /d°g dp 
N: BEI ’)+ (3 19; 2, 
DET Anderen 7 Mc ltr 
oV’e?/d Ve? V e? /d? R) 
E (=? +iap)+ (pr — p)+ E 8 ) a GEBR PR 
dor r\dr ! pe a „2° 7 Ör 
U'V ”_r.„ © Fe Fu u ; \ 
u IE : EEE RT 2 
r "7er r: \dz2 ' IT redelar v 
| rm e® IV pr? 3 .6,„! 6 u. 6 N u” 
una ie ud +65" —6rr 9) + 
Be. a a IE . „UU’ 
nV -PETEN- ZIEHT 
e* N) g' 0?’ E e? . DR dtp” ’ 
+ (4c: 2 2.2’) — dia—— 2 — 20? ) 
BE. ar dr? r') r + dr? ‚ 
vr ep, 0° q 0° q °P | 
— - dia -— 6a? - — 4ia? ——.a? 
d.r° r (3 Fu or: dr? >= ? 


Diese Gleichung läßt sich sofort vereinfachen, da die laminare Grundströmung U,V die Navier- 
Stokesschen Gleichungen allein schon erfüllt. Man erhält so aus (1,3) die Gleichung: 


BEE ER ET FD. av’ _ 18V 
wol Lee) (non 0) 


+3 yde 4V dg 3U ak 
| © r dr i dr | 
+„[4Ie_40e"_ dp, , dr] 
Ir dz Or a2’ ' dr) 
i v E q m 3 q „ 39 ’ ' q ’ ] 
=— IpIV 2 + — — 2a? y" 2a? arg) 
a [ a r? r3 Wr r + 2 


elr 3 "ia mM 


2avy B6.1A 0.12 20 2 RATE "UV g' U’ y' 


BE; ra r?a a ra dra 


v[29 „eg Mp,. „ap 
| - +b6ba | 
dr? 


+i 


a Y Vv» 18V2d4 U’dg Ve 18V’ 

zilzaedt _ go _10V00 War Von, area 
il Ör rad2? ardzdr ardte ad adrdr 
HD” dg U o° [72 


+ * ‚dp Udy  Udg” 
aodr? Or ardz 


when 





a dr aodr 


Die weitere Vereinfachung gelingt nun wie beim ebenen Problem wieder dadurch, daß man die Größen- 
ordnung der einzelnen Glieder in (1,4) abschätzt. Übertragen wir die dort angestellten Dimensions- 
betrachtungen sinngemäß auf das rotationssymmetrische Problem, indem U,,! und ö5 jetzt die in 
Kap. II der vorliegenden Untersuchung angegebene Bedeutung haben, so erhalten wir die Größen- 
ordnung der einzelnen Glieder der Störungsgleichung (1,4), wenn N wieder die Größenordnung der 


Störungsamplitude bezeichnet: 





(1,3). 


(I, 4). 
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Iren 2 r U 
r ’ 4 > 7 I NL x * 
(l c)|p” — a? - - - Rez?, 
v| x In" 39" 30’ ee 1 
| „IV -g4 ı 2% >49 WR WR 7 Bun | -a”g ) 
K Tr ——— 2a Tr + a! 
a |’ r m ag ® r vr 
‚2aVgp 379" , 3V/o , Vgp” EEE 
ı| Veg’— | 
r Di: ar + ar: ' 7% ’ dr a) 
+il2ac?? + 07% _U”dg 3aud*4 _Ugp Udg") _ NU, Re’ 
| a2 'er3z a dx $z era: a dr 1? oc 
? v’ 18V la p de’ 4Vdg | SV 0°] 
P ee En —. 2 —-( c 
"\dr r d2 | dr a dr? 
(1.5) 
4 dm’ Be ’ NÜU (1,5). 
Bee Be a BE 
Ir dr dr Ör| I® 
iv[2 08° 0°" 0? p i 8? V 
ü - | u 3” 7_ +90 * er p’ 
3). a|lr dx? N dr? a dr? 
En: i | vo _ 1 oVdyg ‚0 , 8W’dg 4: ‚gl _ NU, 
'“i r 32 r 092092 d2? ' 82982 °' dr? Re, ' 
0° 4 NU, 


v a —# 
d2° I?Rei? 
ivdtop NU„ 
a dz! I12Re,’ 





Da Re* als groß vorausgesetzt wird, berücksichtigen wir nur diejenigen Glieder, die von der 
NU NU j a : i A 
Größenordnung RB Re‘, und —,* Re, sind, so daß sich die Störungsgleichung (1,4) vereinfacht zu 














l 
der Form der Gl. (7): 
/ R “ U’ 
(I uw: ey) |(L R q 
\ F 
T- > [273 [23 “ ’ ’ 
ıv | Pan’ 3 q 3q ’ en = | 
a | v —2 r - r? a + 2a? adg 
En ‚,39 39" 2atg\ U | 
} '7 | rr ii se PR 
Ei ’ En a dx ' 
i| U\dP UdBg’ dg''] 
2ac—3aU —- U” ) _ +U 
+ a E n: + r/)dr ra dr 
aus der man in dimensionsloser Schreibweise die Gl. (10) und für sehr große Re‘, die reibungslose 
Störungsgleichung (14) erhält. 
Anhang Il. 
1). 
Ableitung der Differentialgleichung für die Übergangssubstitution unter Zu- 
grundelegung der allgemeinen Störungsdifferentialgleichung. 

Um die Reibungskorrektur für g, in der Nähe der kritischen Stelle U=c, an der die reibungs- 
lose Lösung 9, eine logarithmische Singularität besitzt, zu berechnen, nähern wir in Analogie zu der 
beim ebenen Problem von Tollmien benutzten Methode die laminare Grundgeschwindigkeit dort durch 
eine Parabel an: 

: _[dU 
U—ıx TR (11,1). 
Führen wir die neuen Veränderlichen ein: 
l- r 1 2 
ä E fl ö .„/[pBL Is. are, r ar 
- > mit ER Sa Re‘, Er, aı Eu jan 1 ser (11, 2) 
und 
r [>=2, 
Tr 
so erhalten wir zunächst die Transformationsformeln: 
X 
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on 2 N, 2 
ar gr. yl-en; ar: er 
> n dE\doJe=0 , 2 dre 
Im. B EEE ni 7, BE TR 
dr 3 (2), 3er. dE 
Ö 
F 2 ud i 5603. Me 
a 2 , ( 00 /» I dre 
ordz 3er, } 1-0 er) t5gmV1-en de 
00 Ja 
a (a0) 
en _ 2 de \do) s dr. 
ar dz Ber >) 3erd dE 
d0/, 
mit deren Hilfe wir die in der Störungsgleichung auftretenden Differentialoperatoren berechnen können: 
d _9n d 
ör drdn 
d° 9) »° h 0°? n d 
dr? dor) dn?  dr?dn 


dn0°?n 


0° (7) d° l 
dor or? dn? 


or: \or/ an 


+3 


o* -e”) (ar dv el d° 


























Om dr on? '\dr or? dn? dr? dm? 
(I, 4). 
EL on » 
dor »E dron 
” ann ® On #* | nd 
dzdr DBrdrdnm | drdndE | dArdedn 
u ‚-() 0 , 92n70dn on nm 0? 
drdr? 3r) dr d9 "ar dr 9dn? re dn? 
[dn\? » vn o? d n d 
+ (ar) OmdE | dr?dndE or?:dr0n 
Die Geschwindigkeit V ergibt sich aus der Kontinuitätsgleichung durch Integration 
4 
x 1 EN 
en uch dr m. 
0 
4; BU 1+2en/®U\]dre Fre de , 1+end /dU\ r 
a ’IimT26 i NER . (11,5) 
ay(-—em acht 3 (Fe), de ayü-enlget 35 j:lar), 
‚1+ten+en: d (| 
6 ds \d0o”)o| 
Die allgemeine Störungsgleichung (10) lautet somit in den Koordinaten £, n vollständig, wenn wir noch 
ör C 
1 — ne setzen: 
| aU\ 1 „ „[U\]| .« Ü | 4 UN 
(5), u 2 (de), ler Tre r? we ler er 
i [16 ‚"y_ 32 > 9 __ 8a? , 
. | en — A ER 4 
2 Re; leir ‚(1 en)? Im nie (1—en) = ie zu (1—en) tert 
| DU\ | 2en »U\\dr. 1+end [BU 1+en ten d (RU 
3 | = ) + 1+ ? ), e+fe = Me 7 - Ar )+ + + en sl Re 
SB [\ “rs 3 do®/,f de Fr 2 dE\do) 6 d&\do®), 
8 a(1—en)dp | 2a’p] 
20 1—en)- = 
&3 a ‘ 2+ ETe „er re | 
id ke hu BU d en. 2 dr. fjdU ) 
a 4 4 Ei E; au u —g 6 
TTRE er dn FriErJ enaslae), re de de = Aa 1 Der : N 
ic ( er 1 ‚el 4 U 
Mn 2 2 —8 2 _ 22 AA 2 . 
ToRes, f u. > ten (z, + Dad >71) „1 DIR 
| d =) 
dp ‚,dp|n d&\do 2 4% 2 
DU und z 8 en) 


E a dn| 3 ( 
el 
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1 /dl 
+ ic (tens) + ge" (),) er | Wentieen n?dr. , 4(1 _ närlae).| 


a Rey, do? 'K P. Inte dE "3er en | 
90), 
d e) 
Be; d°gq 0% J10(1—en)dE 80) 82 ee 
# TR BER. di ES an . 
+ a 1 tee 38er en 3er “ 
00 /o 
® 
Multiplizieren wir diese Gleichung mit En und berücksichtigen wir wegen der voraussetzungs- 
(5) 
‚00/o 


gemäß großen Reynoldsschen Zahl Re} höchstens die in e linearen Glieder, so reduziert sich diese 
allgemeine Störungsdifferentialgleichung zu der einfacheren Form: 


u U 
1 en? ‚bei, 4 _ (ae) A 
4 — 
Ge" Do) 
— i(1—2en) tie 2477 ink 6} (2) Ge 


de 1 d (dU d (e®U Ceredre 
+refasts ala) +5 + 3 12 


ni-en +, 


N (11,7). 


On? 








Um die Reibungskorrektur in der kritischen Schicht zu erhalten, entwickeln wir wie beim ebenen Sta- 
bilitätsproblem die Lösung 9, nach Potenzen von e: 


en Re 2 re EEE EEE (11,8), 
die sich in einiger Entfernung von der kritischen Stelle U=c an die reibungslose Lösung anschließt. 


Setzen wir (II,8) in (II,7) ein, so wird mit 90% =1 und Multiplikation mit — - 








&.|_&9 
d* Qa u a d° Pa “ e 00°), Sr 00°), Bi 
d nt (i en)—i Im n1—en)- ‚eÜ) en?|= et (1+89,) 
. "\do Es 
£ Ya 2 1°) U \dre jde 1 d (BU (11,9). 
ar 7 + Nee)t tal, fa: ER \aet3 ieloo), 
1 d/L | ‚ Cere ] 
ur: er) 7 MM 
Vernachlässigt man auch hier wieder die Glieder mit &, so erhält man als Differentialgleichung für y,;: 
(Tl 
Ya . 0? Y9 (de: e=V 
— =— ri ee 
N i av) (II, 10) 
5 0e—V 


Das ist nun aber dieselbe Differentialgleichung, die Tollmien für die Reibungskorrektur im ebenen 
Störungsproblem unter der Annahme erhielt, daß U und von der Bogenlänge x unabhängig sind. Es 
ist also auch für die Berechnung der Übergangssubstitution des rotationssymmetrischen Problems diese 
Abhängigkeit von der Bogenlänge z zu vernachlässigen und die ebene Lösung ohne weiteres für das 
rotationssymmetrische Problem zu übernehmen. 


Anhang Ill. 


Berechnung der Reibungslösung 9, unter Zugrundelegung der allgemeinen 
Störungsdifferentialgleichung. 


Die Differentialgleichung für 9, erhält man aus Gl. (II,7) [siehe Anhang II] für &>0: 


, 0%, dp 
Im tr am 


Das ist wiederum dieselbe Gleichung, die Tollmien beim ebenen Problem für y, erhielt, indem er U 
und 9 als nur vom Wandabstand abhängig voraussetzte. Auch bei der Berechnung der Reibungslösung 9, 
und somit in allen Teilen des rotationssymmetrischen Störungsproblems ist also die z-Abhängigkeit der 
Grundströmung und der Störungsamplitude zu vernachlässigen. 286 
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Bestimmung der Verzerrungsgrößen eines räumlich 
gekrümmten Stabes mit Hilfe des Prinzips von Castigliano. 


Von K. Marguerre in Berlin-Adlershof. 


Das Verfahren, die Verzerrungs-Verschiebungsgleichungen eines elastischen 
Gebildes mit Hilfe des Prinzips der virtuellen Kräfte zu bestimmen, das 
E. Trefftz für die Bestimmung der Schalengleichungen benutzt hat, wird 
in der vorliegenden Arbeit auf den räumlich gekrümmten Stab angewandt. 
Das Ziel der Rechnung ist dabei außer der Gewinnung der betr. Gleichungen 
in möglichst übersichtlicher Form ($ 5), Sinn und Reichweite des „Prinzips“ 
herauszuarbeiten ($ 2, und außerdem zu untersuchen, ob und inwieweit die 
Ergebnisse für einen Ring „endlicher“ Dicke verfeinert werden müssen ($ 7). 


1. Einleitung. 


Drei Gleichungssysteme bilden den Ausgangspunkt der Elastizitätstheorie: Die rein 
mechanische Forderung des Gleichgewichts für die Spannungskomponenten 


00x OT ‚or 
x uxı x 


2 27]} ER 
B’ntn RA (1.1) 


die rein geometrischen Beziehungen zwischen den Verschiebungen und Verzerrungen 


ou 
— Se re a ea ee 7 Ve ee 


EL =I—_ ....: 
x dx 


und das Elastizitätsgesetz, durch das die Spannungen mit den Verzerrungen verknüpft werden: 
o=ofe),, 2 B. o=F#e...... EI ESENET EIERN © © 


Das Elastizitätsgesetz ist empirischer Natur und kann wesentlich in keiner anderen 
Form, als (1.3) sie andeutet, angeschrieben werden — die beiden anderen Gleichungssysteme 
lassen sich dagegen auch durch gewisse Energieprinzipe ersetzen. So kann an die Stelle der 
Gleichgewichtsbedingungen das Prinzip der virtuellen Verrückungen treten, mit dessen Hilfe 
man je nach der Ausgangsform für die Variation der Formänderungsenergie (6 /J,;=ode 
oder ö/J;=d(1/2 Ee?}) die Gleichgewichtsaussagen ausgedrückt durch die Spannungen oder 


6 „1, [9° ou . m y 
durch die Verschiebungen erhält (= +++ oder ar Und an Stelle der Verzerrungs- 
R 2 


Verschiebungsbeziehungen kann man das Prinzip der „virtuellen Kräfte“ (Prinzip von Casti- 
gliano) benutzen, mit dessen Hilfe man (wieder je nach der Ausgangsform für den Energie- 
ausdruck) entweder die sog. Verträglichkeitsgleichungen bekommt (das sind die Gleichungen 
zwischen den Verzerrungen, die sich aus (1.2) bei Elimination der Verschiebungen ergeben) 
oder die Verzerrungs-Verschiebungsgleichungen selbst. Beide Prinzipe, insbesondere das der virtu- 
ellen Verrückungen, haben in den Anwendungen vor allem als Ausgangspunkt für Nähe- 
rungsrechnungen ihre Nützlichkeit erwiesen; indessen sind sie ihrem eigentlichen Sinne 
nach den entsprechenden Differentialaussagen (1.1) bzw. (1.2) exakt gleichwertig. Man be- 
dient sich ihrer mit Vorteil, wenn für das besondere elastische Gebilde, das man zu unter- 
suchen wünscht, der Ausdruck für die Formänderungsenergie sich leichter beschaffen läßt, 
als eines der Gleichungssysteme selbst — die „Mühle der Variationsrechnung“ liefert dann 
auf rein formalem Weg die gesuchten Gleichungen. 


Für Stabilitätsprobleme und für Probleme großer Verformungen kann es z. B. von Vorteil 
sein, sich (nach Beschaffung der kinematischen Gleichung auf „anschaulichem“ Wege) die 
Gleichgewichtsbedingungen am verformten Element auf diese Weise herzuleiten; bei anderen 
Problemen, wie dem der beliebig gekrümmten Schale‘) oder dem des räumlich gekrümmten 
Ringes, sind es die kinematischen Beziehungen, die weniger leicht anschaulich gewinnbar 
sind; dort liefert uns — wenn wir die Gleichgewichtsaussagen haben — das Castigliano- 


sche Prinzip einen formalen Weg zu ihrer Gewinnung. Allerdings ist dieser Weg nur gangbar 
8 g 


bei Beschränkung auf „kleine“ Verschiebungen, weil das Castiglianosche Prinzip nur inner- 
halb der linearisierten Elastizitätstheorie sinnvoll ist. 


1) E. Trefftz: Ableitung der Schalenbiegungsgleichungen mit dem Castiglianoschen Prinzip. Z. angew- 


Math. Mech. Bd. 15 (1935), S. 101 bis 108. 
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2. Zur Differentialgeometrie der Raumkurve; Bezeichnungen. 
Die Raumkurve denken wir uns gegeben in Parameterdarstellung, z. B. in der Form 


zes vere IE  . .. 2: = ro... 
mit der Bogenlänge s als Parameter. Wir bestimmen den Einheitsvektor der Tangente 


“ t=2. +9, +i%, . . . (22) 
(Punkte bedeuten die Ableitung nach s) 





%,>0 und wählen in der Normalenebene irgend 
t zwei aufeinander senkrechte Einheits- 
vektoren 
u. 
k,>0 h 
t die mit Hauptnormale und Binormale n, b, 
der Kurve zusammenfallen können, aber 
keineswegs müssen?). Den Vektor der 
MALER) nz „Drehgeschwindigkeit“ o dieses Dreibeins 
Bild la. Kurve mit begleiten- Bild 1b. Zum Vorzeichender zerlegen wir nach den Achsen t, n,, n, selbst: 
dem Dreibein. Krümmungen k,, ka. 


o=ht+kın,+k,n,. . (23); 
dann sind k,, k, die Komponenten des Vektors tXt der Kurvenkrümmung (der Winkel- 
änderung der Tangente t): ; 

kh=n,tXxt=-—t-n,] 

hr Be a er 

k,=n,.tXt=t-n, | 
während Ah, die t-Komponente von o, eine Art Windung ist (Vektor der Winkeländerung einer 
Senkrechten zu t, z. B. der Normalen n,, gesehen in der t-Richtung): 

Bir 3 ie et 
Wegen 

t:n=0, t-:,=0, n,:n,=0 


läßt sich statt der Formeln (2.4) auch schreiben 


h=—ü,:.n, k, .t, k, ü,-t. ie Ber |, ° 
und mit Hilfe von (2.4) und (2.4’) bestätigt man weh die Badkiligen 
BEE Kur Ei - : + RN 


die wir in Abschnitt 8 brauchen werden. 


Weitere Bezeichnungen. An einem Element des Balkens (2.1) mögen Schnitt- 

kräfte und Schnittmomente 

N=-Nt+Q,1n,+09,n,, \ 

M=Mrt+M,,+M,n, | ' 
angreifen, ferner eine äußere Last 

P=mi+rp,n,. +2 - » - 5 
Unter Einwirkung dieser Kräfte erleidet der räumlich Eätente Stab For een und 
zwar eine Verschiebung 


(2.5) 


aaa .., es E er Be > 
eine Verdrehung 

ae Fe 5 
eine Verkrümmung 

ee SV . (2.9) 


und schließlich Längen- und Winkelverzerrungen &„, Y,, 7;: Von diesen 12 RER 
größen sind die 6 ersten Verschiebungsgrößen, und unabhängig voneinander vorgebbar. 
Die 6 anderen sind die Verzerrungsgrößen, deren physikalischen Sinn, zugleich mit 
ihrer Abhängigkeit von den Verschiebungen zu bestimmen der Zweck der vorliegenden Unter- 
suchung ist. Dabei werden von Interesse insbesondere diejenigen nee sein, die 
sich bei Einführung der in der Balkenbiegungslehre üblichen Annahme y,=y,=0 (Vernach- 
lässigung der von den Querkräften unmittelbar hervorgerufenen ac ergeben. 


2) Für die praktische Rechnung empfiehlt sich das natürliche Dreibein t, n, b im allgemeinen nicht, da die Be- 
stimmung der Binormalen die Kenntnis der dritten Ableitungen von (2.1) voraussetzt, während man bei Benutzung eines 
beliebigen Systems mit den zweiten Ableitungen auskommt. Der Vorteil des natürlichen Systems, daß eine der Krüm- 
mungskomponenten verschwindet, fällt dagegen nicht ins Gewicht. 

3) h ist nur dann die „Windung der Kurve“ im Sinne der Differentialgeometrie, wenn als Bezugsdreibein 
speziell das natürliche Dreibein t, n, b gewählt wird. 





220 Marguerre, Verzerrungsgrößen eines räumlich gekrümmten Stabes er ns 





3. Die Gleichgewichtsbedingungen am Balkenelement. 

Die Anwendung des Castiglianoschen Prinzips, dessen Sinn am Beispiel des ebenen 
Spannungszustandes im Abschnitt 4 erläutert werden soll, erfordert die Kenntnis der Gleich: 
gewichtsaussagen am Element. In vektorieller Form lassen 
sich diese nach Bild 2 für das Balkenelement ds unmittel- 
bar hinschreiben. Das Gleichgewicht der Kräfte fordert 


NE ee © 
das der Momente 
RHHXR=I. . .....0 * 


Setzt man (2.5, 6) und (2.4”) in (3.1) und (3.2) ein, so erhält man 
die Gleichgewichtsaussagen für die Koponenten, die nıan ge- ” 

hr - £ . er Fa > ce Bild 2. Zum Gleichgewicht des Balken- 
wöhnlich die Gleichungen von Kirchhoff-Clebsch nennt: siemantse ds. 


N-,Q, +k,Q,=-p, 

DR Re nl... 2222. 
Q,—kN+hQ, =-p, 

Mr—k,M, +k,M, = —0, 

M,—hM, +k,Mr=-+0,,. 

M,—k,Mm+hM, =—-9, | 


4 Das Prinzip der virtuellen Kräfte. 
Zwischen dem Prinzip der virtuellen Verrückungen und dem Prinzip der virtuellen 
Kräfte (wie wir dieses gewöhnlich nach Castigliano benannte Prinzip zweckmäßiger 
nennen wollen *)) besteht im Falle kleiner Verschiebungsableitungen ee eine weitgehende 


dx 


Entsprechung. Wir geben hier für das zweite Prinzip eine Herleitung, die diese Entsprechung 
deutlich macht (um Schrejbarbeit zu sparen, für den zweidimensionalen Spannungszustand). 
Das Prinzip der virtuellen Verrückungen leitet man her, indem man die Gleich- 
gewichtsaussagen 
00, , Or dr doy 
"3 Fam det 3y 
mit virtuellen Verrückungen n öv multipliziert, durch Teilintegration umfornt, und zeigt, 
daß die so entstehenden Terme als „Variationen“ der Formänderungsenergie II; oder des 
Potentials der äußeren Kräfte //, gedeutet werden können. Auf diese Weise erhält man die 
Aussage 


ne Fe 


ST=8I,+öI,=0, 


die den Gleichgewichtsaussagen (4.1) gleichwertig ist. — Entsprechend gehen wir, um das 
zweite Prinzip zu gewinnen, aus von den geometrischen Aussagen 
du dv du dv 
an By 7Taytaz 
multiplizieren beide Seiten mit virtuellen Spannungen 
Ö0,, doy, Ör, 


und integrieren über den uns interessierenden Bereich xy. Wenn wir auf der rechten Seite 
Teilintegration anwenden, ergibt sich 


(4.2), 


\ ler do, +2,80, + yör)dedy 


—=$(d0,cosa+örsina)uds+ Plörcosa+ do,sina)vds 2), 
do,  Qdr dör en) 
-\\& +- Sy )udedy- + Ey vdedy 


Bild 3. wobei die beiden ersten Integrale rechts zu erstrecken sind über den Rand 
(s. Bild 3), die beiden letzten über das Innere des Bereiches. Genügen die 


dy ds 


dx 


%) Den Ausdruck „Prinzip der virtuellen Kräfte“ benutzt auch Schleusner (Die Variationsprinzipien der 
Elastizitätstheorie, Stahlbau 1938, Heft 24). 
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Spannungen öo,... den Gleichgewichtsbedingungen 
am Rande 
ö0,cosa+törsina =ÖE, 
örcosa+ öo,sina=öH 
im Innern 
800, , ddr 
dx Yy 
N) Ö T N) ß) Oy Fe > 
dx do e. 


=—ÖX, 


so geht (4.2’) über in 
\ lex d0g +2,00, +ydr)dedy— \(WmöEZ+vöH)ds— \\WöX+röY)dedy=OV 


Das erste Integral ist die Variation ö C, der „inneren Ergänzungsarbeit“ 


(x oy 


C(o,d)=|\ \edo,+leydo,+\ydrldaedy, 
vu 0 I) 


das zweite und dritte Integral stellen zusammen die Variation der „äußeren Ergänzungsarbeit“ 
dar, d. h. die Arbeit öC, der äußeren Kräfte (der virtuellen Randkräfte öZ, ö Fl und der 
im Bereich angreifenden virtuellen Kräfte öX, öY) an den elastischen Verschiebungen «, v. 
(4.4) schreibt sich daher als „Prinzip“ kurz 


ec ra N A 


Bei der Deutung der Größe C; ersetzt man meistens mit Hilfe des Hookeschen Gesetzes die 
&, y durch die o, r°) — indes ist dies eine zum Wesen des Prinzips nicht mehr gehörige 
Operation; die Gl. (4.5') ist unabhängig von der Form des Elastizitätsgesetzes. Unter der Vor- 


. „du En s ’ 
aussetzung kleiner Je" ersetzt sie die Verschiebungs-Verzerrungsgleichungen, die man also 


umgekehrt aus ihr herleiten kann. 


Im Einzelfall erhält man die Verzerrungs-Verschiebungsgleichungen, indem man, von 
(4.5°) und den Gleichgewichtsbedingungen (4.4) ausgehend, «die Teilintegrationen sozusagen 
„rückwärts“ laufen läßt, bis man auf die (4.2) entsprechenden Beziehungen stößt. Dabei 
behandelt man die Minimalforderung (4.5’), die unter Beachtung der Nebenbedingung (4.4) 
erfüllt werden soll, zweckmäßig nach Lagrange®°). Schreibt man 


so lautet die Minimalforderung, die jetzt ohne Rücksicht auf die nach (4.4) zwischen den 
Spannungen bestehenden Bindungen erfüllt werden soll, 


öl; -cJ)+ Fi, @,dedy=0........ 0m 


Die physikalische Bedeutung der Lagrangeschen Faktoren ergibt sich bei der Teilintegration 
des Terms 24,6G,. Dabei treten Randintegrale auf, die sich (unter Beachtung von (4.3)) 
gegen das Randintegral in (4.5) wegheben müssen; daraus folgt in unserem Beispiel 


h,=u, L=v, 


und da man das Gebiet durch „Ränder“ beliebig auftrennen kann, so muß das auch im ganzen 
Innern des Bereiches gelten. Die Multiplikatoren der Gln. (4.6) sind also die elastischen 
Verschiebungen selbst; wenn man sich die Herleitung der Gln. (4.5”) bzw. (4.7) aus (4.2’) ver- 
gegenwärtigt, ist das nicht überraschend, 


5) In dieser Form (und mit d C„=®) tritt (4.5*) als Castiglianosche Formel in den Anwendungen gewöhnlich 
auf. — Auch in der zitierten Arbeit von Trefftz wird das Castiglianosche Prinzip in diesem engeren Sinne benutzt; 
das Verfahren hat den Nachteil, daß die zwei Schritte, die zu den Spannungs-Verschiebungsgleichungen führen (das 
Elastizitätsgesetz und die rein geometrischen Verzerrungs-Verschiebungsgleichungen) nicht deutlich getrennt erscheinen. 

6, S, E. Trefftz: a. a. O. Z. angew. Math. Mech. Bd. 15 (1935), S. 106. 
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5. Die Bestimmung der Verzerrungsgrölsen des räumlich gekrümmten Stabes. 
Kennzeichnend für das Verhalten des dünnen Stabes ist, daß die von den Querkräften 
herrührenden Verformungen y,, y, so geringfügig sind, daß man sie neben den Momenten- 
und Längskraftverformungen vernachlässigen kann. Die innere Ergänzungsarbeit öC; enthält 
daher nur vier Terme 
86; =\W®5Mr +x,68M, +2,05 M, +m65N)ds. . ..... 1): 
in Vektorform lautet sie 
DEREIERR EL BB . 2 3 rer 
Die äußere Ergänzungsarbeit ist 
0, =\v-öpds+P- IHM +[P. NR]. . . 2. 2 22.2. .(82), 
wobei die eckigen Klammern andeuten sollen, daß von den darinstehenden Ausdrücken die 
Randwerte an dem irgendwo abgeschnittenen Stab zu nehmen sind. Die Nebenbedingungen 


(3.1) und (3.2), geschrieben in den virtuellen Lasten, Schnittkräften und Momenten, multipliziert 
mit den zugehörigen vektoriellen Lagrangeschen Multiplikatoren 4,, Ay liefern 


22,0, ÖR HIN HL OEM HEN IM ds. . . . .. (3). 
Integriert man (5.3) partiell, so kommt 
Ö Ey G,=[Arö M]+ [Ar Rt] Hin ö Mt li AnxXt-ON+A:-öp) d =’; (5.4). 
Setzt man nun in (4.7) ein, so müssen sich die ausintegrierten Bestandteile gegen die ent- 
sprechenden Anteile in (5.2) wegheben; daraus folgt, daß A} mit dem Drehvektor d, A, mit 
dem Verschiebungsvektor dv identisch ist; die Arbeit der äußeren Kraft öp fällt dann 
ebenfalls heraus. Es bleibt also 
\EOMHEnt-NJds= \d-IM+(d -DXE)-ON]ds. 
Nach dem üblichen Schluß der Variationsrechnung (wir dürfen 9 M und N, weil die Neben- 
bedingungen (3.1) und (3.2) durch den Lagrangeschen Ansatz abgegolten sind, als vonein- 
ander unabhängig ansehen) folgt daraus 
t=$, | 
me dxXt| 
wir brauchen nun diese Gleichungen nur noch in Komponenten anzuschreiben, um die 
gesuchten Beziehungen zwischen den vier Verschiebungsgrößen u, v, w, y und den vier Ver- 
zerrungsgrößen £, #,, %,, 9 zu erhalten. Aus der ersten Gleichung ergibt sich 
=y —k,v, +k,»,, 
„=n hr, +k,x, | De re a 


(5.5): 


=», —ky+h», 
aus der zweiten 
EB. ee wi 2 


v, =—w—hv +ku, | 
; Ra ee hl AR 
„„=v—hw+k,u | 


Die Gln. (5.8) enthalten die Annahme y,=y,=0; sie hat zur Folge, daß die Winkel »,, v, (im 
Gegensatz zu y) keine selbständigen Verschiebungsgrößen mehr sind (noch weniger natür- 
lich Verzerrungsgrößen). Trotzdem ist ihre Beibehaltung als Zwischenveränderlicher (Ab- 
kürzungen für die rechten Seiten (5.8)) zweckmäßig, weil auf diese Weise die GIn. (5.6), die 
sonst eine ganz undurchsichtige Form annehmen würden, unmittelbar anschaulich deutbar sind. 

In den Anwendungen pflegt man meistens auch noch die Längskraftarbeit neben der 
Biegearbeit zu vernachlässigen. Das würde bedeuten, daß man e„=0 zu setzen hätte, und 
damit bestände nach Gl. (5.7) auch noch eine Bindung zwischen den Komponenten «, v, w 
des Vektors der Verschiebung d. D. h. ein Balken, dessen Mittellinie als undehnbar ange- 
nommen ist’), erleidet nur drei Verzerrungen (9, x#,, #,) und sein Element hat daher auch 


?, Die Gln. (5.5) lauten in diesem Fall: =, 
b+tXd=l; 
==®. 


man beachte die formale Analogie zu den Gleichgewichtsbedingungen sr KR 
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nur drei Freiheitsgrade der Bewegung: z. B. v, w, y. Entsprechend könnte man bei einem 
Stabe, bei dem das Trägheitsmoment J, sehr viel größer ist als J,. und als der Torsions- 
widerstand Jr für Näherungsüberlegungen #,=0 setzen. Bei einem solchen Stab wären dann 
auch y, ve und w aneinander gekoppelt. 

Will man umgekehrt das vollständige System der Beziehungen zwischen den sechs 
Verschiebungsgrößen dv, d und den sechs Verzerrungsgrößen f, &, y,, y, erhalten, so hat man 
auszugehen von einem entsprechend erweiterten Ausdruck für die innere Ergänzungsarbeit ö C;. 
Bezeichnen y, bzw. y, die n,- bzw. n,-Komponenten der durch die 
Querkräfte Q,, Q, bewirkten Abnahme des ursprünglich rechten Winkels 
zwischen dem körperfesten Tangentenvektor und der Schnittfläche s= const 
(s. Bild 4), so wird 


9 - HM +, 6N +60, —y,0Q,)ds, 


und an Stelle von (5.8) tritt 


„= w—hv-+ku, | 


s Er TE Te Bild 4. Schubverzerrtes 
"a v. + — hw-+ k, u | Balkenelement. 


‘ 


6. Zur Deutung der Gleichung f=2. 


Man kann sich nachträglich natürlich leicht die Bedeutung der einzelnen Terme in den 
Gl. (5.6) bis (5.8) klarmachen. Wir betrachten das Beispiel eines ursprünglich ebenen, durch 
Querlasten zu einer Raumkvrve verbogenen Ringes. Erteilen wir einem Element eine 
Verdrehung v, (um die Normale n,), dem Nachbarelement eine Verdrehung v, + dv, , so ist die 
Krümmung der Stabmittellinie um die Achse n, gegeben durch dv,/ds, ganz unabhängig davon, 
daß die Normale n, längs s nicht sich selbst parallel bleibt. Trotzdem aber ist d»v,/ds nicht die 
mechanisch interessierende Krümmungsgröße, die in das Arbeitsprodukt eingeht. Denn 
die Ausgangskrümmung k,=1/R des Ringes hat zur Folge, daß auch Drehungen y eines 
Elementes Längsdehnungen hervorrufen. Dabei wird ein Teil der Fasern gelängt, der 
andere Teil gekürzt; die diesen Verzerrungen entsprechenden Längsspannungen o, bilden 
ein Kräftepaar um die n,-Achse, das proportional ist dem Umstülpwinkel y und der Vor- 
krümmung k,°®). Als die mechanisch-interessierende Krümmungsgröße hat daher zu gelten 


nn 00, 7 0 


in Übereinstimmung mit der zweiten Gl. (5.6). 

Das Wechselspiel zwischen y und », wiederholt sich mit vertauschten Rollen, wenn wir 
nun nach dem Torsionsanteil 9 des Krümmungsvektors f fragen. Ein erster Anteil rührt her 
von der Änderung / des Umstülpwinkels, ein zweiter von der Verdrehung, die zwei Nachbar- 
querschnitte erleiden, wenn jeder um die jeweilige (der anderen nicht genau parallele Nor- 
male) n, gedreht wird. Denn der räumliche Winkel, den zwei ursprünglich parallele Gerade 
in den beiden Stabquerschnitten vom Abstand ds auf diese Weise erhalten, erscheint, ge- 
sehen in Richtung der Tangente, als ein Winkel von der Größe —v, (k,d s) (wenn man gleich 
die in ds quadratischen Glieder wegläßt), liefert also zur Verwindung ® den Beitrag — v,k,; 
insgesamt ist also zu setzen 

ee ET TE 


in Übereinstimmung mit der ersten Gl. (5.6) °). 


7. Die Elastizitätsbeziehungen für den räumlich gekrümmten Stab. 


Die Beziehungen zwischen den Schnittresultanten und den Verschiebungsgrößen lassen 
sich auf Grund des Hookeschen Gesetzes mit Hilfe der GlIn. (5.6/8) angeben. Achtet man 
nicht auf den Unterschied zwischen „Krümmung der Mittellinie“, „Mittlere Krümmung des 
Stabelementes“ usw., so lautet mit 


F=\\dydz, J, @F=\\z’dydz, J,=\\yzdydz, J a F=\\ydydz, 


s) Siehe die ausführliche Herleitung dieses Sachverhaltes bei Biezeno und Grammel: „Technische Dynamik“. 
Berlin 1939, S. 396 ft. 

9) Bei der Berechnung des Umstülpens von Kreisringen läßt Grammel (a.a.O. S. 401) den Term vr; kg (in unserer 
Schreibweise) weg, weil der Ring hauptsächlich Verdrehungen erleidet. Bei kontinuierlich verteilten Umstülpmomenten 
eutsteht daraus kein merklicher Fehler — und die Rechnung wird einfacher. Für die Formel (12) auf 
S. 44 ergibt die genaue Rechnung allerdings an Stelle von 1/7 den Faktor 7/8; d. h. im Extremfalle der Belastung 
durch 2 Einzelumstülpmomente bewirkt die Vernachlässigung des Gliedes r; ka einen Fehler von - 20%, für den Um- 
stülpwinkel. — In den Kreisringuntersuchungen der vorhergehenden %$ des Biezeno-Grammelschen Buches werden 
die vollständigen Gleichungen für die Krümmungskomponenten benutzt (s. S. 329, Formel (3)). 
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und dem (auf Grund der St. Venantschen Torsionstheorie der betr. Querschnittsform zu be- 
stimmenden) Torsionswiderstand Jr das Elastizitätsgesetz 


N=EFe„, Mr=GJr9, 
M=EJ,s,-EJ,„%, M=EJ,,—EJ,x, 


Will man die angedeutete Unterscheidung — wenigstens soweit sie den Zusammenhang 
zwischen den Längsspannungsresultanten N, M,, M, und ihren Verzerrungen betrifft — be- 
rücksichtigen, so muß man von der sorgfältigeren Definition des betr. Anteils der inneren 
Ergänzungsarbeit ausgehen: 


(7.1). 


sc |tfledodaydands. . ......2:. 02, 


wobei unter ds=ds(y,2) die Länge eines Faserelementes verstanden werden muß. 
ds stimmt wegen der Vorkrümmung des Balkens nicht genau überein mit der Länge ds 
der Schwerpunktsfaser; es gilt, wenn mit y und z die in den Richtungen n, und n, ge- 
messenen Abstände von der Schwerpunktfaser bezeichnet werden, 


dit kNlitkddse......°7:,2.. UP) 


Entsprechend hat man €, durch Integrale vom Typus | \uöodydz auszudrücken und die 
Gleichgewichtsbedingungen für ein Element dsdydz zu formulieren. Um von den so sich 
ergebenden örtlichen Verzerrungsgrößen auf die uns interessierenden Mittelwertsaussagen für 
die o„.resultanten zu kommen, muß man irgendwelche Hypothesen einführen über den 
Verlauf der Größen öo und e mit y und 2. Zwei Annahmen (von denen beim ursprünglich 
geraden Balken jede eine Folge der anderen ist) sind insbesondere sinnvoll: daß die Ver- 
schiebungen linear verlaufen mit y und z („ebenbleibende Querschnitte“), oder daß die Ver- 
zerrungen linear verlaufen mit y und 2 („geradlinige Spannungsverteilung‘“). 
Wir beginnen mit der zweiten Hypothese, indem wir ansetzen 


s-N,, öM, BL. | 
F BE (7.3). 
ee +8% —y%, 
Darin sind öN, öM,, öM, definiert durch 
öN=\\öodydz, öM,=\\öozdydz, öM,=—\\doydydz . 5° 


und zwar ist in (7.3) zur Vereinfachung der Schreibweise angenommen, daß n,, n, in die 
Hauptachsenrichtungen des Querschnitts fallen. Setzen wir (7.3) und (7.2’) in (7.2) ein, und 
führen die elementaren Integrationen aus, so ergibt sich 


döCko)=\EON+R,ÖM, +%,0M))ds. . . 2.2... 174), 
mit 
e = +uk,%, +ük,%,, 
„,—=%+k,6, BE ae (7.4'). 
,=h,+tk,e, 


Der Ausdruck (7.4) hat die gleiche Form wie der entsprechende Anteil in (5.1). Die an diese 
Gl. ein Überlegungen bleiben dieselben, nur daß an die Stelle der Größen 
Em; %,, %, Jetzt die Größen €,x,,x%, treten, die durch (7.4’) mit den Größen &,,%,,%, des 
Verzerrungsansatzes zusammenhängen. Anstelle von (5.6) a (5.7) kommt also 


»,=y, —hw,+tk,y, 
»=»,+hv, —k,y, ee Te et 
e=u —k,v+kw | 
Daraus folgen mit o=Ee für die Elastizitätsbeziehungen 
N=\\odydz=E(ledydz 
=EF,=EF(—-ük,%,—i%k,%,), 
M, -E\\ezdyde=EJ,%3,=EJ, (#, -k,&), 





M,=-—-Elleydyde=EJ,,=EJ, (#%, —k,&,). 
Wenn man die Glieder ik}, ik} gegen 1 vernachlässigt (was man bei einem „Balken“ immer 
kann), läßt sich auch schreiben 

N=EF(e—ük,%,—ük,x,), | 


zu R = u 7.6). 
M=EJ(—kd, M=RJ,@,—k,)] is: 
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\ Die Gln. (7.5) und (7.6) sind die Blastizitätsbeziehungen für den krummen Stab, wenn n,,, n, Haupt- 
e- » k . 5 ” - . - en 
achsen sind, und wenn die Längsspannungen sich linear über den Querschnitt verteilen. 
Mit k, =h=0 erhält man für den Sonderfall des ebenen, in seiner eigenen Ebene ge- 
N bogenen Ringes: 
2,=w,=v+(k,u), e=u—k,r, 
n \i 
m M=EJ@+ko+kn), N=EF@-keo)—EJI+(k,u)), 
en bei k,= const = —1/a speziell 
9), M=EJ( + vla?) 
EJ 
1B. N=EF(+vla) — — (# — ula) aa 
ls a (1.6). 
5. are KEnR EJ . ; 
oder (wegeni?Ja <1)=EF(u+v/a) — z &+ v/a?), 
>’ 
-) Die erste Gl. (7.6) ist ein aus der Theorie des ebenen Kreisringes wohlbekanntes Ergebnis, 
lie Gehen wir aus von der anderen oben erwähnten Hypothese über den Formänderungs- 
ch zustand der Querschnitte, so haben wir für den Verschiebungsvektor u=u(y,2) anzusetzen: 
für 
en u=d+dxX(yn,+tzu,). 
ch Dieser Ansatz liefert für die in (7.2) vorkommende Größe eös einen in 9,2 exakt linearen 
FR Ausdruck, den wir in der Form 
r- 
Br Di » : rien A 
schreiben können, worin & die Dehnung der Schwerpunktsfaser ist. Für die virtuellen 
Spannungen, die (wenn das Hookesche Gesetz gelten soll) nun nicht mehr exakt linear 
3). verlaufen können, müssen wir allgemeiner ansetzen 
1 % öM, öM, ,„6öL, „OL, By 
do=z(öÖN+ÖL,+öL,)+z T, FE —2 J, —yY gi Pi SE 
3°), Dabei sorgen die Zusatzterme in der ersten Klammer dafür, daß die Beziehung 
die öN=\\öodydz 
ınd 3. 
erhalten bleibt; und durch die Punkte ist angedeutet, daß im allgemeinsten Fall noch 
4 Glieder mit yz, y’z,y 2? usf. aufgenommen werden müßten, die aber auf die drei Resultanten 
4), öN,öM,,öM, keinen Einfluß haben, weil 4,2 Querschnittshauptachsen sein sollen, und 
weil Querschnittsfunktionen von höherem als dem zweiten Grade in der Biegetheorie des 
„dünnen“ Stabes nicht berücksichtigt werden. Setzen wir (7.7) in (7.2) ein, so erhalten wir 
4’). (bis auf Glieder, die höhere Querschnittsfunktionen enthalten) 
8C0;()=|(6N+%,8M, +2,08 M,)ds 
ese 
ßen Es übernehmen also jetzt die Größen &,%,,%, die Rolle der &„,%,,%, des Ansatzes (5.1), und 
des anstelle von (5.6) und (5.7) tritt daher 


=, —h „+h,y 


5) =, +hr, —k,y tan) BE NEE BT re 
> e=ü—k,v+k w 
Die Beziehungen zwischen den Schnittkräften und den Verschiebungen erhält man 
daraus, indem man (7.7) und das Hookesche Gesetz o=Ee in die Definitionsgleichungen von 
N,M,,M, einführt. Aus (7.7) folgt 
x 2 GE + u a f \ er 
e=(E +2%, -YyR) 57 Ü-kyü+k =(£+2%, —y%,)(1l+k,)(1-—k, 2); 
daher ist 
ner N E\\edydz=EF( —Ük,%, —ük,%,) | 
M,=El\ezdyda=EJ (#&, —k, 6) (7.9). 
7.6). M, E\ \eydydz=EJ,(%, — k, €) | 


Yııaa 
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Die Gleichungspaare (75), (7.6) und (7.8), (7.9) sind bis auf die Bezeichnung der 
Zwischenveränderlichen &,x identisch. Es ergibt sich also, daß die Elastizitätsbeziehungen 
des krummen Stabes (abgesehen von Termen, die in (ik) von höherem Grade sind) genau 
dieselben sind, ob man nun von der Annahme linearer Spannungs- oder linearer 
Verschiebungs- Verteilung ausgeht. 

Von der Einschränkung, daß n,,n, Hauptsachsen sein sollten, kann man sich natürlich 
auch befreien. Die Überlegungen bleiben genau die gleichen, nur muß man die d o-Ansätze 
(7.3) bzw. (7.7’) so erweitern, daß auch bei J ,=#0die Definitionsgleichungen (7.3’) bestehen bleiben. 
Die Ergebnisse bezüglich der Krümmungsgrößen ändern sich nicht, weil das Skalarprodukt 
öM,x,+öM,x, eine Drehungsinvariante ist; d. h. an die Stelle der zweiten und dritten 
Gl. (7.6) [oder (7.9)] treten die Beziehungen 


M=EJ, EJ,5, M=EJI,S,—-EIsä, 


2, Dich 
bs R 14.9). 
N, ®, en, — kt 


mit x x 


1 1 


Und bei e, tritt ein zahlenmäßig belangloser Term von der Forni 


(J,./F) (k, 2, + k, %;) 
zusätzlich auf. 

Die ganze in diesem Abschnitt gegebene Verfeinerung der Formeln (7.1) ist vom 
praktischen Standpunkte aus unwichtig. Die Unterscheidung zwischen einer Krümmungs- 
größe #%, der das Moment proportional ist und einem mittleren Kontingenzwinkel pro Längen- 
einheit x oder %, an dem das Moment Arbeit leistet, fällt in sich zusammen, wenn man die 
technisch fast immer zulässige Annahme macht, daß die Schwerpunktsfaser ungedehnt bleibt. 
Und auch bei 2=#0 ist die Korrektur vom Typus k & unwesentlich: denn der zweite Term 
im Ausdruck für die Biegedehnung 


seti2)=(iR) (ih) 


ist beim dünnen Stab (bei dem ix von mindestens der gleichen Größenordnung ist wie ik) 
gegen den ersten klein wie e gegen 1, und beim diekeren Balken (der z. B. eine Stützlinien- 
last trägt, und daher wenig verbogen wird) sind zwar beide Anteile von derselben Größen- 
ordnung, dafür aber beide so klein, daß eine Berechnung der Momente sich erübrigt. Ent- 
sprechendes gilt von den beiden Anteilen, aus denen die Längskraft N sich zusammensetzt. 
Außerdem kranken alle Versuche, den dickeren Balken durch eine Näherungstheorie unmittelbar 
(d. h. ohne Zurückgreifen auf die dreidimensionale Theorie) erfassen zu wollen, daran, daß 
über den Verlauf von e=e(y,z) eine willkürliche Annahme gemacht werden muß, und daß; 
die ganze Verfeinerung, wie wir sie hier gegeben haben, dem Schubanteil der Schnittspan- 
nungen überhaupt nicht zugute kommt. Man müßte dann auch dort auf eine Torsionstheorie 
des gebogenen Stabes zurückgreifen — ein Aufwand, der sich für die Zwecke der An- 
wendungen nur selten lohnen würde. 






















8. Anhang. Eine andere Herleitung der Gleichung td. 
Mit Benutzung des Arbeitsbegriffes war im Abschnitt 5 gezeigt worden, daß die Kom- 






ponenten des Vektors d gerade die mechanisch interessierenden (d. h. den Momenten proportional 





zu setzenden) „Krümmungs“größen sind. Wir wollen anhangsweise für die Beziehung f=d 
noch eine andere Herleitung geben, die zwar etwas Vektor-Formalismus erfordert, dafür aber 
den Vorteil hat, daß bei ihr an die Stelle der im Begriff „Formänderungsarbeit“ etwas ver- 
kappt steckenden mechanischen Überlegung eine unmittelbar anschauliche Aussage über 
das Elastizitätsgesetz tritt. 

Nach (2.4”) ist der Krümmungsvektor einer Kurve mit dem begleitenden Dreibein t,n,,n, 
derjenige Vektor, dessen äußeres Produkt mit t,n,,n, die Ableitungen t, it, it, liefert: 













Ba EN, . nee te A 





Multipliziert man diese Gl. äußerlich mit t,n,,n, und addiert, so erhält man für o eine 
Darstellung in expliziter Form: 






EIKE KU HF, KK EX) +, XlOXn)+n,X(oXn,) 
=0 -t(t-0) +0. min, -o)+n -n,(n,-0)=20; 





l . ö 
v=-Zltxt+rn, x, +m,Xn,) 
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Die eine oder die andere Form der Krümmungsdefinition kann man benutzen, um den 


Krümmungsvektor ö des verformten Systems zu bestimmen. Bezeichnet t,ft,, ft, das Achsen- 
kreuz, in das t,n,,n, infolge der Verformung übergeht, so gilt 


tet Hd, en, HH Xn, nm, HdXN.. . 222222088), 
und, da d nach Voraussetzung eine kleine Größe sein soll, umgekehrt 
tet IX, nen -dXn, neu IX... 2.2.02. 88. 


Wir bilden nun formal die Größe t, (f,,n, folgen daraus ohne Rechnung) mit Hilfe von (8.3) 
und (8.3”). Es ist 


t-t+dXt 4 DXL=0XEH+HdDXE Hd (0X). 


Ersetzen wir rechts t durch t (unter Streichung der in d quadratischen Terme), so kommt 


t=0XÜ-DXDHDKIHDIKINE) 
— (+) XL 0oXPXIY)+dX (0X?) 
AETRKARE ir ra ee 
t läßt sich ausklammern; die Größe 
D=0+dXo+d 


ist also schon die gesuchte „neue Krümmung“, und wir haben das Ergebnis nur noch mecha- 
nisch zu deuten. Zu diesem Zweck machen wir ein Gedankenexperiment. Biegen wir den 
krummen Stab durch geeignete Momente elastisch gerade, so sind diese Momente propor- 
tional den Komponenten der Krümmungsänderung, d.h. den Komponenten von od, bezogen 
auf die körperfesten Achsen t,n,,n,. Biegen wir den Stab ausgehend von der gestreckten 
in die durch t, it,, ft, gekennzeichnete Lage, so müssen wir Momente aufbringen, die wieder 
proportional sind den Komponenten der Krümmung, bezogen auf die neuen körperfesten Achsen 
t,it,,ft,. Die Differenz beider Momenten-Tripel ist das wirkliche elastische Moment, die Differenz 
der körperfesten Krümmungskomponenten also die mechanisch interessierende Krümmungs- 
größe. Um diese Differenz bilden zu können, müssen wir daher die t,ft,,it,- Komponenten 
des Vektors 0+d%X0-+d bestimmen. Bezüglich d ist keine Rechnung nötig, denn von dort 
her kommen beim Übergang von t,ı,,ı, auf t,fi,,f, nur in d quadratische Terme. Die 
Komponenten von 0+dX.o bezüglich t... aber sind gerade wieder die alten Komponenten 
h, k,,k, von o bezüglich t.... Denn wenn man die 3 Gln. (8.3) der Reihe nach mit h,k,,k, 
multipliziert und addiert, so erhält man 


BEE EER EN E TFT  R 
Die drei Komponenten der Krümmungsänderung f sind daher 
h+d-.Yd—h=d-t, (k,+d-n)—k,—=d-n,. (k,+d-n,)— k,=d-n, 
und f=#9t+x,n,+x%,11, ist also die totale Ableitung des Drehvektors d: 
=». 


Ähnlich wie o nach Gl. (8.2) kann man übrigens auch d aus drei unmittelbar „anschau- 
lichen“ Winkelvektoren aufbauen. Nach (8.3) ist 


Dt=t--t=dXt, Dn,=dxXn, Dn,=dxXu, 


und daraus folgt durch äußerliche Multiplikation mit t,n,,n, und Addition 


D- (tx Din XDn,+n,XDn) . ar [ae ee A 
284 
15* 
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Die Spannungserhöhung durch Kreisbogenkerben *). 


Von E. Weinel in Jena. 


Der ebene Spannungszustand in einem von zwei Kreisbogen begrenzten 
Bereich wird mit Hilfe von Bipolarkoordinaten durch Fourierintegrale 
dargestellt. Als Anwendung wird die Spannungserhöhung in Kreisbogen- 
kerben und -Löchern durchgerechnet. 


Der Spannungszustand in einer Halbebene mit Kreisbogenkerbe wurde kürzlich von 
C. Weber!) mit Hilfe eines Abbildungsverfahrens untersucht. Unbeschadet der Tragweite 
des von C. Weber entwickelten Verfahrens fortschreitender Näherung kann dieses Problem 
auch sehr einfach elastisch-genau gerechnet werden, wenn man die dort vorkommenden Ko- 
ordinaten 4, d, die man unschwer als bipolare Koordinaten erkennt, als angepaßtes Bezugs- 
system verwendet. Darüber hinaus lassen sich auf diesem Wege alle Spannungsaufgaben, 
für Bereiche, die von zwei Kreisbogen begrenzt werden, unmittelbar behandeln. Diese Rechı- 
nungen schließen sich Schritt für Schritt an die früheren Anwendungen bipolarer Koordinaten?) 
an und sind dadurch weitgehend vorgezeichnet. Der Aufsatz von C. Weber gibt mir Anlaß, 
diese Ergebnisse, auf die bereits a. a.O. in der Schlußbemerkung hingewiesen wurde, nun- 
mehr mitzuteilen. 


1. Der Spannungszustand in einem von zwei Kreisbogen begrenzten Bereich. 
Durch die komplexe Funktion 


ar 
2 +iy=alotg er 


Sinı 
2=ud aa —, 
GojA — cosg 

sin c 

ya ’ 


Coj A — cosg 
wird ein von zwei Kreisbogen begrenzter Bereich der x, y-Ebene auf einen Parallelstreifen 
der 4, g-Ebene abgebildet. 


Für die Abmessungen des Bereichs (Bild 1) gelten 
die geometrischen Beziehungen: 


r,la=1/sing,; tj/a=ctg (g,[2); “7 
t,ja= h 



























r,ja=1/sing,; = ctg (p,/2). ug 
Mit 9,=0 ergibt sich der Sonderfall der Halbebene IN | 
mit Kreisbogenkerbe. 4 ie | 
An Stelle der Airyschen Spannungsfunktion .. Pi 7 \ fi 

F(i,g) wird entsprechend dem Vorgang a.a.0. eine 1 5 
Hilfsfunktion f (A, @) gebildet, die mit der Funktion F* Ä \ F 
bei ©. Weber wesensgleich ist. Die Beziehung zur Er a o 
Airyschen Funktion wird durch 11. 

2 flag) Sn | & 

ru0=e Coj A — cosgq 
vermittelt, und die Spannungskomponenten berechnen 
sich aus: Bild ı. 








*; Der Verf., der zur Zeit im Felde steht, hat keine Korrektur der Arbeit gelesen. Zu dem ihın nach Eingang 
des Manuskriptes übersandten Aufsatz von Isibasi: Stresses in a Semi-infinite Plate with a Cireular Note under Uniform 
Tension, Memoirs of the Faculty of Engineering, Kyushu Imperial University, Fukuoka, Japan IX, 2 (1940), S. 131 
bis 143, der das gleiche Problem in ähnlicher Weise behandelt, schreibt Herr Weinel: 

„Die Arbeit von Isibasi deckt sieh sachlich vollständig mit der von mir unter Ziff. 2 angegebenen Lösung 
des von Herrn Weber behandelten Problemes. Isibasi nimmt indessen die-übliehe Aufspaltung der Spannungsfunktion 
vor, auf deren Entbehrlichkeit ich hingewiesen habe; I. gelangt infolgedessen zu einer ziemlich unübersichtlichen 
Formel, die natürlich bei entsprechender Umformung in die von mir gewonnene einfachere Form übergeht.“ Willers. 
















1,0. Weber: 





Halbebene mit Kreisbogenkerbe. Z. angew. Math. Mech. Bd. 20 (1940), S. 262 bis 270. 











2, E. Weinel: 
S. 276 bis 287. 


Über einige Randwertaufgaben der Elastizitätstheorie. Z. angew. Math. Mech. Bd. 17 (1937), 
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a a EEE 
o,= (Cji—cosy) u Sini,, sing 5 Eofi fii,g), 
EN: a BE | SO of 7 
o,= (KojA— cosg) ApE Sini,, sing + coso fl, p), 
re o”f 
„=—iR OS z a 
7, (Eoj A — cosg Eyry 


Während nun bei der Behandlung anderer Aufgaben a.a. 0. die Funktion f(4, 9) als Fourier- 
Reihe auftrat, erscheint sie im vorliegenden Fall als das entsprechende Fourier-Integral ?) 
Es 
fii.p) = \ cosi t[ Alt)cosy . Cojy un Bihsing ° Sing t+-((t)cosp: Zing t+ Dit)sin q GCoiy td t 
" 
E. 


+ \sinät[A(hcosg Cojgt+Bitsing Sing t+Ü(deosp Sinyt+ Didsing Kojptlat. 
0 
Wegen der Wesensverwandtheit von Fourier-Reihe und -Integral ist die weitere Behandlung 
durch die entsprechenden Rechnungen a. a. O. vollständig vorgezeichnet. Das gilt insbesondere 
auch für die zu beachtenden Eindeutigkeitsbedingungen für den Verschiebungszustand. 
Bei der Anpassung an die Randbedingungen sind dabei an Stelle der a. a. 0., S. 281, 
benutzten Reihen die folgenden Integraldarstellungen zu verwenden: 


sing _ „Sina —pit , 
Cof A — cosg -2| Sinat oositdi, 
0 
' GCojA+ cosg (a/2—g)tcosit w 
NSi-csp ” Cojr/2t re 


0 

Die Einzelheiten der Rechnung sind durch die entsprechenden Ausführungen a.a.O. fest- 
gelegt, so daß sich ein weiteres Eingehen darauf erübrigt. Es sollen daher sogleich spezielle 
Ergebnisse soleher Rechnungen mitgeteilt werden. 

2. Die Halbebene mit Kreisbogenkerbe (Bild 2). 
Der Ausdruck für die Spannungsfunktion f(4,y) lautet: 
Es 

3.1 „[Sing,tsing Sina, —p)t—t-sing,sin(p, —p)Sinpt _, 
fi, g)- P\ Sin’g,i_Faing, cosAitdt. 
0 


Die Spannung o am Kerbgrund berechnet sich aus 
E- 
“ tsing, Cojg,t — cos p, Sin gp,t 


Sin’go,t —-#-sin’g, a, 


o—4p- sin?’ (p,/2)\t 
N 
und die numerische Auswertung liefert: 
o=10° 15° WW 5 0 
o/p=1.W 2,088 3,063 3,747 3,990 (genau), 
1.00 — 3,07 — 400 nach €. Weber. 


Diese Zahlenwerte bestätigen die Brauchbarkeit der Weberschen Methode. 


R R 
,  Pertrhenten HALRLLANEEN 


* ws 


y,"135° 9%? 9:45? %o=0° 


KRZERE EEET ET ZZ ea EEE EIER EI EEE EIER EEETET 
KzZ R R 























E 


Bild 2. Halbebene mit Kreisbogenkerbe. Bild 3. Gelochte Ebene. 


3) Die Aufspaltung der Funktion /f in zwei Anteile, wie sie a.a.O. vorgenommen wurde, hat sich als unnötig 
erwiesen, da auch das Verhalten im Unendlichen durch den Fourieransatz miterfaßt werden kann. Die vorkommenden 
Reihen und Integrale sind dann zwar nur bedingt konvergent, ohne daß sich daraus Schwierigkeiten für die Reehnung 


ergeben. 
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3. Die Ebene mit symmetrischem, von Kreisbogen gebildetem Loch (Bild 3). 
Die Spannungsfunktion f(4,g) hat gültig für 2 >0 - die Gestalt: 


” 


Coj y,tsin p Sin(p, — p)t tsing, Cojytsiniy, 9) 


füp=p \ Sinp,tCojg,t+tsing,cosgp, cositdi 
0 
i Sin g,tcos gq Sing eo pt—E sin 9, Coj q t sin (p,—g) cosit dt 
+a\ Sing,t&ofp,t+tsing,cosp, tP+1) 
u 
E Sing,tsingpk&oj(p, - p)t— sing, Sinptcos(p, - y) cosit di 
ip Sing,t&ojp,t+tsing,cosp, +1) 


0 
Der Beiwert q ist aus der Eindeutigkeitsbedingung für den Verschiebungszustand zu be- 
stimmen; diese Vorschrift lautet: 
= “ . 
t sin? Yo ER \ Sin’p,t -Psin’o, _ dt 
p le Sinp,t&ofg,t - sing, c0sp, .. Sing,t&ojg,t+tsing,ecosp, fR+1) 
0 


Für die aa am Grund der Kerbe folgt: 


„einst: sing, —k&ojyyt: c0S ps 


! 
Sin y,tCofp,t+tsinp,cosgy, . 


.: 4p - sin? (Pol al 
0 
+4 q- sin? (9,2) \ 


sin @,° Sinp,t e. di 
Sing,tCofp,t+tsinp,cosp, 





Für 9, =®%®° (Kreisloch) ergibt sich daraus das bekannte Ergebnis o/p =3. 

Im Grenzfall 9,=0 (Bild 3) erhält man o/p = 3,84. 

Die Spannungserhöhung in der gelochten Ebene ist also nur unwesentlich kleiner als 
bei der gekerbten Halbebene. 261 


Halbebene mit Kreisbogenkerbe. 
Von C, Weber in Dresden. 


Die von mir in dem gleich überschriebenen Aufsatz des vorigen Bandes 
dieser Zeitschrift angewandte Methode wird als Singularitäten-Methode 
charakterisiert, im Gegensatz zur Methode der Fourier-Entwicklung. 
Die Berechnung der Spannungserhöhung wird für verschiedene Kerbtiefen 
mit größerer Genauigkeit durchgeführt und die Ergebnisse zusammengestellt. 


In meinem Aufsatze in Heft 5 dieser Zeitschrift, Bd. 20 (1940), S. 262 bis 270, zeigte 
ich einen gut konvergierenden Weg zur Berechnung der Spannungen der Halbebene mit 
Kreisbogenkerben, insbesondere der Spannungserhöhung im Kerbgrunde. Da die dort an: 
gegebenen Zahlenwerte ohne Rechenmaschine gefunden worden sind, will ich ergänzend hier 
die Berechnung mit größerer Genauigkeit und auch für weitere Kerbtiefen angeben. Die 
untersuchten Fälle sind in Bild 1 dargestellt. 


pe 171937 
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Vorweg möchte ich noch einiges zu der von mir verwandten 
| Methode sagen. Hat man eine belastete Halbebene nach Bild 2a 
ne a oh oder b oder einen belasteten Streifen nach Bild 2e oder d zu unter- 
Er suchen, so stehen hierfür zwei Methoden zur Verfügung. Da für 
N cos-förmige Belastungen der Halbebene und des Streifens die 
b Lösungen bekannt sind, kann man die Lösungen bei periodischer 
| Belastung (Bild 2a und ce) durch Fourier-Reihen, bei nicht- 
y y periodischer Belastung (Bild 2b und d) durch Fourier-| ntegrale 
— N darstellen. 

! f Für die Halbebene wird man statt dessen geschlossene Lösungen 
bevorzugen, hierzu benutzt man die bekannten Potentialfunktionen, 
die in den Angriffspunkten der Kräfte Singularitäten aufweisen. 
Ich will diese Methode als Singularitäten-Methode bezeichnen, 
die vorherige als Methode der Fourier-Entwieklung. Auch 
a bied. für den Streifen läßt sich die Singularitäten-Methode anwenden, wie 
ich in meinem Aufsatze „Streifen mit Einzellast“') zeigte. Hierbei 
treten außer den Singularitäten in den Angriffspunkten am Streifenrande noch weitere Singu- 
laritäten im Gebiete außerhalb des Streifens auf. Von diesen Singularitäten wird man bei 
numerischen Berechnungen nur die nächstliegenden berücksichtigen und den nachbleibenden 
Spannungszustand entweder nur näherungsweise bestimmen oder, falls dieses nicht genügt, 

durch Fourier-Entwieklungen, die jetzt wesentlich besser konvergieren, berechnen. 

Nach der Singularitäten-Methode hat Herr H. Swoboda, Dresden, die kreisrunde 
Platte mit Einzellast und mit auf einer kleinen Kreisfläche verteilten Last untersucht). 

Wird bei der Ebene mit zwei kreisrunden Löchern oder mit einem Loche, das durch 
zwei Kreisbögen begrenzt ist (hierzu gehört die Halbebene mit Kreisbogenkerbe), das Gebiet 
in bekannter Weise konform auf einen Streifen abgebildet, so erhält man einen Streifen mit 
periodischen Randbedingungen bei zwei Löchern und einen Streifen mit nicht- 
periodischen Randbedingungen bei sich berührenden oder sich schneidenden Kreisen. 
Für sin- oder cos-förmige Randbedingungen lassen sich in beiden Fällen die Lösungen an- 
geben, so daß sich die Methode der Fourierentwicklung ohne weiteres formal an- 
wenden läßt. In meiner Arbeit „Spannungserhöhung in Blechen mit mehreren kreisrunden 
Löchern“) habe ich die Scheibe mit zwei kreisrunden Löchern behandelt, wobei das Gebiet 
auf ein Ringgebiet abgebildet wurde, das durch konzentrische Kreise begrenzt ist. Bildet 
man den Kreisring auf den Streifen ab, so gibt die von mir angegebene Lösung eine Fourier- 
Reihe. Für sich berührende Löcher geht die Fourier-Reihe in ein Fourier-Integral über. 
Schon für Kreise, die dicht nebeneinander liegen, konvergierte die Fourier-Reihe schlecht. 
Auch konnte ich keine Fehlerabschätzung finden. Noch schwieriger wurde die Berechnung 
für das Fourier-Integral, das man für sich berührende Kreise erhält, so daß ich das nur for- 
male Ergebnis nicht veröffentlicht habe. Auch bei den Fourier-Integralen, die sich bei sich 
schneidenden Kreisen ergeben, wird die rechnerische Arbeit der Auswertung und Febhler- 
abschätzung (falls eine solche überhaupt möglich ist), m. E. entweder sehr umfangreich oder 
die Genauigkeit und Sicherheit lassen viel zu wünschen übrig. In welchen Fällen die 
Singularitäten-Methode und in welchen die Methode der Fourier-Entwicklung vorzuziehen ist, 
kann nur ein Vergleich der rechnerischen Arbeit zeigen. 

Ich habe darum mit voller Absicht für diese Scheibenprobleme die Methode der 
Fourier-Entwicklung nicht angewandt, sondern die Singularitäten-Methode entwickelt. Als 
Beispiel für die Darstellung nahm ich die Halbebene mit Kreisbogenkerbe. Es stellte sich 
heraus, daß schon wenige Glieder genügen, um brauchbare Ergebnisse zu erzielen, wobei die 
einzelnen Näherungen sich dem Endwerte von unten nähern. Verbessert man die Lösungen, 
indem man den Einfluß der Restbelastung abschätzt, so erhält man eine zweite Reihe von 
Näherungslösungen, die sich dem Endwerte von oben nähern, eine Reihe, die außerdem 
schnell konvergiert. Aus dem Verhalten beider Reihen läßt sich auch der Fehler abschätzen. 

Für die Berechnung der Spannung im Kerbgrunde ohne und mit Verbesserung wird 
für alle Kerbtiefen ein und dieselbe Reihe von Zahlen Cu, Cu Cam - - „ benötigt. 
die wie folgt definiert sind: 




















ne 
ara de ä ? 


Bild 





l a 7\" Zy. nn —1)n\? ? 
Con 9.1]. ll z) ) er (' > tes), ’ 
(& n)!o%, - } Ih =U 
1 2" +? 
( I; gen+e 
UNE (2n+1)!0%,? 
1) Z. angew. Math. Mech. Bd. 16 (1936), S. 372 bis 375. 


2) Z. angew. Math. Mech. Bd. ?0 (1940), S. 336 bis 350. 
3) Z. angew. Matlı. Mech. Bd. 12 (1922), S. 267 bis 273. 
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Zur Bestimmung der Ableitungen verwendet man am zweckmäßigsten die Potenzreihen- 
entwicklungen von tg®, und cetg d,. 
Diese Zahlen werden: 


w=l, en» = — 0,666667, Cu — 0,532599, cu = — 0,456035, 
eu = 0,405058, ce = — 0,368027, cm —0,339572, cn = — 0,316826, 
| 0.29813, Co=— 0,28235, Cu» = 0,2689, c — 705571. 
Das Verhältnis zweier aufeinanderfolgender Zahlen nähert sich allmählich dem Werte — 1. 
Im Punkte M’, Bild 2 des Hauptaufsatzes’), erhält man für die Spannungen folgende 
Näherungen ohne Verbesserung: 


(11 


> ( > 
Ka Sarg 
TEN E 


2 sin a- 


2/4 za b 
Int 3 a 3 egay Pop“ 


[#7] =(0,% lo - 
(5 (3) + @®\5 


Als Verbesserung der Näherung 9ant+n muß man die Spannung Ao.n+», die durch 
die Kraft T,,„.+,, hervorgerufen wird, hinzufügen: 


. 7 T IT 
snaz az (coSsay 
- _ —_ 


er, aT 
A 0an4 u Tont+y 


z\° ER 
ay)—-sin’ay 


Hierin ist: 


zı\an+ı 


r 7 9 P 
sın a 3 zsın a 3 
1 Re Con+n»- 


7 
eNnNntD 


a? ra 


\ 


Die Spannungen mit Verbesserung nähern sich rasch dem Endwerte oyr. Aus oyr er- 
hält man oy im Kerbgrunde mit o&,=2 und hieraus die Spannungserhöhung 


o „an 
a —=20yp sin? 4 


Für flache Kerben genügen schon die Näherungen o,, und o,, ohne Verbesserung, um die 
Spannungserhöhung bis zur vierten Dezimalstelle zu bestimmen. Für «=1, der Halbkreis- 
rille, sind die Spannungserhöhungen für jede Naherungslösung ohne und mit in Tafel 1 an- 
gegeben. 

Tafel 2 gibt die Spannungserhöhungen für verschiedene Kerbtiefen nach Bild 1. Die 
eingeklammerten Dezimalen sind unsicher°). 


0x 


Tafel 1. Spannungserhöhung im Kerbgrunde für die Halbkreisrille für verschiedene Näherungen. 





3 R 7 9 1 Geschätzter 


Index der Näherung | 1 ) Endwert 





ohne Verbesserung . . . | 2,54648 | 2,91290 | 3,01458 | 3,04747 , 3,05885 | 3,06300 \ 


ev : BE et i£ 3,065(4) 
mit Verbesserung . . . . | 3,12494 | 3,07330 | 3,06703  3,06577  3,06546 | 3,06544 ] 


Tafel 2. Spannungserhöhung im Kerbgrunde. 


714 3/: 1 sa | 12 1/4 0 








0 0,07612 | 0,29289 | 0,61731 1 1,38268 | 1,70711 | 1,92388 2 
Spannungs- 1,000 | 1,5% 2,070 2,616 3,065 3,461 3.75) | 8,987) | 3,99(8) 


erhöhung 
264 
4) Z. angew. Math. Mech. Bd. 20 (1940), S. 263. 
5) Druckfehlerberichtigung: Hauptaufsatz, diese Zeitschrift Bd. 20 (1940), S. 269, Zeile 15 von oben: 
a 3 82 F* 1 
te te 
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Über besondere Seilkurven. 


Ein Beitrag zur graphischen Analysis. 





Von Karl Federhofer in Graz. 






Graphische Ermittlung der Gleichgewichtsform eines undehnbaren dünnen 
Fadens bei alleiniger Wirkung von Zentralkräften nach dem Polarsub- 
normalen- und Krümmungskreisverfahren. Alle in die Klasse der Sinus- 
spiralen gehörenden Kurven haben die Eigenschaft, Seilkurven darzustellen 
für Zentralkräfte, die von der Entfernung des Kraftzentrums vom Faden- 
element in beliebiger Potenz abhängen. 

Ergänzung der Newendorffschen Methode der zeichnerischen Lösung 
von gewöhnlichen Differentialgleichungen beliebiger Ordnung in Polar- 
koordinaten durch eine einfache Konstruktion der Krümmungsmittelpunkte 
der Integralkurve und ihrer abgeleiteten Kurven. 

Graphische Integration einer beliebig gegebenen Funktion in Polar- 
koordinaten. 















1. Einleitung. 

Die Gleichgewichtskurven eines vollkommen biegsamen, undehnbaren dünnen Fadens, 
der nur der Wirkung von Zentralkräften unterworfen ist, sind ebene Kurven, deren Differential- 
gleichung für den besonderen Fall, daß die Zentralkräfte bloße Funktionen der Länge r des 
vom Kraftzentrum OÖ aus gezogenen Polstrahles sind, die Form hat!) 







dr 
dq = ö . . . . . . r . . . (1). 


rVr? al frn)dr +” —1 






Bezeichnet P=kf(r) die auf die Längeneinheit des Fadens 
bezogene Zentralkraft (Bild 1), S die Fadenspannung an 
der Stelle (r,9), © den Winkel zwischen Falırstrahl und 
Kurventangente, so gelten die Gleichgewichtsbedingungen 


A EEE - - \ 
(3), 













Srsin 9=( 






aus der ersten folgt 







Buh\fmärt+e,. - .-. ... BO 
rY A2e27] 
und aus der zweiten wegen sin 9 = - =; ld. 
g vr+r® Bild 1 
b R 
k\findr+0,= 51 en 2 Ve En: © 





letztere liefert mit den Abkürzungen 


k C, j 2 
— a, U ae Fe 







die obige Gl. (1). 

Gilt für die Zentralkraft das Potenzgesetz P—=kr", wo n eine beliebige positive oder 
negative Zahl bedeutet, so ist Gl. (1) unter der Annalıme #=0 durch elementare Funktionen 
integrabel mit Ausnahme des Falles n= —1; im Sonderfalle » =— 2 (Newtonsches Kraft- 
gesetz) ist eine exakte Integration auch ohne die einschränkende Annahme 5A=0 möglich’), 
ebenso in den Fällen » —=-+1*) und »=0**), die mit Hilfe elliptischer Funktionen zu inte- 
grieren sind. 

Hingegen ist einegeschlossene Lösung der für den Fall » = — 1 gültigen Differentialgleichung 


is Sr 
cc yr’(a Inr +” —1 











dy 





nicht möglich. 





1) P. Appell: Trait6 de me6canique rationelle, Bd. I, Paris 1909, 3. Aufl, S. 200. — R. Marcolongo: Theo- 
retische Mechanik, Bd. 1, S. 300; Aufgabe 4, Leipzig u. Berlin 1911. 

2) Vgl. Aufgabe 22, S. 234 bei P. Appell (vgl. Fußnote 1). — H. Egger: Z. angew. Math. Mech. Bd. 19 (1939), 
S. 319 bis 320. 

*) R. Marcolongo: Rendiconti dell’Accad. d. science Fis. e mat. Neapel (?) Bd. 6, 1892, S. 71. 
*), P. Appel: vgl. Fußnote 1), S. 200. 
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Die entsprechenden Seilkurven lassen sich aber aus den in 2 entwickelten Grund- 
gleichungen in einfacher Weise mit Hilfe bekannter zeichnerischer Verfahren konstruieren, 
wobei jede einschränkende Annahme für » oder 5 entfällt. In 3 ist für die Fällen=+1 
eine Schar solcher Seilkurven dargestellt: dort wird auch bewiesen, daß allen in die Klasse 
der Sinusspiralen gehörenden Kurven die mechanische Eigenschaft zukommt, Seilkurven 
darzustellen für Zentralkräfte, die nur von der Länge des Radiusvektors r in beliebiger 
Potenz abhängen. Der Fall eines beliebigen Zentralkraftgesetzes P=kf(r,y) ist in 4 kurz 
behandelt. Schließlich wird in 5 eine Ergänzung zur Neuendorffschen Methode der 
zeichnerischen Lösung von gewöhnlichen Differentialgleichungen beliebiger Ordnung in 
Polarkoordinaten mitgeteilt; hierbei wird die für jeden Integrationsschritt notwendige Aus- 
rechnung des Krümimungshalbmessers durch eine einfache lineare Konstruktion ersetzt, wo- 
durch das ursprünglich halb zeicherische, halb rechnerische Verfahren in ein bequem zu 
handhabendes rein graphisches Integrationsverfahren übergeführt wird. Die dabei verwendete 
Konstruktion des Krümmungsmittelpunktes einer ebenen Kurve ermöglicht, wie in 6 gezeigt 
wird, auch die einfache graphische Integration einer beliebig gegebenen Funktion. 


2. Grundgleichungen. 


Ist M, der Fußpunkt der vom Pole O zur Seillinie gezogenen Normalen und bezeichnen r,. 


S, die zu M, gehörigen Werte von r und S, so gilt wegen 9, = = gemäß Gl. (3) und (4) 
Se IE —E BERzn EN Te = ET a Er ER 
nn aan ET u EEE © 


'=rog 


Es empfiehlt sich für die weitere Untersuchung die Einführung der Dimensionslosen & und y(&) 
entsprechend der Festsetzung 


PEN. (9) 
a ag 
N, er dt 2. SDR 


die zu untersuchenden Seilkurven sind dann charakterisiert durch einen wesentlichen Para- 
meter 
kn _kr 





Hiermit geht Gl. (5) bei Beachtung von (8) über in 
/E LE ? 
Vers =1+y\giod: (12) 
:Z1 
woraus folgt a 8 
dz ’ | / E f 
?==+:V @li+y \giäadel —1 (A). 
a EA 








Zur Zeichnung der dieser Differentialgleichung entsprechenden Integralkurve bietet sich un- 
mittelbar das Subnormalenverfahren von R. Neuendorff?) an, welches auf der Ausnutzung 


\  dE : ; . 
der Beziehung z =Polarsubnormale der durch die Polarkoordinaten & und 9 festgelegten 
/ 


Seilkurve beruht. 

Vom Punkte M,(&,=1, &=0) ausgehend, trägt man (Bild 2) auf der Tangente t, des 
Punktes M, die kurze Strecke M,M, auf, liest die Polarkoordinate &, ab, berechnet den 
Wert & aus Gl. (A) und trägt dieses Maß als Polarsubnormale vom 
Pole O aus auf der Senkrechten zu OM, auf, so daß OM/ = &; nun ver- 
bindet man M; mit M,, errichtet eine kurze Strecke M,M, senkrecht 
auf M; M, und wiederholt die beschriebene Konstruktion. Wenn hierbei 
der Endpunkt M; der Polarsubnormalen unzugänglich wird, so konstruiert 





man die Kurventangente M, unmittelbar aus 80 =, wo © den Win- 
= 
kel zwischen Radiusvektor £& und der Kurventanrgente bedeutet. 
Bei starker Krümmung der Seilkurve versagt dieses Verfahren, man 
ersetzt es dann durch das Krümmungskreisverfahren von Lord Kelvin, 
das bei Verwendung von Polarkoordinaten besonders einfach wird. 





3») R. Neuenudorff: Z. angew. Math. Mech., Bd. 2 (1922), 8. 135. 
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Aus Gl. (12) folgt durch Ableitung nach 4 
et a 7 0a j Di” 2 SEE |: 
Ist £ der Krümmungshalbmesser der Seilkurve an der Stelle (r,y) und o jener der auf r, 
reduzierten (£, 9) Kurve, so gilt 


Eh I FE 2, EN 


somit wegen (12) und (13) 


0 — tr Iatde) RE a A alu 5: 
ygiE) i 

Man ersetzt nun die Seilkurve durch kleine Kreisbogenstücke, die ohne Knick aneinander- 
schließen und deren Krümmungshalbmesser aus (B) berechnet werden. Für die Seilspannung 5 
liefert Gl. (4) nach kleiner Umformung 


Ss=S,1+y\giäds) EV Zee ee ER 
| | 


3. Die Seilkurven für das Kraftgesetz y(&) = :”. 





Für die Polarsubnormale ergibt sich nach Gl. (A), wenn der Fall u = I ausgenommen 
wird, 
de y |? 
ei - ML... ann. 
dg | n+ 1 | 
für den Krümmungshalbmesser nach (B) 
gi-n u 12 
0—= — 14 ER FE Et Fr 


und für die Seilspannung nach (Ü) 
N S, e. 2 { 
n+1 
Ein besonders bemerkenswerter Sonderfall ergibt sich, wenn der Parameter » der Seilkurve 
und der Exponent n des Kraftgesetzes der Beziehung 


er 


(17). 


ea a a a a 
eehorchen; dann folgt aus (15): 

des 2 Eon 1 

dg , ai 


oder 
"rcoss(n+2)gp=1. 
Dies ist aber die Polargleichung einer Sinusspirale vom Index —(n +2), wonach der Satz gilt: 
Die Sinusspirale vom Index » hat die Form der dem Zentralkraftgesetze y(&) =&E-"+” 


; ß kr z a 
entsprechenden Seilkurve; der zugehörige Parameter y= S ° der Seilkurve hat den Wert 
‚Io 
+1). 
(Gemäß Gl. (16) wird 
kim 
0 ” al Ta EN a BE er DE RE ER ER F AA 


es erfreut sich daher die Sinusspirale der Eigenschaft, daß ihr Krümmungshalbmesser ver- 


hältnisgleich der (1 —»)ten Potenz des Fahrstrahles & ist. Für die Seilspannung S dieser 
besonderen Seilkurven liefert Gl. (17) den Wert 


a N ee ee ee 5 ne 
die Seilspannung ist der —(1-+»)ten Potenz des Fahrstrahles £ verhältig. Bezeichnet o, 
den Krümmungshalbmesser dieser (£, 9) Kurven im Scheitel M, |wo 5, =1, &=0, Fr iem ); 
: 1 \ u 
so ist 0, =— — unabhängig vom Index », und da gemäß (19) 
? 
une”, 
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so folgt aus (20) z 
irv 


Y 9, \'-* 
s-s,(%) . 


Die Seilspannung S ist also für diese besonderen Seilkurven der e „ten Potenz des 
Krümmungshalbmessers verhältnisgleich. 

Der Kurvenklasse der Sinusspiralen gehören viele bekannte Kurven an, von denen einige 
mit Angabe ihres Index r, der zugehörigen Potenz n des für ihre Eigenschaft als Seilkurve 
maßgebenden Zentralkraftgesetzes und des Parameters y der Seilkurve in nachstehender 
Tafel zusammengestellt sind. 



































e * Seil- |», ‚Seil- | 
Para- Para- 
Index |} ER Potenz | span- Index Der Potenz | span- i 
u Kurve n nung meter naar Kurve a nung meter 
ss | ° S/S, 
1 Kreis —3 I —2 er ur 
# Lemniskate, ei 5) Gleichseitige 0 F1 1 
m v. Bernoulli > : Hyperbel u 
3 Kiepertsche Kurve —5 rl —4 — = = 
1 eh 5 8 3 1 3 1 
= Kardioide sr > "5 > Parabel 3 3 3 
1 ? 7 4 4 1 ei 2 5 2 2 
3 Cayleys Sextik 3 he 3 Trisektrix v. Catalan | — 3 Es 3 
D Fußpunktkurve 8 , 5 > |. Negative Fußpunkt- 4 e 1 
3 der Lemniskate nn L 3 3 kurve der gleichseit. 3 |F® 3 
; für den Mittelpunkt g ; 1 Hyperbel k 


Negative Parameter y deuten, da die Seilspannung S, nur positiv sein kann, an, daß 


zentrale Abstoßungskräfte (k,< 0) wirksam sind. Für den vorhin ausgeschiedenen Sonder- 


—_ 


fall n—=—1, dem das Kraftgesetz g(&)=- zugehört, liefert die Auswertung der Gleichungen 


> 
“ 
S 


(A—() die Formeln 


=+:yPl+yYaP—1, 


=- (lH eatlitrndt, 2222... 
‘ 

Sor/ 

S=$,1+,m9)=?7] & 


co 





Einer geschlossenen Integration wird dieser Sonderfall »—=—1 zugänglich, wenn die An- 
nahme konstanter Seilstärke aufgegeben und die Form jenes Seiles von veränderlicher 
Dicke gesucht wird, das überall gleich beansprucht wird (Seil gleicher Festigkeit). 

Mit o als konstanter Seilbeanspruchung und ö als veränderlicher Seildicke wird S=06, 


en a ' kö . a 
womit die Integration der Gl]. (2), in der P= - einzutragen ist, für die Diekenänderung des 


Seiles das Gesetz 


k 
r\o 
=.) 
liefert; aus Gl. (3), die in 
h Y® Es C 
vet 
06, (2° 


übergeht, folgt nach Integration 
142 k' 1+ - 
r cos (1 + ) o0=r, 


Der gesuchten Seilform entsprechen somit Sinusspiralen mit dem Index — 142) dieses ein- 


fache Ergebnis stammt mit etwas umständlicherer Herleitung von OÖ. Bonnet®). Im Zu- 


4) O. Bonnet: Liouvilles Journ. Bd. IX (1844), S. 97. 
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sammenhalte mit dem früher auf S. 235 bewiesenen Satze gilt demnach folgendes: Die Sinus- 


spirale vom Index » ist sowohl Seilkurve für Zentralkräfte ”+® bei konstanter Seil- 
dicke wie auch für Zentralkräfte ö&°' bei veränderlicher Dicke ö, wobei do=d,2"+»., 


S 
Im Bild 3 ist eine nach dem vorhin beschriebenen graphischen Verfahren ermittelte 
Schar von Seilkurven dargestellt, die dem Sonderfalle a» — — 1 mit den Parametern y= +3, 


+8, 1, + > zugehören. Für die gleichen Parameterwerte yundn=+ 1 sind die zugehörigen 
Seilkurven in Bild 4 eingetragen. Positiven Parametern y entsprechen Seilkurven, die sich 
ins Unendliche erstrecken, ihre Asymptoten gehen durch den Pol O0. Bei dem Parameter v=—l 
fällt (unabhängig vom Werte ») die Seilkurve mit dem Einheitskreis (Mittelpunkt im Pol 0) 
zusammen, wie aus der Anwendung des Krümmungskreisverfahrens ohne weiteres folgt. 





| / Bild 3 (links außen). 
| Ya Seilkurven für das Kraft 
| ü 1 

/ gesetz g(S) — 


Bild 4 (Mitte). 
Seilkurven für das Kraft- 
gesetz 9 (5)=5. 
t=+]) 


Polarochse 





M 





>, Bild 5 (unten). 


"We 





Mg r$ 4 


/ 
/E 
r | 7 
0 














OYaRzS 


Da das Wertepaar y„=+2,n=-+1 der Gl. (18) genügt, so ist die entsprechende be- 
sondere Seilkurve eine Sinusspirale vom Index —3 mit der Gleichung 


Pcos3o—=1, 


womit sich eine wünschenswerte Überprüfung der Zeichengenauigkeit ergab; die Fehler 
wachsen natürlich mit zunehmendem Winkel , bleiben aber durchwegs unter '/, v. H.; die 
Asymptote 0 A schließt mit dem Nullstrahl O M, einen Winkel von 30° ein (Bild 4). Bei 
Benutzung des Krümmungskreisverfahrens wurde mit dem Kontingenzwinkel der einzelnen 
Kreisbogenstücke nie unter 5° hinabgegangen. 

Weitere y-Kurven können aus den in den Bildern 3 und 4 dargestellten Scharen durch 
Interpolation konstruiert werden; eine Erhöhung der Interpolationsgenauigkeit läßt sich durch 
Anwendung des kürzlich von R. Sauer°) angegebenen Verfahrens erzielen. Die gezeichnete 
Kurvenschar wird vor der Interpolation in eine Schar von Kurven affın deformiert, welche 
drei unendlich benachbarte Punkte gemeinsam haben. Da die gezeichneten Seilkurvenscharen 
im Scheitel M, schon in zwei unendlich benachbarten Punkten übereinstimmen, so genügt 


: a u. ae 2 ß i nn ö 
hier die Transformation &=x#,n=ry (Bild 5) mit „= (wo ß eine beliebige Konstante, 
; . . 
y der Seilparameter), um die Gleichheit der Krümmung aller verzerrten Kurven im Scheitel M, 
} a ’ i 1 
herzustellen; sie beträgt ursprünglich, wie aus (16) bzw. (21) hervorgeht, u und nach 


zo 
der Verzerrung —vy, was aber bei der eben getroffenen Wahl für » gleich der beliebig wähl- 
baren Konstanten — ß ist. 





5) R. Sauer: Z. angew. Math. Mech. Bd. 20 (1940), S. 280 bis 284. 
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Wird die zu interpolierende einparametrige Kurvenschar durch rechtwinklige Koor- 
dinaten «,y in der Form 
:=&e(mY) y=y(p,7) 


dargestellt, so sind hierdurch außer den Kurven y=konst. (den Seilkurven) noch die Kurven 

p=konst. (die Polstrahlen durch 0) definiert, welche auf den y-Kurven nach @ bezifferte 

Skalen ausschneiden. Zwischen entsprechenden Skalenpunkten der g-Skalen der affın ver- 

zerrten y-Kurven wird sodann linear oder quadratisch interpoliert. Durch dieses Verfahren 

werden die Seilkurven mit ihren g-Skalen einander angenähert und daher genauere Ergeb- 
nisse bei der Interpolation erreicht. 

N) In Bild 6 sind die nach Vornahme der 

‚7? obigen Transformation der y-Scharen von Bild 4 

VD, 7 entstandenen, wesentlich schmäleren Scharen- 

KH. %  bereiche eingezeichnet, in denen zu interpo- 

L, lieren ist. Hierbei wurde bei positivem 

Y7 Parameter » die Konstante / gleich „1* ge- 

VG; wählt, so daß die dem Parameter y=1 ent- 

PILZ sprechende Seilkurve wegen »=1 unverzerrt 

bleibt. Auch bei negativem Parameter » wurde 

P=1 gesetzt, so daß die dem Parametery=—1 

zugehörige Seilkurve (Kreis) wegen v=—1 

unverzerrt bleibt. Dem Grenzfalle „=0, dem 

die durch M, gezogene Gerade LOM, als 

Seilkurve zugehört, entspricht alstransformierte 

9 Kurve je eine Parabel (mit dem Parameter 

EN p=]1 im positiven und negativen y-Bereich). 


Polorochse 





3 SS a N Die in Bild 4 gezeichnete Schar von Seil- 

IN kurven füra—=-+1 und verschiedene Parameter 

y ermöglicht z. B. die Lösung der folgenden 

Aufgabe: Auf einer um einen festen Punkt 

NN drehbaren waagrechten Scheibe ruhe span- 

\ N\N nungslos ein in den gleichweit von O entfern- 

N N ten Punkten A,B festgehaltener Faden von 

\ der Länge !; nach welcher Form wird er bei 

HN NO Drehung der Scheibe mit konstanter Dreh- 
Bild ö. schnelle » gespannt? 





E 
Ir 
= 





Die den Faden belastenden Fliehkräfte sind mit « als Masse der Längeneinheit des 
Fadens gegeben durch 
P= ur wo’ = (ur, °) B* 


der Vergleich mit (10) zeigt, daß 


k,k,=-ur®®, 
ga)=FE, ae: 
somit ist n=+1 und nach (11) der Parameter = r ” >, Wenn die Fadenlänge ! gerade 
gleich ist der Länge !, des Kreisbogens AB mit dem Mittelpunkte 0 (Öffnungswinkel 2a), 
dann ist die Seilkurve ein Kreis (y=—1) und S,=ur; w’. Für 1 >1, ist die entsprechende 


Form des Fadens durch den zugehörigen Parameter y, der Seilkurve bestimmt. Durch die 
Schnittpunkte der Polstrahlen 9=+ a mit der Schar der y-Kurven (n=+1) werden auf 
diesen Kurven bestimmte Längen abgegrenzt und es ist nun durch Versuch jene y-Kurve 


a ’ I : : . ; 
auszuwählen, die der gegebenen reduzierten Bogenlänge Er entspricht, wobei der Wert r, je- 


o 





0A a 
- bestimmt ist. 


weils durch — 
»1 
4. Seilkurven für das Zentralkraftgesetz P (&,9)=k, y(&,y). 

Die in Ziffer 2 abgeleiteten Grundgleichungen (A — C) behalten ihre Gültigkeit, doch ist das 
darin vorkommende \g(& y)d£ nicht mehr explizit angebbar, da der Zusammenhang &=£(y) 
vorerst noch unbekannt ist. Die bei der Durchführung des Krümmungskreisverfahrens not- 
wendige oftmalige Berechnung des Krümmungshalbmessers o nach Formel (14) wird verein- 
facht, wenn diese Formel unter Ausnutzung der Beziehung (13) in die Form 
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gebracht wird, wofür auch einfach 


oe= e yg(&,g)- MN . ° . . . . . . . . . . . . (22) 
gesetzt werden kann; hierbei wird der Punkt N auf dem Fahrstrahle OM (Bild 7) durch den 
Linienzug ON’ N gewonnen (ON ıL MM’, N N_LOM). 











by 
"7 N t 
u  ° N 
up = ei — 5 
Pi? u MM 
\ H 
£ | r2 Lu 7 \ N 
Ren 9 € x \ UL >Sx Polorachse 
Bild 7. E20 u - 
f u, ww M Ser A 
A \ F u] M=7% 
E 2 er 
M' 2 F ei 
I F Mm a A 
y,. 4 : u 
, [63 A26 
e: HM 3 Ye Bild 9 
| x Polachse 
; il 
gr 
8 [a2e2.2] 


Bild 8 


5. Ergänzung der Methode von Neuendorff. 

Abschließend sei noch folgendes bemerkt: Für die Konstruktion der Seilkurven erwies 
sich das Krümmungskreisverfahren in Polarkoordinaten nach R. Neuendorff zweckmäßig, weil 
die für das Problem charakteristische Differentialgleichung (13) die Darstellung einer sehr 
bequemen Formel für den Krümmungshalbmesser o erlaubt. Wird aber dieses Verfahren zur 
graphischen Integration einer Differentialgleichung von der Form 


&"=fl(p,8) 
angewendet, so macht die oftmalige Berechnung des Krümmungshalbmessers o nach der all- 
gemeinen Formel (14) das Verfahren ziemlich umständlich. Es läßt sich aber, wie ich a. a. O.*) 
bewiesen habe, der Krümmungsmittelpunkt der ebenen Kurve £=£(g) aus den Angaben £, 
und &’ in denkbar einfacher Weise konstruieren. Ist (Bild 88 OM=: und wird <= OM’ 
auf der Normalen zu OM aufgetragen, so gibt MM’ die Normale der Kurve € im Punkte M: 
aus dem gleichen Grunde ist M’M” die Normale der polaren Differentialkurve €’ in M’, 
wenn man OM”=" auf der Senkrechten zu OM’ aufträgt. Wird in M’ die Senkrechte zu 
MM’ errichtet bis zum Schnitte mit OM, und dieser Schnittpunkt B mit dem Mittelpunkte H 
von M’ M” verbunden, so schneidet diese Gerade BH die Normale MM’ im gesuchten 


Krümmungsmittelpunkte 2 der Kurve (' an der Stelle M’). Mit dem Ziehen zweier Geraden 





%, K. Federhofer: Sitz.-Ber. Akad. Wiss. Wien, Bd. 135 (1926), S. 79. 

”) Der für die Richtigkeit dieser Konstruktion a. a. O. (Anm. 6) auf Grund einer einfachen kinematischen 
Betrachtung geführte Beweis kann durch folgenden elementar-geometrischen Beweis ersetzt werden. Sei (Bild 9 
ds= Mu Mı ein Boge ıelement der Polarkurve (,do= mM. Mı das entsprechende Bogenelement der zugehörigen polaren 
Differentialkurve C”, «, v die Winkel der Kurventangenten t, t’ mit ihren Polstrahlen, — im Bild 9 soll daher ® deu 
Winkel zwischen f und O Mo (nieht B M,) angeben — so folgt mit dr als Kontingenzwinkel des Bogenelementes ds aus 
dem Dreiecke Mı M, Q: do ar 
"wo  sinw+v) 
und wegen ds—=QOMeo-dr=odr alurre 

de MR 
ds osinwm+e)' 
Es ist aber 
ds=dg] 72 +r2=dy-MoMo, 


do=dy Yr?+rr2—dg:MoMo, 

somit gilt Da ai 
mM Ma 

mM 3 osin(u+v) 


Zieht man durch M% die Parallele zu V Q bis zum Sehnitt G mit OÖ Mo, dann ist 
gie 
M0 En. 


rG 
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läßt sich also die umständliche Auswertung der Formel (14) vermeiden, wodurch eine nicht 
unerhebliche Vereinfachung des Krümmungskreisverfahrens in Polarkoordinaten erzielt wird. 
Wenn o sehr groß wird, also 2 weitab zu liegen kommt, so ersetzt man das zugehörige 
Bogenelement durch eine kurze Strecke auf der Kurventangente. 

Wenn die Punkte M der Kurve € außerhalb der Zeichenfläche zu liegen kommen oder 
aber sehr nahe an den Pol O heranrücken, ist es zweckmäßig, anstatt &(p) die Funktion 
n(g)=E(p) Fa einzuführen, wo a passend gewählt wird. Man verkürzt oder verlängert also 
jeden Fahrstrahl £ um ein konstantes Zwischenstück a, d. h. man konstruiert anstatt der 
Integralkurve £(g) die durch a bestimmte zugehörige Konchoide. Da 7 (p)= (pn ()=E"'gp),... 
so bedeutet der Übergang zur Konchoide keine Erschwernis der Konstruktion der Krümmungs- 
mittelpunkte. (Vgl. das in 6 behandelte Beispiel.) 

In analog einfacher Weise können auch die Krümmungsmittelpunkte der weiteren aus ( 
abgeleiteten polaren Differentialkurven (der Polarbilder 0” von &”, C”” von &” u. s. f.) kon- 
struiert werden. 

Bei der von R. Neuendorff*) gezeigten zeichnerischen Lösung von gewöhnlichen 
Ditferentialgleichungen beliebiger Ordnung in Polarkoordinaten werden die Krümmungshalb- 
messer der aufeinanderfolgenden Polarbilder €, C’, €”... im Wege der Rechnung bestimmt; 
die eben beschriebene einfache Konstruktion der Krümmungsmittelpunkte 2,02’, 2”... macht 
diese Zwischenrechnungen überflüssig und verwandelt das Neuendorffsche Verfahren in 
ein rein graphisches Verfahren. Die Polarkoordinatendarstellung (r,g) einer Funktion ist 
für die Zwecke der graphischen Analysis deshalb vorteilhaft, weil dabei — im Gegensatze 
zur Darstellung der Funktion in rechtwinkligen Koordinaten — alle Ableitungen von r nach 
p als Längen (Polarsubnormalen der aufeinanderfolgenden Polarbilder €’, C”,...) einfach 
konstruiert und nach obiger Konstruktion auch die Krümmungsmittelpunkte der Bilder €, €’, C”... 
bequem ermittelt werden können. 


6. Graphische Ausführung der Integration. 

Da die in Bild 8 gezeigte lineare Konstruktion des Krümmungsmittelpunktes (2 des 
Polarbildes © umkehrbar ist, so ermöglicht sie auch die schrittweise Zeichnung der einer 
zeichnerisch vorgegebenen Kurve (’ zugehö- 
rigen Integralkurve (', also des Polarbildes der 
Funktion &=\®d@. Bild 10 zeigt die an- 
zuwendende Konstruktion. Die gegebene Kurve 
C’ wird in einzelne Stücke M}, M/, M! M},,... 
unterteilt, dem gegebenen Anfangswert £,obigen 
Integrals entspricht der Anfangspunkt M, der 


der Senkrechten zu OM,;, liegen muß. Die 
Gerade M/ B,_L M/} M, und die Normale M} M}’ 
zur Kurve (” liefern den Krümmungsmittel- 
punkt £,. Beschreibt man um 2, einen kleinen 
Kreisbogen M,M, bis zum Schnitte M, mit 
der in O0 auf OM]} errichteten Senkrechten, 
dann ist M, näherungsweise ein Punkt der 
Integralkurve ©. Ist nach n solchen Schritten 


der Punkt M„ erhalten worden, so gibt die 
In 
Strecke OM,„ den gesuchten Wert\" do. 
% —=(0 
In Bild 10 ist nach diesem Verfahren 








y —nj4 
das bestimmte Integral £=\lncospdg unter 
y—=l 
Bild 10. der Annahme £=0 für 9=0 graphisch er- 
und daher gemäß (a) 
MOMs .. rG 
sın (u-+ v), = ’ 
MoMo M, 


Mo Mo 
“ \ a a Zu z 
Da man aber dem Dreiecke V MoMo die Beziehung 


_ 
Fur M; M; «sin (u+v) 
Ki A sin u 


v5 = M Mo -sin (u+v) 


— sin u: > 
sin u 


oder wegen 


ee, u re —, 
entnimmt, so ist VG=V M, ‚oder &H=H Mo ‚d.h. H halbiert die Strecke Mo Mo - 
») R. Neuendorff: Z. augew. Math. Mech. Bd. 3 (1923), S. 34. 


zu konstruierenden Integralkurve C, der auf 
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x 


mittelt worden. Die exakte Integration liefert ") 


= 
\" ( a 
.) ‘ m“ 
2 u (Zn+1) 
u 


“ T I 
£e=——In2+ 
4 


Nach M. Bresse '’) ist 


Sc yr 
Zn an +1) 


1) 


0: 915966, 


womit & = 0.086413 . 





’ 


Die gegebene Kurve C’ mit der Polargleichung ®—=1Incosg ist in Bild 10 durch die 
Punkte 0,1’,....9 (Maßstabseinheit gleich 400 mm) dargestellt, welche den um je 5° fort- 
schreitenden, von der Polachse 0 4’ aus im Gegenzeigersinn aufgetragenen Polwinkeln  ent- 
sprechen. Wegen E=0 für 9=0 fällt der Anfangspunkt M, der Integralkurve € in den 
Pol 0, so daß die oben angegebene Konstruktion versagt. Man führt daher die Konstruktion 
an der zu C’ gehörigen Konchoide © (Zwischenstück a gleich 50 mm gewählt) durch, wobei 
aufeinanderfolgend die Punkte 1, 2,....7 erhalten wurden: für den Zwischenpunkt M/; sind 
die erforderlichen Konstruktionslinien voll ausgezogen. 


Vom Punkt M,, ab wurde die Integration an der (’-Kurve fortgesetzt; hierbei ist zu 


350 350 


beachten, daß die Integration der Konchoide den Wert „= \""dpo=\F"dgy-+at(gp) liefert, 


0 v0 
so daß die Kürzung der Strecke OM; um das Maß (a 9)? (gleich 30,5 mm) den Punkt M, 
ergibt. Die Strecke OM, beträgt 34 mm, mithin ist der gesuchte Wert £ obigen Integrals 
bei Beachtung des Vorzeichens der &-Werte gleich — 0'085 mit einer Abweichung von rd. 
1,2 v. H. gegenüber dem oben angegebenen Wert. 





Dem Differentiationsprozeß entspricht das Zeichnen der aufeinanderfolgenden Subnormal- 
kurven, die Umkehrung des Verfahrens führt von der Funktion zu ihrem Integral. 

R. Grammel'*) bereicherte kürzlich die graphische Analysis durch ein sehr bemerkens- 
wertes Verfahren, wobei dem Funktionselement (£, 9) nicht der Punkt mit «en Koordinaten 


u s , . 1 + A u 
(£, p), sondern mit den Koordinaten (+) zugeordnet wird. Die Benutzung der Kehrwerte 


. 


£ 
von £& bietet den besonderen Vorteil, daß sich die Konstruktion oder Berechnung der 
Krümmungshalbmesser völlig erübrigt und die Polarbilder nach dem einfachen Tangenten- 
oder Sehnenverfahren konstruiert werden können. Das Grammelsche Verfahren bildet das 
duale Gegenstück zu dem bekannten Verfahren von E. Meißner '*). 262 


9 Vgl. W. Läska: Sammlung von Formeln der reinen u. angewandten Mathematik. Braunschweig 1888 bis 1894, 
S. 260, Formel (7). 


10) M. Bresse: Comptes rendus Bd. LXIV (1867), S. 1139. Der für obige Reihe bei W. Läska (8. 30, Formel 6) 
angebene Summenwert ist unrichtig. 


11) R. Grammel: Ing.-Arch. Band 10 (1939), S. 395. 


12) E. Meißner: Graphische Analysis vermittels des Linienbildes einer Funktion. Zürich 1932. 
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Anwendung der Nationalbuchungsmaschine für 
wissenschaftliche Rechnungen. 
Von Siegfried Kerridge in Stuttgart. 


Bei der Neuberechnung umfangreicher und vielseitiger Funktionentafeln 
treten gewisse Operationen — Bildung von Differenzenfolgen und Summie- 
rung von Differenzen und vor allem die Unterteilung der von den direkt 
berechneten Funktionswerten gebildeten Intervalle (Untertafelung) — massen- 
weise auf, deren Ausführung mit Hilfe der gewöhnlichen Rechenmaschinen 
äußerst mühsam ist. Im folgenden soll daher über ein Verfahren berichtet 
werden, das L. J. Comrie ausgearbeitet hat, um diese Operationen mit 
Hilfe der Nationalbuchungsmaschine (Klasse 3000) möglichst zu verein- 
fachen, und das er in dem heute schwer zugänglichen „Supplement to the 
Journal of the Royal Statistical Society IIT, No. 2, S, 88 — 113 (1936)“ ver- 
öffentlicht hat. 


Beschreibung der Maschine. Die Nationalbuchungsmaschine nach Patenten des Amerikaners 
Ellis ist eine Addiermaschine mit einer 12-spaltigen Tastatur, 6 Addierwerken, Druckwerk 
und einem Wagen ''). Zwei der Zählwerke (Nr. 1 und 3) können ebensogut subtrahieren wie 
addieren. Die Haupteigenschaft der Maschine, auf der ihre wissenschaftliche Nützlichkeit 
beruht, ist, daß eine mit der Tastatur eingestellte Zahl in jedes beliebige Zählwerk (und zwar 
bei Nr. 1 und 3 positiv oder negativ) oder unmittelbar in jede beliebige Kombination von 
Zählwerken gebracht werden kann. Während dieses Prozesses wird eine solche Zahl auto- 
matisch gedruckt. Der Inhalt jedes Zählwerkes kann in jedes beliebige andere Zählwerk 
(und zwar wieder positiv oder negativ in Nr. 1 und 3) oder in jede beliebige Kombination 
der anderen Zählwerke übertragen und zugleich gedruckt werden mit oder ohne Nullstellung 
des Zählwerks (im Folgenden wie auf den Tasten als „totalisiert“ oder „subtotalisiert“ bezeichnet). 


I. Summierung von Differenzen. 


Die zwei Hauptanwendungen der Maschine sind Summierung endlicher Differenzen 
und ihre Umkehrung: Bildung des Differenzengitters. Wir benützen die folgende Bezeichnung: 


| Il III IV V VI 
F- 3 ; 
) 0} f} 
= ”„ 
f 2 l j 2 UZ 
‘ -1,5 „ er IV 
F- 1 1 | | l : 
[ (\ . 
\ - 0,5 ai 0.5 IV — (1, vI 
fv N l,, u 0 V Ay 
R lo; Mr. ln,; jIV IN 
fi ’ I, m Bi, 
z 1; pr I; 
f: i 1, 
1, > 


Es werde angenommen: Bekannt seien die Funktionswerte f-3.... fr und ihre zuge- 
. . o * * . ” . 
hörigen Differenzen; verlangt sei, weitere Funktionswerte aus einer Folge von sechsten 
Differenzen zu berechnen. Wir wenden nacheinander die folgenden Gleichungen an: 


Ws =Alos +4"  =4ı +4 
= +4 A;=4s +4 
As =4As +4 h =fhı +4; 
Bei der Anwendung der Maschine enthalten zunächst die 6 Zählwerke die Werte von 
f> und der bei f2 beginnenden 5 Rück wärtsdifferenzen (415: Ai er A: Pig wird getas- 


tet und in dem Ns enthaltenden Zählwerk addiert; es ergibt sich dort Aps5- Dies wird in 


dem AlY enthaltenden Zählwerk addiert; so entsteht dort A!Y usw., bis sich fa in dem die 


1) Vgl. K. Lenz: Die Rechen- und Buchungsmaschinen. Leipzig 1932, S. 31, 85, 87, 92. 
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Funktionswerte enthaltenden Zählwerk ergibt: die anderen 5 Zählwerke enthalten dann die 
5 bei fs beginnenden Rückwärtsdifferenzen. Alle Einzelheiten des Rechnungsganges kann 
man der Tafel I entnehmen. 








Tafel I. Summierung der Differenzen von der sechsten an. 
Wagen- } Gedruckt Inhalt von Zählwerk ?) 
u Anschlag Operation i 
Stellung wird 1 3 4 2 6 5 
Anfangs Inhalt: ıY 05 er E. l\ A ! i 
1 Nicht add Taste x r # 
2 + Umschalter Taste en y n - r 
3 3 Subtotal 1 NY alN 
j ! Subtotal 3 „ ; : 
> 2 Subtotal 4 E A, = 
6 6 Subtotal 2 1, 1, n 
T > Subtotal 6 = x . f; 
v. N ’ - V IV vr " ’ 
5 Nicht add Subtotal 5 PA lo; | 4 15 lo I) 5 B 
u) tückkehr u 


Der Wagen kehrt in Stellung 1 zurück 


Als einfaches Beispiel betrachte man die Entwicklung von x® für ganzzahlige Werte 


von x mit Hilfe der konstanten sechsten Differenzen von der Größe 6!= 720. Die übliche 
Anordnung des Differenzenschemas ist: 
= 26 nl u jım jIv jv vi 
| 1 
63 
2 64 602 
665 2100 
3 129 2702 3360 
3367 5460 2520 
4 40% s162 5880 720 
11529 11340 3240 
65) 15625 19502 9120 1720 
31031 20460 3960 
6 46656 39962 13050 
70993 33540 
” 117649 73502 
144495 
Ss 262144 


Wenn man die Werte bis £=6 kennt, läßt sich die Folge der fettgedruckten Rück- 
wärtsdifferenzen leicht finden. Mit Stellung 2 beginnend, bringt man die Werte 25%, 5880, 
....., 46656 mit Hilfe der Tasten in die einzelnen Zählwerke. 4 Arbeitsgänge nach Tafel I 
liefern dann: 


x 4V1 AV IV 4” N” N = 
T 720 3240 9120 20460 39962 10993 117649 
8 120 3960 13080 33540 173502 144495 262144 
3 720 4680 17760 51300 124802 269297 53l44l 

10 720 5400 23160 74460 199262 468559 1000000 


Das Differenzenschema erscheint also hier in einer etwas modifizierten Gestalt; die hier 
in einer Zeile stehenden Zahlen sind als die Rück wärtsdifferenzen, je um '/, Zeilenabstand gegenein- 
ander verschoben, zu denken. 

Jeder Arbeitsgang braucht etwa 10 Sekunden; mit Berücksichtigung der Unterbrechun- 
gen für Auswechselung des Papiers usw. können 200 bis 300 Werte in einer Stunde fertig- 
gestellt werden. Negative Differenzen werden als dekadische Ergänzungen gedruckt, was 
auf den ersten Blick als unbequem erscheinen könnte, sich aber in der Praxis jedem anderen 
Verfahren überlegen gezeigt hat. Wenn 4V! das Vorzeichen wechselt, wird in Stellung 2 
der Umschaltanschlag umgelegt (was in wenigen Sekunden geschehen kann), so daß also 
diese Differenz stets in ihrer tatsächlichen Gestalt auftritt, folglich niemals dekadische 
Ergänzungen getastet zu werden brauchen. 


?2) Die Zählwerke sind in der Maschine in dieser unregelmäßigen Weise angeordnet. 


16* 








z 2 ee R n: s oz Z. angew. Math. Mech. 
244 Kerridze. Nationalbuchungsmaschine für wissenschaftliche Rechnungen Bd.21 Nr.4 Aug. 1941 





Il. Differenzenbildung. 

Es ist bekannt, daß eine Reihe von Werten einer Funktion bei gleichen Schritten der 
unabhängigen Veränderlichen auf zufällige Fehler geprüft werden kann durch Bildung suk- 
zessiver Ordnungen von Differenzen, bis die Differenzen verschwinden oder vielmehr im 
Vorzeichen schwanken, und durch Prüfung der Größe dieser Schwankungen. Sind dann die 
Fehler der Funktionswerte kleiner als eine halbe Einheit der 5. Dezimale, so wird die brauch- 
bare Schwankungsgrenze für verschwindende Differenzen annäherungsweise 

Verschwindende Differenz Im IV V VI i 


Grenze +3 +6 + 12 + 22 
Diese werden nur in seltenen Fällen überschritten, wenn nämlich mehrere aufeinanderfol- 
gende Werte in entgegengesetztem Sinne ungefähr um eine halbe Einheit abgerundet werden. 

Ein Fehler von der Größe E in der Funktion beeinflußt » + 1 aufeinanderfolgende 
Werte der nten Differenz um die Beträge EX dem Binomialkoeffizienten von (a — 5b)". So 
beeinflußt ein Fehler von +1 in einem Funktionswert sechs fünfte’ Differenzen mit den 
Beträgen +1, —5, +10, Id, +5, —1. Hiernach können Fehler leicht nach Lage, 
Größe und Vorzeichen bestimmt werden. Die Aufgabe, Differenzen von Hand oder mit einer 
gewöhnlichen Rechenmaschine zu bilden, ist sehr mühsam, indem sie viel Schreibarbeit 
eine der Hauptfehlerquellen — und vielfaches Einsetzen erfordert; deshalb ist es nicht ver- 
wunderlich, daß bei der Prüfung oft ebenso viele Fehler gemacht werden wie bei der Berech- 
nung der Funktionswerte selbst. Überdies wird ein mittelmäßiger Rechner es oft versäumen, 
eventuelle Fehler in den vorderen Ziffern der Funktion aufzudecken. — Differenzenbildung ist auch 
notwendig als das erste Stadium der Interpolation oder Untertafelung. 

Die Nationalmaschine kann in einem Arbeitsgang die Differenzen bis zur fünften Ordnung 
bilden. Dabei wird eine ziemlich wenig gebräuchliche Darstellung der Differenzen benützt, da 
sie praktisch zu demselben Zyklus von Operationen führt wie bei der Summierung: 

IV _ gi 


u 
lo; 

[273 ‚ IV 

h,; | lo; 3: m } 

[23 ‚ IV 
» —(Aı +A5 + As ) 
’ f ”„ 77 IV 
ls ; | I; + 1, j ln; + m } 
. = R ’ [3 7 IV 
I: (f. T I; 5 I, + As + I } 

Die Anwendung dieser Gleichung kann aus Tafel II ersehen werden, die keiner Er- 
läuterung bedarf. Sobald fs getastet ist, wird jene Gleichung angewendet. Die noch übrig- 
bleibenden Operationen, die in Wirklichkeit Vorbereitungen auf die folgende Zeile sind, 
drucken die Zwischendifferenzen und wiederholen den Druck der Funktion. Wenn die 
Funktion negativ ist, wird sie negativ in das Zählwerk 1 eingesetzt, d. h. unter Benützung 
eines Umschaltanschlags und wird in Form der dekadischen Ergänzung am Ende der Zeile 
gedruckt. Bei einem Vorzeichenwechsel wird ein kleiner Knopf auf dem umlegbaren Um- 
schaltanschlag gedreht. 

Tafel II. Bildung der Differenzen bis zu den fünften. 











Ä Ge- Inhalt von Zählwerk 
Anschlag Operation druckt 
wird 1 3 t 2 6 5 
Anfangs-Inhalt: (+ p Er t . | er u ra \ a, 8 
l Nicht add Taste $ B5 
2 Umschalter Taste 7; l; I. 
(+ 1) 0,5 
3 3 Total 1 I, Br 
) 1+4 Subtotal 3 AlV 5 2 Fee 3 
> — 1-+2 Subtotal 4 E | hg ni iv) er e. N, 
D I+6 Subtotal2 4% (3 +4), +41Y) R R : I 
. er . . ’ ’ . „ ZZ IV = 
1 1+5 Subtotal 6 1) ; (5 +4 +rAst4) E e ” % ls 
S | Subtotal 5 hı A+4,+% +45 +4) I 5 5, | 
9 Rückkehr Der Wagen kehrt in Stellung 1 zurück 
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x 





Beispiel zur Differenzenbildung. 








R 7 fix) IV IV E r l’ (x) 

117 108022 999998 YAIISS 999796 1112 20655 891978 
118 86468 3 999991 999787 99 21554 91353? 
119 64237 1) 999991 IIITTS 677 29231 935763 
120 11552 S 999999 999777 154 22685 958448 
121 — 18643 999994 999993 999770 224 22909 981357 
122 + 4263 10 3 999773 999997 22906 1263 
123 26940 999998 1 999774 999771 22677 26940 
124 19166 3 | 999778 999549 22226 49 166 
125 70723 0 | 999782 999331 21557 170723 
126 91460 63 67 999849 999180 20737 91460 
127 111015 999722 99789 999638 998818 19555 111015 
128 129366 551 340 999978 998796 18351 129366 
129 146288 999457 999797 999775 998571 16922 146288 
130 161623 270 67 999842 998413 15335 161623 
131 175233 999953 20 999862 998275 13610 1752833 
132 186895 999995 15 999877 998152 11762 186995 
133 + 196808 7 22 999899 998051 9813 196808 


Wenn man von rechts nach links liest, hat man auf einer Horizontalreihe einen 
Funktionswert und die von diesem ausgehenden Rück wärtsdifferenzen bis zur fünften Ordnung. 
Wenn bei der Arbeit ein Rechenfehler, d. h. kein Einsetzfehler gemacht wird, stimmt der 
zweite gedruckte Funktionswert (entweder in der Zeile, in der der Fehler gemacht worden 
ist, oder in der folgenden) nicht mit dem zuerst gedruckten überein. Dies ist eine überaus 
nützliche Probe und hat den weiteren Vorteil, eine Reservekopie zu liefern, so «aß der eine 
Druck als Druckvorlage benützt werden kann, während der andere mit den Differenzen 
zurückbehalten wird. Der zweite Druck kann beliebig oft wiederholt werden. Durch Wieder- 
holung des ganzen Prozesses können Differenzen bis zur 10. Ordnung gefunden werden, indem 
man die schon erhaltenen fünften Differenzen einführt. 

In dem angegebenen Beispiel steckt ein in die Augen springender Fehler. Um ihn 
festzustellen, suchen wir zunächst die beiden aufeinanderfolgenden fünften Differenzen, 
deren Summe nahezu 0 ist, d. h. die auf den Zeilen 128 und 129. Diese sind die fünften 
Differenzen, die +10fach und — 10fach den Fehler enthalten, so daß «der ursprüngliche 
Fehler in Zeile 126 steht und ungefähr '/, von 551 oder 543, d. h. ungefähr 54 oder 55 
beträgt. Eine Vertauschung von 0 und 6 ist also sofort nahegelegt und stimmt überein mit 
der Richtung des Fehlers; d. h. AY auf Linie 126 ist zu groß, da der (positive) Funktions- 
wert in dieser Linie zu groß ist. Folglich muß der fragliche Funktionswert lauten: 91406. 

Man gewöhnt sich rasch an die Anordnung der Differenzen in dieser Form. Eine 
zentrale zweite Differenz, d. h. eine solche, die in der üblichen Anordnung sich auf gleicher 
Zeile mit dem Funktionswert befindet, steht hier eine Zeile tiefer; eine zentrale vierte 
Differenz steht 2 Zeilen tiefer usw. 

Wenn man bedenkt, daß Differenzen bis zur 5. Ordnung mit einer Geschwindigkeit von 
200 bis 300 in einer Stunde gebildet werden können, kann man sich einen Begriff von dem 
Wert der Maschine für ein Recheninstitut machen. Kein Rechner könnte die Differenzen mit 
der Hälfte dieser Geschwindigkeit schreiben, selbst wenn er sie kennte. 


III. Untertafelung’). 


Maschinen ınit mehreren Zählwerken eignen sich zur Untertafelung, d. h. zur systema- 
tischen Interpolation in feste Intervalle, insbesondere Fünftel und Zehntel. Methoden zur 
Untertafelung, bei denen Ausgangsdifferenzen der zu interpolierenden Werte für jedes Gitter- 
intervall berechnet werden und eine konstante Differenz höherer Ordnung benützt wird, bis 
der nächste Gitterwert erreicht ist, wurden seit dem Ende des 18. Jahrhunderts angewendet. 
Ihre Nachteile liegen in der großen Zahl von Hilfsstellen, in der Notwendigkeit, Ausgangs- 
differenzen bei sehr vielen Intervallen zu berechnen, in der Unsicherheit, die Gitterwerte 
genau wiederzugewinnen, und in der Leichtigkeit, mit der Fehler in den einzelnen Phasen 
der Rechnung unentdeckt bleiben können. Was diesen letzten Punkt anlangt, so sind 
Hayashis Tafeln voll von Beispielen für die Unzuverlässigkeit dieser Methode. 


3) Vgl. Z. angew. Math. Mech. Bd. 14 (1934), S. 333 bis 339. 
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Diese Fehler verschwinden, wenn Maschinen mit großer Additionskapazität und vielen 
Zählwerken zur Verfügung stehen, insbesondere, wenn sie ihre Ergebnisse selbst zu drucken 
vermögen. Eine ganz wesentliche Verbesserung ist jedoch erreicht, wenn man nicht die 
Eingangsdifferenzen in jedem Intervall berechnet, sondern die fehlenden oder „Brücken*“- 
Differenzen zwischen sukzessiven Folgen konstanter Differenzenwerte. Die Methoden, die 
für den Gebrauch der Nationalmaschine entwickelt worden sind, haben stets folgende Vorteile: 


1. Die Bildung interpolierter Werte durch Summierung höherer Differenzen verläuft 
ununterbrochen, wenn einmal die Ausgangsdifferenzen gefunden sind. 

2. Die Gitterwerte werden genau wiedergewonnen. 

3. Die Anzahl der Hilfsstellen (Schutzstellen) ist so klein als möglich. 

4. Die Brückendifferenzen höherer Ordnung werden auf möglichst einfache Art berechnet 


Die Methode wird erläutert bei Unterteilung des Gitterintervalls in Zehntel, 
wobei die vierten Differenzen kleiner als 1000 vorausgesetzt sind. Die Interpolationsformel 
von Everett‘') werde in der Form: 


Y=(41-m)fo+tnfı+ Eö Mi + Er Mi 
angenommen, wobei der „Rückwurf“ der vierten auf die zweiten Differenzen benützt wird: 
M"” = 4" — 0,184 AV. 


Dann sind die vollständigen Ausdrücke für die interpolierten Werte: 


= 10f, 


(= 09, 


+o1f, 


0,0285 M/ 


— 0,0165 MT 





fu. = 0,8 f, + 0,2f, — 0,0480 M) — 0,0320 MY 
fa = f, +0,3f, - 0,0595 M, — 0,0455 M” 
= 06 f, +04 fı — 0,0640 MT — 0,0560 MT 
fl; 05h, +05 fi — 0,0625 M, — 0,0625 1” 
ff, +06 f, — 0,0560 7, — 0,0640 17 
fu; = 0,31, + of, — 0,0455 Mi) — 0,0595 M” 
fs = 927, + 98f, — 0,0320 M,, — 0,0480 77 


fu.s — 0] 4 +09f, 0,0165 M' 
h»= 1,0f, 
f1=09f, +01R, 


f.=08fi +0,2f, 


0,0285 MT 


0,0285 MY — 0,0165 MY 


- 0,0480 M — 0,0320 M' 


Die vierten Differenzen öIV dieser Werte verschwinden außer bei 0,9, 1,0 und 1,1. 
Dort findet man durch Differenzenbildung 
IV IV 
Ö0.9 = DR — 


0,0165 (11Y —- 0,184 47) +0,1 0,184 a1V 


31%, = + 0,0340 (A!Y —- 0,184 AY') — 0,2: 0,184 AlV 
Falls hierin 0,184 4!V abgerundet ist, bedeutet 0,184 AV! die zweite Differenz dieser abgerun- 
deten Werte von 0,184 AIV, 

Sollen nun diese „Brückendifferenzen“ auf drei Dezimalen (in Einheiten der letzten 
Stelle der Gitterwerte) beschränkt bleiben, so muß M” stets gerade sein. Zu diesem Zweck 
wird 0,184 AIV stets so abgerundet, daß 0,184 41V gerade ist, falls 4” gerade ist, und umge- 
kehrt. 

3 4) Fr. A.Willers: 


Prakt. Analysis, Berlin 1928, S. 22.—G. Schulz: Formelsammlung, Berlin 1937, Sammlung 
Gösehen Bd. 1110, S. 77. 
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en Nennt man die so abgerundeten Werte DIV und ihre zweiten Differenzen DV!, dann 
en sind die Brückendifferenzen in Einheiten der dritten Dezimale 
lie 
. == --165(4 — Di )+1W0D =A 
lie 
le: en 
öin= 4 (Ai — Di )-20D'=B 
ıft 


Um die Berechnung dieser 4 und B zu erleichtern, ist eine besondere Hilfstafel mit 
den stets geraden Argumenten (41V — DV!) berechnet worden nebst den zugehörigen DIV, 
wobei DIV gleich dem auf die nächstgelegene ganze Zalıl abgerundeten Wert von 0,184 








et (AV — DV!) ist. Ein Stück aus dieser Tafel ist hier abgedruckt. (Tafel III). 
1], Hilfstafel. (Tafel III). 
ıel — - - 
IV — DVI DIV A B IV — DVI DIV A B 
100 18 150 200 150 28 325 500 
.d: 2 19 217 332 2 28 292 32 
4 19 184 264 1 28 259 364 
6 20 251 396 6 29 326 496 
8 20 218 328 S 29 293 128 
110 20 185 260 160 29 260 360 
2 21 252 392 2 30 327 192 
4 21 219 324 l 30 294 124 
6 21 186 256 6 31 361 556 
S 22 253 388 8 31 328 188 
120 22 220 320 170 31 295 120 
2 22 187 252 2 2 362 552 
| 2 254 384 ! 32 329 184 
6 23 221 316 6 32 296 116 
8 24 288 448 8 33 363 548 
130 24 255 380 180 33 330 180 
2 24 222 312 2 3 397 612 
4 25 289 444 ! 34 364 544 
6 25 256 376 6 34 331 176 
8 25 223 308 8 35 398 608 
140 26 290 440 190 35 365 540 
2 26 257 372 2 35 332 172 
4 26 224 304 ! 36 399 604 
6 27 291 136 6 36 366 536 
S 27 258 368 8 36 333 468 
150 28 325 500 200 37 100 600 


Wenn für einen bestimmten Wert von (A1!Y — DV!) der in der Hilfstafel angegebene Wert für 
DIV nieht mit dem im Differenzenschema der Gitterwerte auftretenden tatsächlichen Wert von DIV 
übereinstimmt, so sind die Angaben der Hilfstafel für | 4} und |B| durch Addition (Subtraktion) 
E von 100 bzw. 200 für jede Einheit zu korrigieren, um die DIV tatsächlich von dem in der 
Hilfstafel hierfür angegebenen Wert nach oben (unten) abweicht. So hat man z. B. 


AW — DVI DIV A B 
162 28 127 92 
162 29 227 292 
ın- 162 30 327 192 «In der Hilfstafel!) 
162 3l 427 692 
en 162 32 527 892 
ck 
ge- Dabei hat A stets das Zeichen von AIV, B das entgegengesetzte. 


Die Anfangsdifferenzen sind, wie man leicht nachrechnet, in Einheiten der dritten 
ung Zusatzdezimale: 
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Anfangs-d' — 100 4,5 — 285 M, - 165M) — 10 +10 fı -85M —165M! . . . A) 
N een en Se 3 MR 

. ZZ Be w ne V s 
Anfangs-ö'” = M,„+M, Ap,5 ne ee Be re 
worin 
N ne a ae ae EI er 
stets gerade ist. 
Tafel IV. 
z I (a) ’ AIV DIV DV psi A B 
N } + %: 4 + 
0 11.014 199 56 375 134 35 3 
10 43 041 285 58833 145 27 1 l 291 436 
20 45 127 204 61 436 157 28 Bi +1 226 296 
30 47 274559 64 196 167 30 E +1 261 356 
40 49 486 110 67123 180 33 a 0 330 480 
50 51 764 784 70 230 195 36 4 +1 399 604 
60 54 113 688 713 532 213 10 3 436 656 
70 56 536 124 77047 224 1 - +4 470 720 
80 59 035 607 80 786 247 16 4 1 508 768 
90 61 615 876 84 772 268 50 ’ 0 578 888 
100 64 280 917 9 026 288 54 


Beispiel für die Anordnung der Rechnung. 

Aus der gegebenen Tafel IV der Gitterwerte mit ihren Differenzen entnimmt man 
zunächst die Werte von JIV und DIV und geht damit in die Hilfstafel III ein (wobei man 
die beiden ersten Spalten als Argumente benützt); man wählt hier für DIV den nächstgelege- 
nen ungeraden oder geraden Wert, je nachdem 4” ungerade oder gerade ist. So findet 
man z. B. (unter Beachtung der Vorzeichenregel) 


ı1V D!\ +4 +RB 

(a) 134 unger. 25 289 444 

(b) 145 unger. 27 291 136 
(ce) 157 gerade 28 326 — 100. — 496 + 200 
d.h. 26 d.h. — 296 


(d) 167 gerade 30 361 — 100 896 + 200 
d. h. 261 d.h. — 356 


(e) 150 unger. 33 330 — 480 
(f) 195 gerade 36 399 604 


usw. Im Fall (d) gibt die Hilfstafel IIl für DIV statt des nächstgelegenen geraden Wertes 
25 den Wert 29; infolgedessen müssen nach der oben gegebenen Anweisung die Tafelwerte 
A=2%, B=-4% 

durch Subtraktion von bzw. 100 und (— 200) korrigiert werden; derselbe Fall tritt bei (e) auf. 
Zur Berechnung der Anfangsdifferenzen hat man nach (4) 

M, = 56375 — 25 = 56350 
> (NB! Ohne Schutzstellen!) 

M, zz 58333 27 = 588306 
und erhält so (nunmehr mit Schutzstellen) nach (1), (2), (3) 


Anfangs-d’ —= 200 132 326 


Anfangs-öd’ = 565 956 
Anfangs-ö’”’ = 2456 
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Die Einteilung dieser Zahlen in dreiziffrige Gruppen entspricht der Anordnung der Farben 
der Tasten. 


Wenn eine oder mehrere dieser Anfangsdifferenzen negativ werden, erscheinen sie 
bei der Berechnung mit den gegebenen Formeln in Form dekadischer Ergänzungen und 
werden so in die Maschine eingesetzt. 








Tafel V. 
N n ölV oe" ö” ö (x) 

0 41014 199 500 

1 200 132 326 41214 331 826 

v 3 565 956 200 698 282 41 415 030 108 
1 3 2456 568 412 201 266 694 41 616 296 802 
2 4 2 456 570 868 201 837 562 41 818 134 364 
3 5 2 456 573 324 202 410 886 42 020 545 250 
4 6 2 456 575 780 202 986 666 42 223 531 916 
5 T 2 456 578 236 203 564 902 42 427 096 818 
6 N 2 456 580 692 204 145 594 42 631 242 412 
7 HM) 2 456 583 148 204 728 742 42 835 971 154 
8 10 2456 585 604 205 314 346 43 041 285 500 
9 11 + 291 2747 588 351 205 902 697 43 247 188 197 
10 12 — 436 2311 590 662 206 493 359 43 453 681 556 
11 13 + 291 2 602 593 264 207 086 623 43 660 768 179 
12 14 2602 595 866 207 682 489 43 868 450 668 
13 15 2602 598 468 208 280 957 44 076 731 625 
14 16 2602 601 070 208 882 027 44 285 613 652 
15 17 2 602 603 672 209 485 699 44 495 099 351 
16 18 2602 606 274 210 091 973 44 705 191 324 
17 19 2602 608 876 210 700 849 44 915 892 173 
18 20 2 602 611478 211312 327 45 127 204 500 
14 21 + 226 2 828 614306 211 926 633 45 339 131 133 
20 22 - 296 2532 616 838 212543 471 45 551 674 604 
21 23 + 2236 2 758 619 596 213 163 067 15 764 837 671 
22 24 2 758 622 354 213 785 421 45 978 623 092 
23 25 2 758 625 112 214 410 533 46 193 033 625 
24 26 2758 627870 215 038 403 16 408 072 028 
25 27 2758 630 628 215 669 031 46 623 741 059 
26 28 2758 633 386 216 302 417 46 840 043 476 
27 29 2 758 636 144 216938 561 47 056 982 037 
28 30 2 758 638 902 217 577 463 47 274 559 500 


Tafel V zeigt das Ergebnis der Maschinenrechnung für die ersten 30 Werte. Dabei 
wurde der Kunstgriff von v. Wrede angewandt, die Gitterwerte (und damit auch die Zwischen- 
werte) um 5 Einheiten in der ersten Zusatzstelle zu vermehren, so daß sie dadurch automatisch 
zu derselben Stellenanzahl wie bei den Originalwerten abgerundet werden. Jeder Gitterwert 
wird sogleich nach seiner Berechnung geprüft; zeigt sich ein Fehler, so wird der Wert von 
ö’, der achtmal seit dem letzten Gitterwert erschienen ist, mit 4”” — DV verglichen; jeder 
Mangel an Übereinstimmung zeigt an, daß die benützten Werte für A und B nicht richtig 
sind. 


Die Rechnung geht mit einer Geschwindigkeit von 300 bis 400 Werten in der 
Stunde vor sich. Es sei noch einmal darauf hingewiesen, daß jede beliebige Anzahl von 
Kopien der letzten Spalte angefertigt werden kann. 


Jede Änderung in den Voraussetzungen (d. h. der Anzahl der Zwischenwerte in jedem 
Gitterintervall oder jede vernachlässigbare Differenz) gibt Gelegenheit, die beste, zur Be- 
arbeitung geeignetste Form der allgemeinen Methode anzuwenden. 300 
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KLEINE MITTEILUNGEN 


Instabilität laminarer Grenzschichten an 
konkaven Wänden gegenüber gewissen 
dreidimensionalen Störungen. 


Die folgende kleine Mitteilung gibt im wesent- 
lichen den Inhalt meines auf der Tagung der Ge- 
sellschaft für angewandte Mathematik und Mecha- 
nik im Mai 1940 in Braunschweig gehaltenen Vor- 
trags wieder. Eine ausführliche Darstellung ist in- 
zwischen in den „Nachrichten von der Gesellschaft 
der Wissenschaften zu Göttingen“ erschienen !), es 
ist mir jedoch von verschiedenen Seiten nahegelegt 


worden. die Ergebnisse in kurzer Zusammen- 
fassung einem weiteren Kreise zugänglich zu 
machen. 

Die Untersuchung der Frage des Wandkrüm- 


mungseinflusses auf die  Turbulenzentstehung. 
welcher im Hinblick auf die Formgebung von Strö- 
mungskörpern größeres Interesse entgegengebracht 
wird, hatte ich zunächst in Angriff genommen, in- 
lem ich in der bisher üblichen Weise vorgegebene 
zweidimensionale laminare Grenzschichten kleinen 
Störungen unterwarf in Gestalt zweidimensionaler 


Wellen. welche sich in Strömungsrichtung fort- 
pflanzen. Es zeigte sich, daß gegenüber ebenen 


Wänden eine schwache Wandkrümmung von außer- 
ordentlich geringem Einfluß auf die Stabilität der 
Grenzschichten gegenüber diesen Störungen bleibt. 
Bemerkenswerterweise ergab sich entgegen den 
Erwartungen ein stabilisierender Einfluß konkaver 
und ein anfachender Einfluß konvexer Wände ?). 
Dieser schwach stabilisierende Einfluß  zwei- 
(imensional-konkaver Wände gegenüber den ge- 
nannten zweidimensionalen Störungen bleibt jedoch 
bedeutungslos angesichts einer kräftigen Instabilität 


von Grenzschiehtströmungen an diesen Wänden 
und nur an solehen) gegenüber gewissen (drei- 


«limensionalen Störungen. Es handelt sich bei die- 
sen Störungen um Wirbel in der Grenzschicht. 
deren Achsen parallel zur Hauptströmung verlaufen 
(siehe Veranschaulichung in Bild 1). von ganz 
analoger Natur also wie jene Wirbel, die in der be- 
kannten Untersuchung von G. I. Taylor über die 
Stabilität der Flüssiekeitsbewegung zwischen rotie- 
renden Zylindern theoretisch berechnet und experi- 
mentell nachgewiesen wurden ?). 

Die Erschwerung der Rechnung durch den kom 
plizierteren Störungsansatz ließ sich bewältigen. 
dabei aber auch die Reibung vollständig berück- 
sichtigen, wobei es uns jedoch dem Stande der ex- 
perimentellen Untersuchungen entsprechend nur 
auf größenordnungsmäßige Aussagen über den EI- 
fekt ankam. Der Rechnungsgang werde im folgen 
den kurz zeschildert. 

Zunächst zu den Voraussetzungen der Theorie: 
Wie bei allen früheren Stabilitätsuntersuchungen 
mittels der Methode der kleinen Störungen werden 
aus Gründen der rechnerischen Durchführbarkeit 
Veränderungen der ungestörten zweidimensionalen 
Grundströmung in Strömungsriehtung — es sei dies 
die -Riehtung ( Bogenlänge der Wand) — ver- 
nachlässigt und es wird angenommen, daß sich diese 
Grundströmung in jener Zeit, in der die Störungen 
verfolgt werden sollen, nur sehr wenig ändert. Die 
Geschwindigkeit der Grundströmung — eine durch 


1), H. Görtler: Über eine dreidimensionale Instabilität 
laminarer Grenzschichten an konkaven Wänden. Nachr. 
Ges. Wiss. Göttingen, Math.-Phys. Kl. Neue Folge, Fach- 
gruppe I, Bd. 2 (1940), S. 1 bis %. 


®?, H. Görtler: Über den Einfluß der Wandkrümmung 
auf die Entstehung der Turbulenz. Z. angew. Math. Mech. 
Bd. 20 (1940), S. 138 bis 147. 


3) G. I. Taylor: Stability of a visecous liquid contained 
between two rotatingeylinders. Phil. Trans. Roy. Soc. Lon 
don, Series A, Bd. 225 (1923), S. 289 bis 343. 


irgendeine Vorgeschichte auf Grund der Zähie 
keitswirkung -zur Ausbildung gekommene Grenz 
schichtströmung werde somit durch 4 = u, (4) 
zegeben, wo Y der senkrechte Wandabstand be- 
deutet. Zu ihr gehört der ungestörte Druck p% (4). 
Was die Wandkrümmung I/R betrifft. so wird dies« 
der Einfachheit halber als konstant angenommen 
und als hinreichend klein vorausgesetzt. um in Po 
tenzentwicklungzen nach ö/R (ö Grenzschichtdicke) 
mit den linearen Gliedern abbrechen zu können. 
Es wird R an zur Strömung konkaven Wän- 
(len positiv, an zur Strömune konvexen Wänden 
negativ gesetzt. 
Unser Störungsansatz lautet nun 


u=w(y) + u (y)cosazert 


v»— v,(y) cosa2 edt 

5 1). 
vw w,(y)sinazet 
P=m(y) + pı(ly) eosazert 


Darin bezeichnen z die zur .r- und y-Achse senkrechte 

Koordinate (also in Richtung der Erzeugenden des 

Zylinderstücks. aus dem die Wand gebildet wird). 

! die Zeit, uw, v, w die Geschwindigkeitskomponen 

ten in Richtung der Koordinaten . y. 2. Die reell 
Ir 

Größe «a ist gleich 3 


anzenommene ‚wo 4 die 


Wellenlänge der Störung bedeutet. und auch für 
‚tie Größe £ liefert die Rechnung nur reelle Werte: 
das Vorzeichen von dentscheidet dann über zeitliche 
Anfachung oder Dämpfung der Störung. 

Auf dem Wege über die Navier-Stokesschen 
Gleichungen und die Kontinuitätsgleichung gelangt 
man zu einem System linearer Differentialgleichun- 
sen für u, und ”, allein: w, und ?, berechnen sieh 
dann aus va, und ”, durch einfache Differentiations 


prozesse. Es ergibt sich unter Verwendung iler 
Dimensionslosen 
Y r__ % 
= Fe T Zu (U,= U(Ö)). 
0 
und der Gröben 
U,8\"3 
u | ° u, v=v, 
v 


Differentialgleichungssystem 


(las 





?u ö°\ U 
Ar ut 
dn? -v dy 
| d? Id? »* 88% |] 
| < - az de] |" « — (a: +! —\v’ (2), 
do ||d n? v | i 
U,6\? ö 
- 2a? 8? WW. 
wrr | 5 R 
r R 1 dv ‘ 
und es ist ferner w, . An Randbedin- 
addn 
‚ ‚ dv 
eungen hat man zu fordern: x (0)=» (U) nr 7 (O0 
/ 


wegen des Haftens der Flüssigkeit an der 
Wand: ferner ist das Abklingen der Störung für 
"> zu verlangen. Insgesamt hat man das fol- 
gende Eigenwertproblem eines homogenen Ranrl- 
wertproblems zu lösen: Bei gegebenem U (n) und 
ö/R ist zu jeder Wellenlänge 4 der Störung und zu 
U,ö 


jeder Reynoldsschen Zahl die Anfachungs- 


eröße zu ermitteln. Für die neutralen Störungen 





XUM 





141 


17 
4) 


y). 
Drz 
en 
’ 

io) 


1], 


ei 


te 
Es 
dl). 
»n 

ell 


lie 
ür 
te: 

he 


en 
gt 
INn- 
en 
18 

ler 


(2), 


lin- 


der 
für 
fol- 
nıl- 
und 
zu 





Z. augew. Math. Mech. 
Bd. 21 Nr.4 Aug. 1941 


Kleine Mitteilungen 251 





Bild 2 (rechts). Gesamt verlauf der kritischen 
e Ur B-, i n 
Zahl z über @«* und die ersten An- 
3.42 
fachungskurven = eonst. Dureh die 
Vv 
parallelen Geraden werden in diesem Dia 
gramm einzelne Wirbel (A=econst.) in einer 
sich langsam verdickenden Grenzschicht bei 
konstanter äußerer Strömungsgeschwindig- 
keit U. dargestellt. 


Bild 1 (unten). Wirbelstörungen in der 
Strömungeiner Flüssigkeit aneinerkonkaven 
Wand. Die Wirbelachsen verlaufen parallel 
zur Hauptströmungsriehtung. 





A257 


(3=0) insbesondere hat man zu jeder Störungs- 
wellenlänge 4 eine .„kritische* Reynoldssche 
Zahl zu bestimmen. Man erkennt übrigens unmitte] 
bar aus (2), daß die Revnoldssche Zahl und 
die Wandkrümmung nur in Gestalt der Dimensions- 


U,ö #2 . 
losen . A eineehen (im Rahmen unserer 


Theorie der ersten Näherung). 

Der Gang der praktischen Durchrechnung ist nun 
der folgende: Das vorliegende Differentialglei- 
chungsproblem nebst Randbedingungen wird auf 
ein äquivalentes Integralgleichungssystem zurück- 
eeführt, sodann werden die auftretenden Integrale 
nach dem Vorbild der Fredholm schen Theorie 
lurch endliche Summen ersetzt, so daß die Eigen- 
wertberechnung sich auf die Auflösung einer Deter- 
minantengleichung reduziert. Auf die Eigenwert- 
berechnung allein kommt es uns an. 

Die Rechnung war in der geschilderten Weis« 
praktisch in der in erster Linie interessierenden 
Gegend des Minimums der kritischen Reynolds- 
schen Zahl über «ö mit hinreichender Genauigkeit 
durchführbar und ließ sich durch asymptotische 
Aussagen für eroße und kleine «6 ergänzen. Eine 

U,91/9 ’ 
Auftragung von über a» wobei # 

12 

die Impulsdieke der Geschwindigkeitsprofile be 
zeichnet — ergab ungefähre Übereinstimmung der 
Kurven für ganz verschiedene Geschwindigkeits 
profile gleicher Impulsdicke. Dabei wurde auch ein 
Profil mit Wendepunkt (Druckanstieg in Strö- 
mungsriehtung) durchgerechnet. In dieser Auftra- 
zung spielt also die jeweilige Gestalt des Geschwin- 
digkeitsprofils keine wesentliche Rolle. Die Rech- 
nung lieferte Instabilität nur für R>0O (konkave 
Wände). 


4, H. Blasius: Grenuzsehichten in Flüssigkeiten mit 


kleiner Reibung. Diss. Göttingen 1997, erschienen in Z. 


Math, und Physik Bd. 56 (1908), S. 1 bis 37. 


02 





- deren Berücksichtigung eine 


94 _06 08 1 


4 + 


ı i i | i J 
025 0% 07510 % 5@ 


Die Ergebnisse der Rechnung sind für «las 
Blasiussche Plattengrenzschichtprofil *) in Bild 2 
im doppeltlogarithmischen Netz aufgetragen. Die 
Kurve #/—0 gibt den Verlauf der kritischen Di- 
mensionslosen 2 | " über a # wieder. Darüber 
v R 
: . En bu" 
hinaus sind einige Anfachungskurven const. 
>0 angegeben, soweit die Rechnung einigermaßen 
zuverlässig blieb. 

Als Stabilitätsgrenze ergibt sich 4 ® 058. 
(Bei den Taylorschen Wirbeln zwischen rotieren- 
lem inneren und ruhendem äußeren Zylinder ergab 
sich hierfür der Wert 2,81.) Durch dieses Minimum 
werden jene Wirbel festgelegt. die mit wachsende: 
Reynoldsscher Zahl als Erste angefacht werden. 
Ihre Wellenlänge ist gegeben durch4= 45%: es 
sind damit Wirbel. welche Zellen von einer Breite 
in z-Richtune von ungefähr 25 erfüllen. Sie grei- 
fen in y-Richtung natürlich auch über die Grenz 
schicht hinaus. (Im obengenannten Taylorschen 
Fall ergaben sich zwischen den Zylindern etwa 
quadratische Zellen.) 

Es muß betont werden. daß es sich bei diesem 
ersten Auftreten von neutralen und bei etwas 
höherer Revnoldsscher Zahl angefachten Wir- 
beln durchausnicht umein Turbulentwerden der Strö- 
mung handelt. es wird vielmehr nach wie vor eine 
durchaus geordnete Strömung erfolgen. (Dabei ist 
natürlich der Fall ausgenommen. daß bei sehı 
eroßem R/öd bereits gewöhnliche Plattenturbulenz 
auftritt.) Erst bei wesentlich höheren Reynolds- 
schen Zahlen, bei welchen die Störungen eines 
sanzen Wellenlängenbereichs «eine hinreichende 
Anfachune erfahren. kann der Umschlag in turbu- 
iente Strömune erwartet werden. Ganz em 
sprechend verhielt es sich auch bei den Taylor- 
schen Beobachtungen an der Strömung zwischen 
rotierenden Zylindern. 

Die Ergebnisse der Theorie werden auch bei Zu- 
lassune einer schwachen Veränderlichkeit der 
Grundströmung in z-Richtung — rechnerisch würde 
wesentliche Er- 
schwerung bedeuten jedenfalls als gute Nähe- 
ıung herangezogen werden dürfen. Wir fragen 
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‚laher nach dem Schicksal eines Wirbels fest vor- 
segebener Wellenlänge A beim Wandern in Strö- 
mungsrichtung durch eine sich bei konstanter 
äußerer Strömungsgeschwindigkeit U, langsam ver- 


: : U,» 3 \ 
diekende Grenzschicht. In dem | . | me »} 
v 


Diagramm wird sein Weg durch eine Kurve der 
- U,01/» . ; 
Gestalt " — const. (a Ö)’!2 wiedergegeben. 
v : 
in der Auftragung des Bildes 2 also durch Geraden 
der für verschiedene feste 4-Werte eingezeichneten 
Art. Man kann demnach der Figur entnehmen. 
“ Be 1 Sy 
daß Wirbel mit — 
R 
den, woraus für jene Wirbel, die überhaupt vor- 


< 50 nie angefacht wer- 


übergehend angefacht werden, die Forderung 
7 AU es 
R? 13,6 > } folgt. Die Zahlenangabe bean- 


sprucht jedoch nur, als grobe Aussage bewertet 
zu werden. Wirbel, welche dieser Forderung ge- 
nügen, geraten beim Wandern durch die Grenz- 
schicht früher oder später in ein kritisches Ver- 
hältnis zu der sich verdickenden Grenzschicht und 
werden dann so lan&e angefacht, bis ein zweites 
kritisches Verhältnis eintritt. wonach sie wieder 
zedämpft werden — letzteres gilt jedoch nur, wenn 
die Störung auf dem zuvor durchsehrittenen An- 
fachungsweg nicht bereits den unsere Theorie be 
schränkenden Rahmen ..kleiner“ Störungen über 
schritten hat. 

Was nun noch «(ie im Hinblick auf den Umsehlag in 
turbulente Strömung gefährlichsten Wirbei betrifft. 
so kann man jedenfalls aussagen, daß es jene Wir- 
bel sind, welche sich zu Anfang ihres Anfachungs 
weges im Gebiet der Minima der Anfachungskurven 
39? 22 

eonst. halten. Für 


genauere Aussagen 


müßte die jeweilige Gesamtanfachung beim Durch- 
laufen desInstabilitätsbereiehsdurch Integration er- 
mittelt werden, allein hierfür reiehen unsere weni- 


zen bei größeren Werten von unzuverlässig 


werdenden Anfachungskurven nicht aus. 

Es ist in Bild 2 noch ein Eintrag gewisser Meß 
ergebnisse von M. und F.Ulauser?°) über das Auf- 
treten turbulenter Strömung an konkaven Wänden 
gemacht. Umgerechnet ergeben sich die gestrichelt 
U,d [ 0 

R 
und 8,6). Wagt man eine Extrapolation der An- 
fachungskurvenschar, so findet man, daß beim Um- 
schlag die Grenzschichtdicke, grob gesagt, bis zur 
Größenordnung der Breite der am stärksten ange- 
fachten Wirbel angewachsen wäre, welche ihrerseits 
U,R 
v 


angegebenen Werte von (nämlich 10,6 


9 
n he 
durch R 50 ) gegeben wären. 


H. Görtler. 2355 


Göttingen. 


5) M. u. F. Clauser: The efleet of eurvature on the 
transition from laminar to turbulent boundary layer. 
NACA Techn. Note No. 613, 1937. 


Zur Spannungserhöhung bei gelochten 
gezogenen Streifen. 

Ein Streifen von der Breite 1 mit kreisrundem 
oder elliptischem Loch in der Mitte werde in Längs- 
richtung gezogen, Bild 1. Im Querschnitte. deı 
durch das Loch geschwächt ist, tritt die durch- 
schnittliche Spannung 0, auf. Die wirkliche größte 
Spannung o,,., in den Punkten 1 und 2 ist ein 
Mehrfaches davon. Bei kleinen Löchern kann man 


das Verhältuis 0,,.,/0g näherungsweise angeben. - 


da man hierfür als Lösung die gelochte Vollebene 
nehmen kann. 








La » 4 
A M Ph / 
| ee: — 1 
a | | as r— I 





2.2 9 
Bild 1. Bild 2. 


Ich will für den anderen Grenzfall, daß die Loch 
breite fast gleich der Streifenbreite ist, eine ein- 
fache Lösung angeben. Nach Bild 2 sei in der 
engsten Stelle die nachgebliebene Streifenbreite auf 
jeder Lochseite gleich a und der Krümmungshalb- 
messer des Loches gleich o. Es wirke in diesem 
Querschnritte die Kraft P und das Biegemoment 
Mo. Im Abstande z ist die Breite y—= a+ 22/20 
und das Moment M— Ma — Px?/4o, da die Kraft 
P am Hebelarme 2?/4o wirkt. Für die Durch- 
biegung v erhält man nach der elementaren Biege- 
theorie, die für a>o anwendbar ist: 


E d’v_ M,— Pı2]4o 
Fri 1 r: : 
d „utzRo® 
% 
. ER ‚iG%y : 
Aus Symmetriegründen mub F = -d.r=0V sein. 
) 
Hieraus erhält man Mo = !/, Pa. 


Man erhält im engsten Querschnitte eine lineare 
Spannungsverteilung, die vom Werte null außen 
bis zum Werte o,,.x — 20,; innen ansteigt. 


Dresden. C. Weber. 236 


Potenzreihenzerlegung von tg .: und cig.. 


Ich benötigte die Potenzreihen dieser Funktionen; 
in den Handbüchern fand ich nur wenige Glieder. 
Die Reihen sinz und cosr gegenseitig zu teilen, 


sin 


oder Incos.r und In in Reihe zu entwickeln 
4 


und zu differentiieren, war mir zu umständlich. Am 
zünstigsten erwies sich folgender Weg: 


ER N ni 
\ k,z", etgr= + \ hc". 
H 7 






gg: = 
1.8,» 8 
etg® tgz—=2ctg2z 2( + N h,2" a"). 
r de 
1,3, 9 
Hieraus 
. =) 
a + y Ay En — N kn" 
3,5 1,3,» 
1 u 
ar re h 9An +1 r", 
T ae Rn, 
> 
Dieses gibt folgende Beziehung zwischen den 


Koeffizienten beider Reihen: 


= — (er 1)A,. 
Weiter gilt tgr tg x =1 
/ Ya z ER r 
| (2r +1 —1) h, ="), ++ Am m) ,; 


1,3,5 1,3,5 
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Fürrn+m=0: ,=—1/, 

für n+-m--s>0, s gerade, n und m ungerade: 
" v 

(+? —1)R4ı = " P. 2 N) Au hs 


n» 


Also für s=2, 4, 6 usw.: 


3.56, =— IM, 6 h=— ”. 
2 
5 2 2 
79h, = - (döhla +15, h,). A\=— 3.5.7 
15-174; (3h,l,+ 1573-4634, h,). h,= 93 n — 
. 9.4 
In gleicher Weise erhält man: 
> 
= — i 
: 3°.5-7-11 
RER 1382 
Bu 36.53.72.11-13° 
h=- - 


36.52.7-11-13° 


Besteht der Koeffizient von A,;, aus genügend 
hohen Primzahlen, so muß sich der Zähler der 
rechten Seite der Gleichungen durch diese teilen 
lassen (Beispiele: A,, Ag, As. Gegenbeispiel A,,). 
Dieses dient zur Prüfung der Zahlenrechnung. 

Dresden. C. Weber. 258 


Bemerkung zu einem Versuch von Mach. 

Wohlbekannt ist der folgende Versuch von 
Mach!) zur Erläuterung des Schwerpunktsatzes: 
Auf einen kleinen Wagen ist ein Pagemotor ge- 
lagert. Die hin- und hergehenden Teile des Mo- 
tors bewirken. da der Gesamtschwerpunkt in Ruhe 
bleibt, daß der Wagen Bewegungen ausführt. Diese 
Versuchsanordnung wurde von Professor Dr. 
Lechner benutzt, um die Einwirkung der be- 
wegten Teile auf das Fundament zu zeigen. Hält 
man den Wagen fest, so spürt die Hand die 
Kräfte, welche das Fundament aufzunehmen hat. 
Dabei zeigt es sich. daß die Drehzahl des Page- 
motors abnimmt. Im folgenden soll versucht wer- 
den zu zeigen, ob ein Drehzahlunterschied bei frei 
beweglichem und festzebremstem Fundament vor- 
handen sein kann. 

Vorausgesetzt sei, daß die Kurbelstangenlänge 
des Pagemotors unendlich lang sei, d. h. daß eine 
Kurbelschleife gegeben sei. dann ist. wenn z die 
Koordinate im beweglichen System und « den 
Kurbelwinkel bedeuten 


z=rcosp und F=—rwsing 
3edeutet © die Wagengeschwindigkeit. so wird 


die Absolutzeschwindiekeit der hin- und hergehen- 
den Teile 


n=V—Tr@Sing 


und die kinetische Energie L 


L= 5 m, (v—rosinp)”— : mMsv?—+ - I. 
Darin bedeutet m, die Masse der bewegten Teile 
des Pagemotors, m, die mit der Geschwindigkeit © 
sich bewegenden Teile und J, das Trägheitsmo- 
ment der rotierenden Massen. Wendet man den 
Schwerpunktsatz an, so wird 


1) E. Mach: Die Mechanik in ihrer Entwicklung, 
2. Aufl., Leipzig 1912, S. 229. 


m, 


m, + m, 


0) sin q : 
damit wird 


L= E 


2 
oe. 


1 mim, \ 1 
=> 1 2 ©?’ sin?g u = J; 

2 m, + m, 2 

Die Anwendung der Lazrangeschen Glei 
chung 


(+) _dL_ 


oO: 
4; ei 


°gi 
ergibt bei Unterdrückung einiger Zwischenrech- 
nungen 

m, My Ei fi 


1 5 
n @o— > ©? sin? G 1-J 
Mm, TMa\ - e 


9) ' P : °: 


Das Moment M steht nun mit der Radialkraft P, 
in einer einfachen Beziehung. es ist 


P,= = eotgY. 


Zur Berechnung der Radialkraft kann G. (1) 
verwendet werden, dazu sei vorausgesetzt, daß das 
Schwungrad so bemessen ist. daß die auftreten- 
den Geschwindigkeitsänderungen klein sind, man 
annehmen darf, daß die Winkelgeschwindigkeit 
konstant sei. Damit wird 


Mm; My 


= r0?cos®’p. 
 m,+tm, , 


Die Reibungsarbeit A,. die während einer Um 
drehung geleistet wird. kann man berechnen. 
wenn man die Radialkraft kennt und Coulomb- 
sche Reibung voraussetzt. Die Reibungsarbeit 
wird, wenn angenommen sei, daß in der Entfer- 
nung r die gleitenden Teile sich berühren und / 
der zugehörige Reibungskoeffizient ist 


In 
fr? @? [eos: pda 


0 


ar. mM, 
RTTute, 


m 
= — ' __f?o®!n. 
m, 
M, 
Ist die zugeführte Energie konstant, so ist der 
Wert der Reibungsarbeit je Umdrehung unver- 
änderlich. Damit wird 


R 
ww“ 


I+ m, 


M; 
d. h. eine vergrößerte Fundamentmasse bedingt 
eine verkleinerte Umdrehungszahl. 


Wien. Fritz Söchting. 282 


= Con8t, 


Einfache Bestimmung des Viereck- 
schwerpunktes. Hierfür gab kürzlich R. Berger 
in der Z. VDI. Bd. 85 (1941) S. 258 eine Konstruktion 
ohne Benutzung der Schwerpunkte der Dreiecke, 
in die das Viereck durch Ziehen je einer Viereck- 
diagonale zerlegt wird. Eine ähnliche und ebenso 
einfache Konstruktion!) ist die folgende, die ich 
schon seit mehr als 2 Jahrzehnten in meinen Vor- 
lesungen mitteile: 

Man trägt (Bild 1) auf der Viereckdiagonalen 
AC von 4 aus die Strecke MC bis C’ auf, ebenso 
auf der Diagonalen BD von B aus die Strecke 
MD bis D’. Sind H, und HA, die Mitten der beiden 
Viereckdiagonalen, so liegt der gesuchte Schwer- 
punkt S im Schnitte der Geraden H,D’ und H,(”. 


1) Die gleiche Konstruktion haben, wie R. Berger in 
einem während der Drucklegung erschienenen Nachtrage 
zur obenangeführten Notiz angibt (Z. VDI Bd. 85 (1941) 
S, 450), die Herren G. Geissenhöner und R,. Hennig 
unabhängig voneinander gefunden. 
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Bild 1. 








Pr B 


Bild 2. Zum Beweise der Konstruktion in Bild 1. 


Beweis: Sind S, und S, die Schwerpunkte der 
Teildreiecke ABD und BCD (Bild 2), F, und Fa 
deren Flächeninhalte, so ist 5, S, | AC und es muß 
der Schwerpunkt S die Strecke S, 5, im Verhält- 
nisse 


| 


RE 


» 


55 HF, 





teilen. Da 


so folgt 








S,S__ML 

Ss, AM 
oder wegen MC — AO. AM =(C' 

a 


Diese Beziehung ist erfüllt, wenn der Schwer- 
punkt S auf der Geraden C’H, liegt; ebenso ist 
zu beweisen, daß die Gerade D’H, eine Schwer- 
linie ist. Daraus folgt aber, daß S im Sehnitt- 
punkt von €’ H, und D’ H, liegen muß. 


(iraz. K. Federhofer. 278 
Sechs Berichtigungen 

zu X M Legrendes Tafeln der Elliptischen 
Integrale (Photographischer Nachdruck 1931 bei 
Wittwer in Stuttgart), aufgefunden und mitgeteilt 
von Prof, Dr. Gustav Witt, Berlin: 

F (9°) = 0,34 923 statt 0,34 293 
S. 364 (29), g = 9°, F (23°) — 1,63 651 statt 1,63 631 
.368 (33), = 85°, F (27%) — 1,56 480 statt 1,56 840 
.380 (45), @ = 80°, F (44°) = 1,59 806 statt 1,58 906 
.396 (61), 9 = 79°, F (64°) — 1.84 793 statt 1.84 693 
S. 416 (81), 9 = 49°, F (87°) = 0,98 328 statt 0,98 238 


8. 350 (15), 7 = 2°, 


A N 


MR 


Stuttgart. Fritz Emde. 277 


BUCHBESPRECHUNGEN 


Hydromechanische Probleme des Schiffsantriebs, 
Teil II, Veröffentlichung der Vorträge, die anläß- 
lich des 25jährigen Bestehens der Hamburgischen 
Schiffbau-Versuchsanstalt am 14. Juni 1939 xe- 
halten wurden, herausgegeben von Dr.-Ing. G. 
Kempf. 236 S. m. 102 Abb. München 1940, 
Verlag R. Oldenbourg. Preis geb. 12,50 M. 

Die in diesem Buch enthaltenen 10 Vorträge ent- 
stammen sämtlich dem Bereich der Hamburgischen 
Schiffbauversuchsanstalt (HSVA.). Sie geben ein 
Bild der Vielseitigkeit der Aufgaben, mit denen 
sich das Institut befaßt und zeigen, daß trotz der 
Überlastung mit laufenden Arbeiten der wissen- 
schaftliche Geist im Schiffbauversuchswesen er- 
freuliche Fortschritte macht. 

Lerbs berichtet über Erfahrungen mit der be- 
kannten Kavitationsanlage, über Neuerungen in 
(der Apparatur und über systematische Versuche, 
die sich auf Einfluß der Kavitationszahl, der Rey- 
nolds- und Froudeschen Zahl und des Luft- 
eehalts des Wassers beziehen. Lerbs und Bau- 
mann teilen die Ergebnisse von Modellmessungen 
über Trägheitsmoment und Dämpfung belasteter 
Schiffsschrauben mit und vergleichen diese mit 
theoretischen Überlegungen: die Schlußfolgerungen 
sind nicht eindeutig. weisen aber darauf hin, daß 
die Drehschwingungen des Systems Antriebs- 
maschine, Wellenleitung. Schraube unter üblichen 
Verhältnissen ausreichend zedämpft sind. Auf den 
Vortrag Kempf „Ergebnisse naturgroßer Schrau- 
benversuche“ ... seien alle nachdrücklich hinge- 
wiesen, die sich für das Problem der Übertragung 
von Modellergebnissen auf schraubengetrienene 


Großausführungen interessieren. Durch eine geist- 
volle Versuchsanordnung war es möglich, die zahl- 
reichen Einflüsse zu trennen und zu konkreten 
Folgerungen zu gelangen. Beachtenswert eine ein- 
fache Formel für den Einfluß der Öberflächen- 
reibung auf den Wirkungsgrad des Propellers. - 
Sehmierschalski bringt die im Tandem- 
propeller liegenden Möglichkeiten in Erinnerung 
Vorzüge hinsichtlich Kavitation!); Heckscher 
veröffentlicht umfangreiches statistisch durchge- 
arbeitetes Material über Schiffsform und Wider- 
stand; Graff behandelt das etwas vernach 
lässigte Gebiet (des Manövrierens von Schiffen 
(Kursstetigkeit. Bewegung auf gekrümmter Balın. 
Stabilität im Drehkreis). Die Arbeit von Helm 
„Tiefen- und Breiteneinfluß von Kanälen auf den 
Sehiffswiderstand“ bedeutet einen aufschlußreichen 
Beitrag zur Versuchstechnik auf diesem Gebiet: 
verschiedene Unstimmigkeiten werden geklärt und 
Material über den Einfluß der Breite auf den 
Widerstand beigebracht, das von anderer Seite ge- 
vebene theoretische Überlegungen bestätigt. 
Zwei Abhandlungen über «den Schiffsbewuchs 
Kühl) und Eintluß des Bewuchses auf die Ge- 
sehwindigkeit (Sehubart) beschließen das Sam- 
melwerk. 
Dessau-Roßlau. Weinblum. 272 
Bauelemente des Flurzeugs, nach 
Vorlesungen von Dr.-Ing. HERBERT WAGNER, 
ehem. Professor und Leiter des Flugtechnischen 


instituts an der Technischen Hochschule Berlin, 
bearbeitet von Dipl.-Ing. Gotthold Kimm. 
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296 S. m. 280 Bildern. München und Berlin 1940, ix 
Verlag R. Oldenbourg. Preis geb. 12 M. Fi= - [erıw Re 
. = ı a i 2ni, 
Kennzeichnend für das vorliegende Buch ist das Ei 


Bestreben, die im Tragsystem der Flugzeuge auf- 
tretenden Beanspruchungen mit einfachsten Mitteln 
und mit knapp gehaltenem Text darzustellen. Da- 
bei tritt als bemerkenswert der Umstand zutage, 
daß bei dieser Art der Behandlung ihrer „Bau- 
elemente“ so ziemlich die ganze Festigkeitslehre 
vorgetragen wird. Nach Erklärung der grund- 
lergenden Verfahren der Statik für die aus Stä- 
ben und Scheiben zusammengesetzten Tragsysteme 
und der einfachen Beanspruchungen ist besonders 
hervorzuheben die Berechnung «der gekrümmten 
Träger, der Rahmen, Hohlkörper und Schalen, die 
im Flugzeugbau vielfache Verwendung finden: 
ferner die Betrachtung der Schubverbände von 
Bierungsträgern in verschiedenen Ausführungen. 
die Lehre vom Schubmittelpunkt, die „Einleitung“ 
von Schubkräften durch besonders geformte Bau- 
elieder, schließlich die zusammengesetzten Bean- 
spruchungen, wie Biegungsverdrehung, Knickung 
und Knickbiegung u. dgl. — Ein wesentliches Ge- 
heimnis für den unerreicht hohen Stand des deut- 
schen Flugzeugbaues liegt in der tiefgehenden 
theoretischen Durchdringung seiner Baugedanken. 
Daß die statische Beherrschung des Tragsystems 
mit einfachen Mitteln auch weiteren Kreisen ver- 
inittelt werden kann, dafür zeugt das vorliegende 
Buch, das bei aller Beschränkung der verwendeten 
mathematischen Hilfsmittel die wesentlichen Merk- 
male der statischen Methoden und ihre weit- 
reichende Anwendbarkeit in vorbildlicher Weise 
darlegt. 


Karlsruhe. Th. Pösehl. 279 


Professor Dr. phil. Dr.-Ing. e. h. KARL WILLY 
WAGNER, Operatorenreehnung nebst 
Anwendungen in Physik und Tech- 
nik, nach Vorträgen, veranstaltet durch den Be- 
zirk Berlin des Verbandes Deutscher Elektrotechni- 
ker — vormals Elektrotechnischer Verein e. V. — 
in Gemeinschaft mit dem Außeninstitut der Tech- 
nischen Hochschule Berlin. XIV + 448 S. m. 126 
Abb. Leipzig 1940, Verlag Johann Ambrosius Barth. 
Preis brosch. 27.60 M. 


Der Verf. des vorliegenden Werkes hat sich be- 
kanntlich schon früher um die Popularisierung des 
Heaviside-Kalküls in Deutschland Verdienste er- 
worben. Dieser kühne Kalkül, der auf der Ablösung 
des Differentiationssymbols von der Funktion und 
der Behandlung als multiplikative Größe beruht 
und auf Grund gewagter Deutungen der dabei ent- 
stehenden Ausdrücke manchmal zu überraschenil 
einfacher Lösung von Differentialgleichungen führt, 
stellt ein im mathematischen Bezirk sonst ganz un- 
ewohntes Experimentieren dar, über dessen wirk- 
lichen Erfolg niemals Sicherheit besteht. Dem Be- 
‚dürfnis nach einer „Begründung“, d. h. nach An- 
abe einer exakt mathematischen Theorie, die ge- 
wissermaßen zum Heaviside-Kalkül isomorph ist. 
also dasselbe mit anderen Worten sagt, dabei aber 
renau die Grenzen des Erlaubten aufzuweisen ge- 
stattet, hat Wagner schon im Jahre 1916 durch eine 
sehr bekannt gewordene Arbeit zu genügen ver- 
sucht. Da die Exponentialfunktion das im Hea- 
viside-Kalkül benutzte Postulat, daß die Differen- 
tiation mit einer algebraischen Multiplikation äqui- 


F 2a "3 d . i 
valent sein soll. tatsächlich erfüllt: e0% — g.e** 
dir 


lage es nahe, eine beliebige Funktion als Super- 
position von Exponentialfunktionen aufzufassen, 
was auf die „Spektraldarstellung‘“ durch ein Fourier- 
Integral hinauskommt: 


In dieser Darstellung äußert sich die Differen- 
tiation von Ft), d. h. die symbolische Multiplika- 
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tion mit in einer gewöhnlichen Multiplikation der 
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zugeordneten Amplitudenfunktion f(s) mit der al- 
eebraischen Größe s: 
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Damit wird sogleich klar. daß die Unmöglichkeit, 
den Heaviside-Kalkül aus sich heraus zu recht- 
fertigen, «darin begründet ist, daß (dieser bei der 
Umdeutung des Differentiationssymbols als mult;- 
plikative Größe an der Funktion F kleben bleibt, 
während eine effektive Multiplikation im üblichen 
mathematischen Sinn gar nieht an F, sondern nur 
an einem „Abbild”, einer zugeordneten Funktion / 
vorliegt. (Man muß also neben dem Bereich der 
„Originalfunktionen“ F noch einen Bereich von 
„Bildfunktionen“ f betrachten.) In der Folge hat 
sich nun gezeigt, daß eine wirklich völlig befriedi- 
ende Behandlung des Heaviside-Kalküls nur mög- 
lich ist, wenn man gleich von Anfang an von F zu / 
übergeht. Dazu braucht man eine Darstellung von / 
durch F (die Umkehrung von (1)), die sich zwang 
!os aus der Darstellung von F durch ein Fonrier- 
sches Doppelintegral ergibt und die Form hat: 


f(s) [« de 1 EEE 3 


Die rechte Seite stellt die bekannte Laplace-Trans- 
formation dar, wobei die Integrationsgrenzen sich 
nach dem der Betrachtung von F zugrunde liegen- 
dien Intervall richten, also z. B. bei einem „Dauer- 
zustand“ — oo und oo lauten oder 0 und oo bei 
einem im Zeitpunkt ?=0 einsetzenden .„Einschalt- 
vorgang“. 


Der Verf. hat sich nun auch vorgenommen, seine 
für Physiker und Techniker bestimmte Darstellunz 
der Operatorenrechnung auf die bewährte Theorie 
der Laplace-Transformation zu gründen, wie das in 
dem mehr für rein mathematische Kreise geschrie- 
benen Buch des Referenten „Theorie und Anwen- 
dung der Laplace-Transformation“ (Berlin 1937) ge- 
schehen ist. Es ist zu bedauern, daß davon in zwei 
Richtungen abgewichen wurde: Einmal benutzt 
Wagner. um in“engere formale Übereinstimmung 
mit Heaviside zu kominen, statt (2) die mit s multi- 
plizierte Laplace-Transformation, was zu mancher- 
lei Umständlichkeiten führt. Zum anderen Knüpft 
er aus alter Vorliebe für seine Methode aus dem 
Jahre 1916 häufig statt an die Laplace-Transfor- 
mation an das komplexe Integral (1) an, was ihn 
an manchen Stellen an einer befriedigenden Durch- 
leuchtung der Materie hindert. Das komplexe Inte- 
gral hat auch nach Ansicht des Referenten durch- 
aus seine ihm in der Gesamttheorie zukommende 
Stellung einzunehmen, diese liegt aber nicht am 
Anfang. sondern am Ende der Prozedur, nämlich 
wenn es sich darum handelt, von der auf symbo- 
lischenn Wege, d. h. im .„Bildbereich“ heraus- 
xerechneten Funktion f(s) zu der eigentlichen 
Öriginalfunktion F(t) zurückzugelangen. Hier bie- 
tet das komplexe Integral manchmal, wie eine ge- 
nauere Analyse zeigt: nicht immer, die Möglichkeit. 
aus / die Funktion F zu berechnen. evtl. unter 
Heranziehung funktionentheoretischer Methoden 
Residuenrechnung), wobei zu bemerken ist, daß 
die meisten, die das komplexe Integral benutzen, 
sich hinsichtlich der Bedingungen, unter denen es 
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tatsächlich verwendet werden kann, allzu großer 
Sorglosigkeit hingeben. Stellt man das komplexe 
Integral aber an den Anfang, d. h. setzt man die 
gesuchte Funktion von vornherein in der Form (1) 
an, so stößt man vor allem auf die Schwierigkeit. 
der Funktion die richtigen „Anfangsbedingungen“ 
aufzuprägen, was im Falle der primären Anwen- 
dung der Laplace-Transformation (2) auf die ge- 
gebene Differentialgleichung gerade besonders 
leicht ist. Alles was in diesem Punkt zugunsten 
des komplexen Integrals vorgebracht worden ist. 
2. B. von Bromwich in einer oft zitierten, mit der 
Wagnerschen fast gleichzeitigen Arbeit, hält einer 
mathematischen Kritik nicht stand. Ein schlagen- 
des Beispiel für die Überlegenheit der primären 
Anwendung der Laplace-Transformation über den 
Ansatz des komplexen Integrals ist die Ableitung 
des Heavisideschen „Expansion Theorem“ für ge- 
wöhnliche Differentialgleichungen. Für die Laplace- 
Transformation ist der Beweis dieses Satzes fast 
eine Trivialität, die nur die elementarsten Eigen- 
schaften der Laplace-Transformation benutzt und 
ganz im Reellen verläuft, während die Ableitung 
auf dem Weg über das komplexe Integral, wie sie 
Wagener in Anlehnung an seine frühere Arbeit vor- 
zieht, sehr umständlich ist, wenn man sie exakt 
durehführen will. und allerlei Hilfsmittel aus der 
Funktionentheorie benutzen muß. (Die Ausdehnung 
des Expansion Theorem auf partielle Differential- 
eleichungen, wo die bei gewöhnlichen Differential- 
eleichungen auftretende gebrochen rationale Funk- 
tion durch eine meromorphe Funktion zu ersetzen 
ist, bleibt übrigens bei Wagner genau so unbefrie- 
digend wie bei den meisten anderen Autoren, da 
die Schwierigkeiten dieses Falles bedeutend unter- 
schätzt werden.) 

Nach diesen Ausstellungen ist es an der Zeit, 
nun auch die Vorzüge des Buches gebührend her- 
vorzuheben. Diese liegen voı allem in der bewun- 
dernswerten Mannigfaltiekeit der Beispiele und der 
Art ihrer Einzelbehandlung, woran derjenige, dem 
eine strenge systematische allgemeine Anleitung 
nieht liegt und der lieber an konkreten Einzel- 
fällen nach und nach alle Besonderheiten der Me- 
thode lernen will (das ist ja der den meisten Tech- 
nikern sympathischste Weg), das für die Anwen- 
dungen Wichtige sich aneignen kann. Wenn das 
Buch nach der mathematischen Seite noch nicht 
allen Anforderungen, die man billigerweise auch an 
ein für Techniker bestimmtes Werk stellen muß, 
renügt. so ist es doch als eine sehr wertvolle 
Etappe auf dem Wege von der früheren sorglosen 
Jonglierkunst «des Heaviside-Kalküls zu einer die 
modernen Fortschritte berücksichtigenden Darstel- 
lung, die die Schwierigkeiten nicht geflissentlich 
übersieht, sondern ihnen zu Leibe geht und den 
Ingenieur zu verantwortungsbewußtem, aber nicht 
zimperlichem Handeln erzieht, sehr zu begrüßen. 


Freiburg i. B. {z. Zt. bei der Luftwaffe). 
Doetsch. 270 


Dr. HUGO SIRK, Dozent a. d. Universität Wien, 
ehemal. Prof. a. d. Universität Laibach, Mathe- 
matik für Naturwissenschaftler und 
Chemiker, eine Einführung in die 
Anwendungen der höheren Mathema- 


tik. XI+268 S. m. 126 Abb. und einer Aus- 
schlagtafel. Dresden und Leipzig 1941. Verlag 


Theodor Steinkopff. Preis geb. 12 M. 


Wie der Verfasser im Vorwort schreibt, ist der 
leitende Grundgedanke seines Buches die Erarbei- 
tung der mathematischen Begriffe an Beispielen 
aus den Naturwissenschaften. Das klingt für die 


jungen Chemiker und Naturwissenschaftler, an die 
sich der Verfasser in erster Linie wendet, gewiß 
sehr verführerisch; sie werden es begrüßen, daß sie 
ihr mathematisches Rüstzeug nur an Beispielen 
ihres engeren Interessengebietes kennenlernen 
sollen und sich nicht um die sog. graue Theorie 
zu kümmern brauchen, daß sie sozusagen im Son- 
nenbad praktischer Beispiele mathematisch schwim- 
men lernen werden. Nun soll auch zugegeben sein, 
daß der Verfasser eine große Zahl lehrreicher Bei- 
spiele aus den verschiedensten Anwendungs- 
gebieten der Mathematik zusammengetragen und 
geschickt dargestellt hat. Trotzdem scheint es 
sehr fraglich. ob das Buch geeignet ist, auch dem 
mathematisch weniger Vorgebildeten ein klares und 
verständliches Bild selbst von den schwierigsten 
Gebieten der Mathematik zu vermitteln, wie es auf 
dem Umschlag heißt. Niemand wird von einer 
solchen Einführung erwarten, daß sie den Leser in 
alle Feinheiten der Infinitesimalrechnung einweiht: 
aber was der Leser von der Differential- und In- 
tegralrechnung erfährt, das sollte doch unbedingt 
einwandfrei und zuverlässig sein! Dieser Forde- 
rung genügt nun das Sirksche Buch leider kaum. 
Denn dort, wo der Verfasser ins Freischwimmer- 
bassin mathematischer Theorie gerät, wo er also 
den begrifflichen Schwierigkeiten nicht aus dem 
Wege geht, zeigt das Buch beträchtliche Schön- 
heitsfehler (f(x) >0 ist nieht notwendig dafür, 
daß /(z) zunimmt: es gibt doch ein F (x). das der 
pe 
x”? 
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die Konvergenzprobe “ '<a<1 für die Reihe 
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e} a, gilt in dieser Form nur für Reihen mit 
positiven Gliedern und ist auch keineswegs 
notwendig für die Konvergenz; z2=f (r,y) 


Dy/ Of 
d« —+ dy dy. 


dz E- 
wenn sich 2 um d@,y um dy ändert: aus a<b 
folgt nieht immer am <bım, sondern nur dann, 
wenn m positiv ist!). Derartige Ungereimtheiten 
können einem jungen Studenten der Chemie, der 
seinem mathematischen Lehrbuch sicherlich nicht 
sehr kritisch gegenüberstehen wird, leicht gefähr- 
lich werden. Wer dagegen schon zwischen mathe- 
matisch Gut und Böse zu unterscheiden weiß, der 
mag aus den abwechslungsreichen Beispielen des 
Sirk schen Buches mancherlei Anregung schöpfen. 


ändert sich nieht um 


Was den Inhalt betrifft, der von den Anfängen 
der Differentialrechnung bis zu den einfachsten ge- 
wöhnlichen und partiellen Differentialgleichungen 
reicht, so wäre noch eine Darstellung der zeich- 
nerischen Integration erwünscht. 
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VDI-Jahrbuch 194. 
1941. VDI-Verlag G. m.b.H. 
für VDI-Mitglieder 3,15 M. 


Mit dem Jahrbuch 1940 liegt der siebente Band 
der VDI-Jahrbücher vor. Durch die Zeitumstände 
bedingt, erscheint es mit etwa einjähriger Ver- 
spätung. (Gegenüber den bisherigen Jahrgängen 
weist der vorliegende als wesentliche Verbesse- 
rung die auf, daß bei den rund 13 000 Schrifttums- 
nachweisen die Verfassernamen angegeben sind. 
Eine besondere Empfehlung des Jahrbuches, das 
heute wohl überall den weiterarbeitenden In- 
genieuren eine erwünschte Hilfe ist, erübrigt sich. 
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Preis brosch. 
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